[360] | 1 | /* -*- mode: C++; indent-tabs-mode: nil; -*- |
---|
| 2 | * |
---|
| 3 | * This file is a part of LEMON, a generic C++ optimization library. |
---|
| 4 | * |
---|
| 5 | * Copyright (C) 2003-2008 |
---|
| 6 | * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
---|
| 7 | * (Egervary Research Group on Combinatorial Optimization, EGRES). |
---|
| 8 | * |
---|
| 9 | * Permission to use, modify and distribute this software is granted |
---|
| 10 | * provided that this copyright notice appears in all copies. For |
---|
| 11 | * precise terms see the accompanying LICENSE file. |
---|
| 12 | * |
---|
| 13 | * This software is provided "AS IS" with no warranty of any kind, |
---|
| 14 | * express or implied, and with no claim as to its suitability for any |
---|
| 15 | * purpose. |
---|
| 16 | * |
---|
| 17 | */ |
---|
| 18 | |
---|
| 19 | #ifndef LEMON_NAUTY_READER_H |
---|
| 20 | #define LEMON_NAUTY_READER_H |
---|
| 21 | |
---|
| 22 | #include <vector> |
---|
| 23 | #include <iostream> |
---|
| 24 | #include <string> |
---|
| 25 | |
---|
| 26 | /// \ingroup io_group |
---|
| 27 | /// |
---|
| 28 | /// @defgroup nauty_group NAUTY format |
---|
| 29 | /// |
---|
| 30 | /// \brief Read \e Nauty format |
---|
| 31 | /// |
---|
| 32 | /// Tool to read graphs from \e Nauty format data |
---|
| 33 | |
---|
| 34 | /// \ingroup nauty_group |
---|
| 35 | /// \file |
---|
| 36 | /// \brief Nauty file reader. |
---|
| 37 | namespace lemon { |
---|
| 38 | |
---|
| 39 | /// \ingroup nauty_group |
---|
| 40 | /// |
---|
| 41 | /// \brief Nauty file reader |
---|
| 42 | /// |
---|
| 43 | /// The \e geng program is in the \e gtools suite of the nauty |
---|
| 44 | /// package. This tool can generate all non-isomorphic undirected |
---|
| 45 | /// graphs with given node number from several classes (for example, |
---|
| 46 | /// general, connected, biconnected, triangle-free, 4-cycle-free, |
---|
| 47 | /// bipartite and graphs with given edge number and degree |
---|
| 48 | /// constraints). This function reads a \e nauty \e graph6 \e format |
---|
| 49 | /// line from the given stream and builds it in the given graph. |
---|
| 50 | /// |
---|
| 51 | /// The site of nauty package: http://cs.anu.edu.au/~bdm/nauty/ |
---|
| 52 | /// |
---|
| 53 | /// For example, the number of all non-isomorphic connected graphs |
---|
| 54 | /// can be computed with following code. |
---|
| 55 | ///\code |
---|
| 56 | /// int num = 0; |
---|
| 57 | /// SmartGraph graph; |
---|
| 58 | /// while(readNauty(std::cin, graph)) { |
---|
| 59 | /// PlanarityChecking<SmartUGraph> pc(graph); |
---|
| 60 | /// if (pc.run()) ++num; |
---|
| 61 | /// } |
---|
| 62 | /// std::cout << "Number of planar graphs: " << num << std::endl; |
---|
| 63 | ///\endcode |
---|
| 64 | /// |
---|
| 65 | /// The nauty files are quite huge, therefore instead of the direct |
---|
| 66 | /// file generation the pipelining is recommended. |
---|
| 67 | ///\code |
---|
| 68 | /// ./geng -c 10 | ./num_of_pg |
---|
| 69 | ///\endcode |
---|
| 70 | template <typename Graph> |
---|
| 71 | std::istream& readNauty(std::istream& is, Graph& graph) { |
---|
| 72 | graph.clear(); |
---|
| 73 | |
---|
| 74 | std::string line; |
---|
| 75 | if (getline(is, line)) { |
---|
| 76 | int index = 0; |
---|
| 77 | |
---|
| 78 | int n; |
---|
| 79 | |
---|
| 80 | if (line[index] == '>') { |
---|
| 81 | index += 10; |
---|
| 82 | } |
---|
| 83 | |
---|
| 84 | char c = line[index++]; c -= 63; |
---|
| 85 | if (c != 63) { |
---|
| 86 | n = int(c); |
---|
| 87 | } else { |
---|
| 88 | c = line[index++]; c -= 63; |
---|
| 89 | n = (int(c) << 12); |
---|
| 90 | c = line[index++]; c -= 63; |
---|
| 91 | n |= (int(c) << 6); |
---|
| 92 | c = line[index++]; c -= 63; |
---|
| 93 | n |= int(c); |
---|
| 94 | } |
---|
| 95 | |
---|
| 96 | std::vector<typename Graph::Node> nodes; |
---|
| 97 | for (int i = 0; i < n; ++i) { |
---|
| 98 | nodes.push_back(graph.addNode()); |
---|
| 99 | } |
---|
| 100 | |
---|
| 101 | int bit = -1; |
---|
| 102 | for (int j = 0; j < n; ++j) { |
---|
| 103 | for (int i = 0; i < j; ++i) { |
---|
| 104 | if (bit == -1) { |
---|
| 105 | c = line[index++]; c -= 63; |
---|
| 106 | bit = 5; |
---|
| 107 | } |
---|
| 108 | bool b = (c & (1 << (bit--))) != 0; |
---|
| 109 | |
---|
| 110 | if (b) { |
---|
| 111 | graph.addEdge(nodes[i], nodes[j]); |
---|
| 112 | } |
---|
| 113 | } |
---|
| 114 | } |
---|
| 115 | } |
---|
| 116 | return is; |
---|
| 117 | } |
---|
| 118 | } |
---|
| 119 | |
---|
| 120 | #endif |
---|