/* -*- mode: C++; indent-tabs-mode: nil; -*-

  *

  * This file is a part of LEMON, a generic C++ optimization library.

  *

  * Copyright (C) 2003-2008

  * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport

  * (Egervary Research Group on Combinatorial Optimization, EGRES).

  *

  * Permission to use, modify and distribute this software is granted

  * provided that this copyright notice appears in all copies. For

  * precise terms see the accompanying LICENSE file.

  *

  * This software is provided "AS IS" with no warranty of any kind,

  * express or implied, and with no claim as to its suitability for any

  * purpose.

  *

  */

 ///\ingroup lemon_io

 ///\file

 ///\brief \ref lgf-format "LEMON Graph Format" reader.

 #ifndef LEMON_LGF_READER_H

 #define LEMON_LGF_READER_H

 #include <iostream>

 #include <fstream>

 #include <sstream>

 #include <set>

 #include <map>

 #include <lemon/core.h>

 #include <lemon/lgf_writer.h>

 #include <lemon/concept_check.h>

 #include <lemon/concepts/maps.h>

 namespace lemon {

   namespace _reader_bits {

     template <typename Value>

     struct DefaultConverter {

       Value operator()(const std::string& str) {

         std::istringstream is(str);

         Value value;

         if (!(is >> value)) {

           throw FormatError("Cannot read token");

         }

         char c;

         if (is >> std::ws >> c) {

           throw FormatError("Remaining characters in token");

         }

         return value;

       }

     };

     template <>

     struct DefaultConverter<std::string> {

       std::string operator()(const std::string& str) {

         return str;

       }

     };

     template <typename _Item>

     class MapStorageBase {

     public:

       typedef _Item Item;

     public:

       MapStorageBase() {}

       virtual ~MapStorageBase() {}

       virtual void set(const Item& item, const std::string& value) = 0;

     };

     template <typename _Item, typename _Map,

               typename _Converter = DefaultConverter<typename _Map::Value> >

     class MapStorage : public MapStorageBase<_Item> {

     public:

       typedef _Map Map;

       typedef _Converter Converter;

       typedef _Item Item;

     private:

       Map& _map;

       Converter _converter;

     public:

       MapStorage(Map& map, const Converter& converter = Converter())

         : _map(map), _converter(converter) {}

       virtual ~MapStorage() {}

       virtual void set(const Item& item ,const std::string& value) {

         _map.set(item, _converter(value));

       }

     };

     template <typename _Graph, bool _dir, typename _Map,

               typename _Converter = DefaultConverter<typename _Map::Value> >

     class GraphArcMapStorage : public MapStorageBase<typename _Graph::Edge> {

     public:

       typedef _Map Map;

       typedef _Converter Converter;

       typedef _Graph Graph;

       typedef typename Graph::Edge Item;

       static const bool dir = _dir;

     private:

       const Graph& _graph;

       Map& _map;

       Converter _converter;

     public:

       GraphArcMapStorage(const Graph& graph, Map& map,

                          const Converter& converter = Converter())

         : _graph(graph), _map(map), _converter(converter) {}

       virtual ~GraphArcMapStorage() {}

       virtual void set(const Item& item ,const std::string& value) {

         _map.set(_graph.direct(item, dir), _converter(value));

       }

     };

     class ValueStorageBase {

     public:

       ValueStorageBase() {}

       virtual ~ValueStorageBase() {}

       virtual void set(const std::string&) = 0;

     };

     template <typename _Value, typename _Converter = DefaultConverter<_Value> >

     class ValueStorage : public ValueStorageBase {

     public:

       typedef _Value Value;

       typedef _Converter Converter;

     private:

       Value& _value;

       Converter _converter;

     public:

       ValueStorage(Value& value, const Converter& converter = Converter())

         : _value(value), _converter(converter) {}

       virtual void set(const std::string& value) {

         _value = _converter(value);

       }

     };

     template <typename Value>

     struct MapLookUpConverter {

       const std::map<std::string, Value>& _map;

       MapLookUpConverter(const std::map<std::string, Value>& map)

         : _map(map) {}

       Value operator()(const std::string& str) {

         typename std::map<std::string, Value>::const_iterator it =

           _map.find(str);

         if (it == _map.end()) {

           std::ostringstream msg;

           msg << "Item not found: " << str;

           throw FormatError(msg.str());

         }

         return it->second;

       }

     };

     template <typename Graph>

     struct GraphArcLookUpConverter {

       const Graph& _graph;

       const std::map<std::string, typename Graph::Edge>& _map;

       GraphArcLookUpConverter(const Graph& graph,

                               const std::map<std::string,

                                              typename Graph::Edge>& map)

         : _graph(graph), _map(map) {}

       typename Graph::Arc operator()(const std::string& str) {

         if (str.empty() || (str[0] != '+' && str[0] != '-')) {

           throw FormatError("Item must start with '+' or '-'");

         }

         typename std::map<std::string, typename Graph::Edge>

           ::const_iterator it = _map.find(str.substr(1));

         if (it == _map.end()) {

           throw FormatError("Item not found");

         }

         return _graph.direct(it->second, str[0] == '+');

       }

     };

     inline bool isWhiteSpace(char c) {

       return c == ' ' || c == '\t' || c == '\v' ||

         c == '\n' || c == '\r' || c == '\f';

     }

     inline bool isOct(char c) {

       return '0' <= c && c <='7';

     }

     inline int valueOct(char c) {

       LEMON_ASSERT(isOct(c), "The character is not octal.");

       return c - '0';

     }

     inline bool isHex(char c) {

       return ('0' <= c && c <= '9') ||

         ('a' <= c && c <= 'z') ||

         ('A' <= c && c <= 'Z');

     }

     inline int valueHex(char c) {

       LEMON_ASSERT(isHex(c), "The character is not hexadecimal.");

       if ('0' <= c && c <= '9') return c - '0';

       if ('a' <= c && c <= 'z') return c - 'a' + 10;

       return c - 'A' + 10;

     }

     inline bool isIdentifierFirstChar(char c) {

       return ('a' <= c && c <= 'z') ||

         ('A' <= c && c <= 'Z') || c == '_';

     }

     inline bool isIdentifierChar(char c) {

       return isIdentifierFirstChar(c) ||

         ('0' <= c && c <= '9');

     }

     inline char readEscape(std::istream& is) {

       char c;

       if (!is.get(c))

         throw FormatError("Escape format error");

       switch (c) {

       case '\\':

         return '\\';

       case '\"':

         return '\"';

       case '\'':

         return '\'';

       case '\?':

         return '\?';

       case 'a':

         return '\a';

       case 'b':

         return '\b';

       case 'f':

         return '\f';

       case 'n':

         return '\n';

       case 'r':

         return '\r';

       case 't':

         return '\t';

       case 'v':

         return '\v';

       case 'x':

         {

           int code;

           if (!is.get(c) || !isHex(c))

             throw FormatError("Escape format error");

           else if (code = valueHex(c), !is.get(c) || !isHex(c)) is.putback(c);

           else code = code * 16 + valueHex(c);

           return code;

         }

       default:

         {

           int code;

           if (!isOct(c))

             throw FormatError("Escape format error");

           else if (code = valueOct(c), !is.get(c) || !isOct(c))

             is.putback(c);

           else if (code = code * 8 + valueOct(c), !is.get(c) || !isOct(c))

             is.putback(c);

           else code = code * 8 + valueOct(c);

           return code;

         }

       }

     }

     inline std::istream& readToken(std::istream& is, std::string& str) {

       std::ostringstream os;

       char c;

       is >> std::ws;

       if (!is.get(c))

         return is;

       if (c == '\"') {

         while (is.get(c) && c != '\"') {

           if (c == '\\')

             c = readEscape(is);

           os << c;

         }

         if (!is)

           throw FormatError("Quoted format error");

       } else {

         is.putback(c);

         while (is.get(c) && !isWhiteSpace(c)) {

           if (c == '\\')

             c = readEscape(is);

           os << c;

         }

         if (!is) {

           is.clear();

         } else {

           is.putback(c);

         }

       }

       str = os.str();

       return is;

     }

     class Section {

     public:

       virtual ~Section() {}

       virtual void process(std::istream& is, int& line_num) = 0;

     };

     template <typename Functor>

     class LineSection : public Section {

     private:

       Functor _functor;

     public:

       LineSection(const Functor& functor) : _functor(functor) {}

       virtual ~LineSection() {}

       virtual void process(std::istream& is, int& line_num) {

         char c;

         std::string line;

         while (is.get(c) && c != '@') {

           if (c == '\n') {

             ++line_num;

           } else if (c == '#') {

             getline(is, line);

             ++line_num;

           } else if (!isWhiteSpace(c)) {

             is.putback(c);

             getline(is, line);

             _functor(line);

             ++line_num;

           }

         }

         if (is) is.putback(c);

         else if (is.eof()) is.clear();

       }

     };

     template <typename Functor>

     class StreamSection : public Section {

     private:

       Functor _functor;

     public:

       StreamSection(const Functor& functor) : _functor(functor) {}

       virtual ~StreamSection() {}

       virtual void process(std::istream& is, int& line_num) {

         _functor(is, line_num);

         char c;

         std::string line;

         while (is.get(c) && c != '@') {

           if (c == '\n') {

             ++line_num;

           } else if (!isWhiteSpace(c)) {

             getline(is, line);

             ++line_num;

           }

         }

         if (is) is.putback(c);

         else if (is.eof()) is.clear();

       }

     };

   }

   template <typename Digraph>

   class DigraphReader;

   template <typename Digraph>

   DigraphReader<Digraph> digraphReader(Digraph& digraph, 

                                        std::istream& is = std::cin);

   template <typename Digraph>

   DigraphReader<Digraph> digraphReader(Digraph& digraph, const std::string& fn);

   template <typename Digraph>

   DigraphReader<Digraph> digraphReader(Digraph& digraph, const char *fn);

   /// \ingroup lemon_io

   ///

   /// \brief \ref lgf-format "LGF" reader for directed graphs

   ///

   /// This utility reads an \ref lgf-format "LGF" file.

   ///

   /// The reading method does a batch processing. The user creates a

   /// reader object, then various reading rules can be added to the

   /// reader, and eventually the reading is executed with the \c run()

   /// member function. A map reading rule can be added to the reader

   /// with the \c nodeMap() or \c arcMap() members. An optional

   /// converter parameter can also be added as a standard functor

   /// converting from \c std::string to the value type of the map. If it

   /// is set, it will determine how the tokens in the file should be

   /// converted to the value type of the map. If the functor is not set,

   /// then a default conversion will be used. One map can be read into

   /// multiple map objects at the same time. The \c attribute(), \c

   /// node() and \c arc() functions are used to add attribute reading

   /// rules.

   ///

   ///\code

   /// DigraphReader<Digraph>(digraph, std::cin).

   ///   nodeMap("coordinates", coord_map).

   ///   arcMap("capacity", cap_map).

   ///   node("source", src).

   ///   node("target", trg).

   ///   attribute("caption", caption).

   ///   run();

   ///\endcode

   ///

   /// By default the reader uses the first section in the file of the

   /// proper type. If a section has an optional name, then it can be

   /// selected for reading by giving an optional name parameter to the

   /// \c nodes(), \c arcs() or \c attributes() functions.

   ///

   /// The \c useNodes() and \c useArcs() functions are used to tell the reader

   /// that the nodes or arcs should not be constructed (added to the

   /// graph) during the reading, but instead the label map of the items

   /// are given as a parameter of these functions. An

   /// application of these functions is multipass reading, which is

   /// important if two \c \@arcs sections must be read from the

   /// file. In this case the first phase would read the node set and one

   /// of the arc sets, while the second phase would read the second arc

   /// set into an \e ArcSet class (\c SmartArcSet or \c ListArcSet).

   /// The previously read label node map should be passed to the \c

   /// useNodes() functions. Another application of multipass reading when

   /// paths are given as a node map or an arc map.

   /// It is impossible to read this in

   /// a single pass, because the arcs are not constructed when the node

   /// maps are read.

   template <typename _Digraph>

   class DigraphReader {

   public:

     typedef _Digraph Digraph;

     TEMPLATE_DIGRAPH_TYPEDEFS(Digraph);

   private:

     std::istream* _is;

     bool local_is;

     std::string _filename;

     Digraph& _digraph;

     std::string _nodes_caption;

     std::string _arcs_caption;

     std::string _attributes_caption;

     typedef std::map<std::string, Node> NodeIndex;

     NodeIndex _node_index;

     typedef std::map<std::string, Arc> ArcIndex;

     ArcIndex _arc_index;

     typedef std::vector<std::pair<std::string,

       _reader_bits::MapStorageBase<Node>*> > NodeMaps;

     NodeMaps _node_maps;

     typedef std::vector<std::pair<std::string,

       _reader_bits::MapStorageBase<Arc>*> >ArcMaps;

     ArcMaps _arc_maps;

     typedef std::multimap<std::string, _reader_bits::ValueStorageBase*>

       Attributes;

     Attributes _attributes;

     bool _use_nodes;

     bool _use_arcs;

     bool _skip_nodes;

     bool _skip_arcs;

     int line_num;

     std::istringstream line;

   public:

     /// \brief Constructor

     ///

     /// Construct a directed graph reader, which reads from the given

     /// input stream.

     DigraphReader(Digraph& digraph, std::istream& is = std::cin)

       : _is(&is), local_is(false), _digraph(digraph),

         _use_nodes(false), _use_arcs(false),

         _skip_nodes(false), _skip_arcs(false) {}

     /// \brief Constructor

     ///

     /// Construct a directed graph reader, which reads from the given

     /// file.

     DigraphReader(Digraph& digraph, const std::string& fn)

       : _is(new std::ifstream(fn.c_str())), local_is(true),

         _filename(fn), _digraph(digraph),

         _use_nodes(false), _use_arcs(false),

         _skip_nodes(false), _skip_arcs(false) {

       if (!(*_is)) {

         delete _is;

         throw IoError("Cannot open file", fn);

       }

     }

     /// \brief Constructor

     ///

     /// Construct a directed graph reader, which reads from the given

     /// file.

     DigraphReader(Digraph& digraph, const char* fn)

       : _is(new std::ifstream(fn)), local_is(true),

         _filename(fn), _digraph(digraph),

         _use_nodes(false), _use_arcs(false),

         _skip_nodes(false), _skip_arcs(false) {

       if (!(*_is)) {

         delete _is;

         throw IoError("Cannot open file", fn);

       }

     }

     /// \brief Destructor

     ~DigraphReader() {

       for (typename NodeMaps::iterator it = _node_maps.begin();

            it != _node_maps.end(); ++it) {

         delete it->second;

       }

       for (typename ArcMaps::iterator it = _arc_maps.begin();

            it != _arc_maps.end(); ++it) {

         delete it->second;

       }

       for (typename Attributes::iterator it = _attributes.begin();

            it != _attributes.end(); ++it) {

         delete it->second;

       }

       if (local_is) {

         delete _is;

       }

     }

   private:

     template <typename DGR>

     friend DigraphReader<DGR> digraphReader(DGR& digraph, std::istream& is);

     template <typename DGR>

     friend DigraphReader<DGR> digraphReader(DGR& digraph, 

                                             const std::string& fn);

     template <typename DGR>

     friend DigraphReader<DGR> digraphReader(DGR& digraph, const char *fn);

     DigraphReader(DigraphReader& other)

       : _is(other._is), local_is(other.local_is), _digraph(other._digraph),

         _use_nodes(other._use_nodes), _use_arcs(other._use_arcs),

         _skip_nodes(other._skip_nodes), _skip_arcs(other._skip_arcs) {

       other._is = 0;

       other.local_is = false;

       _node_index.swap(other._node_index);

       _arc_index.swap(other._arc_index);

       _node_maps.swap(other._node_maps);

       _arc_maps.swap(other._arc_maps);

       _attributes.swap(other._attributes);

       _nodes_caption = other._nodes_caption;

       _arcs_caption = other._arcs_caption;

       _attributes_caption = other._attributes_caption;

     }

     DigraphReader& operator=(const DigraphReader&);

   public:

     /// \name Reading rules

     /// @{

     /// \brief Node map reading rule

     ///

     /// Add a node map reading rule to the reader.

     template <typename Map>

     DigraphReader& nodeMap(const std::string& caption, Map& map) {

       checkConcept<concepts::WriteMap<Node, typename Map::Value>, Map>();

       _reader_bits::MapStorageBase<Node>* storage =

         new _reader_bits::MapStorage<Node, Map>(map);

       _node_maps.push_back(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Node map reading rule

     ///

     /// Add a node map reading rule with specialized converter to the

     /// reader.

     template <typename Map, typename Converter>

     DigraphReader& nodeMap(const std::string& caption, Map& map,

                            const Converter& converter = Converter()) {

       checkConcept<concepts::WriteMap<Node, typename Map::Value>, Map>();

       _reader_bits::MapStorageBase<Node>* storage =

         new _reader_bits::MapStorage<Node, Map, Converter>(map, converter);

       _node_maps.push_back(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Arc map reading rule

     ///

     /// Add an arc map reading rule to the reader.

     template <typename Map>

     DigraphReader& arcMap(const std::string& caption, Map& map) {

       checkConcept<concepts::WriteMap<Arc, typename Map::Value>, Map>();

       _reader_bits::MapStorageBase<Arc>* storage =

         new _reader_bits::MapStorage<Arc, Map>(map);

       _arc_maps.push_back(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Arc map reading rule

     ///

     /// Add an arc map reading rule with specialized converter to the

     /// reader.

     template <typename Map, typename Converter>

     DigraphReader& arcMap(const std::string& caption, Map& map,

                           const Converter& converter = Converter()) {

       checkConcept<concepts::WriteMap<Arc, typename Map::Value>, Map>();

       _reader_bits::MapStorageBase<Arc>* storage =

         new _reader_bits::MapStorage<Arc, Map, Converter>(map, converter);

       _arc_maps.push_back(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Attribute reading rule

     ///

     /// Add an attribute reading rule to the reader.

     template <typename Value>

     DigraphReader& attribute(const std::string& caption, Value& value) {

       _reader_bits::ValueStorageBase* storage =

         new _reader_bits::ValueStorage<Value>(value);

       _attributes.insert(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Attribute reading rule

     ///

     /// Add an attribute reading rule with specialized converter to the

     /// reader.

     template <typename Value, typename Converter>

     DigraphReader& attribute(const std::string& caption, Value& value,

                              const Converter& converter = Converter()) {

       _reader_bits::ValueStorageBase* storage =

         new _reader_bits::ValueStorage<Value, Converter>(value, converter);

       _attributes.insert(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Node reading rule

     ///

     /// Add a node reading rule to reader.

     DigraphReader& node(const std::string& caption, Node& node) {

       typedef _reader_bits::MapLookUpConverter<Node> Converter;

       Converter converter(_node_index);

       _reader_bits::ValueStorageBase* storage =

         new _reader_bits::ValueStorage<Node, Converter>(node, converter);

       _attributes.insert(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Arc reading rule

     ///

     /// Add an arc reading rule to reader.

     DigraphReader& arc(const std::string& caption, Arc& arc) {

       typedef _reader_bits::MapLookUpConverter<Arc> Converter;

       Converter converter(_arc_index);

       _reader_bits::ValueStorageBase* storage =

         new _reader_bits::ValueStorage<Arc, Converter>(arc, converter);

       _attributes.insert(std::make_pair(caption, storage));

       return *this;

     }

     /// @}

     /// \name Select section by name

     /// @{

     /// \brief Set \c \@nodes section to be read

     ///

     /// Set \c \@nodes section to be read

     DigraphReader& nodes(const std::string& caption) {

       _nodes_caption = caption;

       return *this;

     }

     /// \brief Set \c \@arcs section to be read

     ///

     /// Set \c \@arcs section to be read

     DigraphReader& arcs(const std::string& caption) {

       _arcs_caption = caption;

       return *this;

     }

     /// \brief Set \c \@attributes section to be read

     ///

     /// Set \c \@attributes section to be read

     DigraphReader& attributes(const std::string& caption) {

       _attributes_caption = caption;

       return *this;

     }

     /// @}

     /// \name Using previously constructed node or arc set

     /// @{

     /// \brief Use previously constructed node set

     ///

     /// Use previously constructed node set, and specify the node

     /// label map.

     template <typename Map>

     DigraphReader& useNodes(const Map& map) {

       checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();

       LEMON_ASSERT(!_use_nodes, "Multiple usage of useNodes() member");

       _use_nodes = true;

       _writer_bits::DefaultConverter<typename Map::Value> converter;

       for (NodeIt n(_digraph); n != INVALID; ++n) {

         _node_index.insert(std::make_pair(converter(map[n]), n));

       }

       return *this;

     }

     /// \brief Use previously constructed node set

     ///

     /// Use previously constructed node set, and specify the node

     /// label map and a functor which converts the label map values to

     /// \c std::string.

     template <typename Map, typename Converter>

     DigraphReader& useNodes(const Map& map,

                             const Converter& converter = Converter()) {

       checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();

       LEMON_ASSERT(!_use_nodes, "Multiple usage of useNodes() member");

       _use_nodes = true;

       for (NodeIt n(_digraph); n != INVALID; ++n) {

         _node_index.insert(std::make_pair(converter(map[n]), n));

       }

       return *this;

     }

     /// \brief Use previously constructed arc set

     ///

     /// Use previously constructed arc set, and specify the arc

     /// label map.

     template <typename Map>

     DigraphReader& useArcs(const Map& map) {

       checkConcept<concepts::ReadMap<Arc, typename Map::Value>, Map>();

       LEMON_ASSERT(!_use_arcs, "Multiple usage of useArcs() member");

       _use_arcs = true;

       _writer_bits::DefaultConverter<typename Map::Value> converter;

       for (ArcIt a(_digraph); a != INVALID; ++a) {

         _arc_index.insert(std::make_pair(converter(map[a]), a));

       }

       return *this;

     }

     /// \brief Use previously constructed arc set

     ///

     /// Use previously constructed arc set, and specify the arc

     /// label map and a functor which converts the label map values to

     /// \c std::string.

     template <typename Map, typename Converter>

     DigraphReader& useArcs(const Map& map,

                            const Converter& converter = Converter()) {

       checkConcept<concepts::ReadMap<Arc, typename Map::Value>, Map>();

       LEMON_ASSERT(!_use_arcs, "Multiple usage of useArcs() member");

       _use_arcs = true;

       for (ArcIt a(_digraph); a != INVALID; ++a) {

         _arc_index.insert(std::make_pair(converter(map[a]), a));

       }

       return *this;

     }

     /// \brief Skips the reading of node section

     ///

     /// Omit the reading of the node section. This implies that each node

     /// map reading rule will be abandoned, and the nodes of the graph

     /// will not be constructed, which usually cause that the arc set

     /// could not be read due to lack of node name resolving.

     /// Therefore \c skipArcs() function should also be used, or

     /// \c useNodes() should be used to specify the label of the nodes.

     DigraphReader& skipNodes() {

       LEMON_ASSERT(!_skip_nodes, "Skip nodes already set");

       _skip_nodes = true;

       return *this;

     }

     /// \brief Skips the reading of arc section

     ///

     /// Omit the reading of the arc section. This implies that each arc

     /// map reading rule will be abandoned, and the arcs of the graph

     /// will not be constructed.

     DigraphReader& skipArcs() {

       LEMON_ASSERT(!_skip_arcs, "Skip arcs already set");

       _skip_arcs = true;

       return *this;

     }

     /// @}

   private:

     bool readLine() {

       std::string str;

       while(++line_num, std::getline(*_is, str)) {

         line.clear(); line.str(str);

         char c;

         if (line >> std::ws >> c && c != '#') {

           line.putback(c);

           return true;

         }

       }

       return false;

     }

     bool readSuccess() {

       return static_cast<bool>(*_is);

     }

     void skipSection() {

       char c;

       while (readSuccess() && line >> c && c != '@') {

         readLine();

       }

       line.putback(c);

     }

     void readNodes() {

       std::vector<int> map_index(_node_maps.size());

       int map_num, label_index;

       char c;

       if (!readLine() || !(line >> c) || c == '@') {

         if (readSuccess() && line) line.putback(c);

         if (!_node_maps.empty())

           throw FormatError("Cannot find map names");

         return;

       }

       line.putback(c);

       {

         std::map<std::string, int> maps;

         std::string map;

         int index = 0;

         while (_reader_bits::readToken(line, map)) {

           if (maps.find(map) != maps.end()) {

             std::ostringstream msg;

             msg << "Multiple occurence of node map: " << map;

             throw FormatError(msg.str());

           }

           maps.insert(std::make_pair(map, index));

           ++index;

         }

         for (int i = 0; i < static_cast<int>(_node_maps.size()); ++i) {

           std::map<std::string, int>::iterator jt =

             maps.find(_node_maps[i].first);

           if (jt == maps.end()) {

             std::ostringstream msg;

             msg << "Map not found: " << _node_maps[i].first;

             throw FormatError(msg.str());

           }

           map_index[i] = jt->second;

         }

         {

           std::map<std::string, int>::iterator jt = maps.find("label");

           if (jt != maps.end()) {

             label_index = jt->second;

           } else {

             label_index = -1;

           }

         }

         map_num = maps.size();

       }

       while (readLine() && line >> c && c != '@') {

         line.putback(c);

         std::vector<std::string> tokens(map_num);

         for (int i = 0; i < map_num; ++i) {

           if (!_reader_bits::readToken(line, tokens[i])) {

             std::ostringstream msg;

             msg << "Column not found (" << i + 1 << ")";

             throw FormatError(msg.str());

           }

         }

         if (line >> std::ws >> c)

           throw FormatError("Extra character at the end of line");

         Node n;

         if (!_use_nodes) {

           n = _digraph.addNode();

           if (label_index != -1)

             _node_index.insert(std::make_pair(tokens[label_index], n));

         } else {

           if (label_index == -1)

             throw FormatError("Label map not found");

           typename std::map<std::string, Node>::iterator it =

             _node_index.find(tokens[label_index]);

           if (it == _node_index.end()) {

             std::ostringstream msg;

             msg << "Node with label not found: " << tokens[label_index];

             throw FormatError(msg.str());

           }

           n = it->second;

         }

         for (int i = 0; i < static_cast<int>(_node_maps.size()); ++i) {

           _node_maps[i].second->set(n, tokens[map_index[i]]);

         }

       }

       if (readSuccess()) {

         line.putback(c);

       }

     }

     void readArcs() {

       std::vector<int> map_index(_arc_maps.size());

       int map_num, label_index;

       char c;

       if (!readLine() || !(line >> c) || c == '@') {

         if (readSuccess() && line) line.putback(c);

         if (!_arc_maps.empty())

           throw FormatError("Cannot find map names");

         return;

       }

       line.putback(c);

       {

         std::map<std::string, int> maps;

         std::string map;

         int index = 0;

         while (_reader_bits::readToken(line, map)) {

           if (maps.find(map) != maps.end()) {

             std::ostringstream msg;

             msg << "Multiple occurence of arc map: " << map;

             throw FormatError(msg.str());

           }

           maps.insert(std::make_pair(map, index));

           ++index;

         }

         for (int i = 0; i < static_cast<int>(_arc_maps.size()); ++i) {

           std::map<std::string, int>::iterator jt =

             maps.find(_arc_maps[i].first);

           if (jt == maps.end()) {

             std::ostringstream msg;

             msg << "Map not found: " << _arc_maps[i].first;

             throw FormatError(msg.str());

           }

           map_index[i] = jt->second;

         }

         {

           std::map<std::string, int>::iterator jt = maps.find("label");

           if (jt != maps.end()) {

             label_index = jt->second;

           } else {

             label_index = -1;

           }

         }

         map_num = maps.size();

       }

       while (readLine() && line >> c && c != '@') {

         line.putback(c);

         std::string source_token;

         std::string target_token;

         if (!_reader_bits::readToken(line, source_token))

           throw FormatError("Source not found");

         if (!_reader_bits::readToken(line, target_token))

           throw FormatError("Target not found");

         std::vector<std::string> tokens(map_num);

         for (int i = 0; i < map_num; ++i) {

           if (!_reader_bits::readToken(line, tokens[i])) {

             std::ostringstream msg;

             msg << "Column not found (" << i + 1 << ")";

             throw FormatError(msg.str());

           }

         }

         if (line >> std::ws >> c)

           throw FormatError("Extra character at the end of line");

         Arc a;

         if (!_use_arcs) {

           typename NodeIndex::iterator it;

           it = _node_index.find(source_token);

           if (it == _node_index.end()) {

             std::ostringstream msg;

             msg << "Item not found: " << source_token;

             throw FormatError(msg.str());

           }

           Node source = it->second;

           it = _node_index.find(target_token);

           if (it == _node_index.end()) {

             std::ostringstream msg;

             msg << "Item not found: " << target_token;

             throw FormatError(msg.str());

           }

           Node target = it->second;

           a = _digraph.addArc(source, target);

           if (label_index != -1)

             _arc_index.insert(std::make_pair(tokens[label_index], a));

         } else {

           if (label_index == -1)

             throw FormatError("Label map not found");

           typename std::map<std::string, Arc>::iterator it =

             _arc_index.find(tokens[label_index]);

           if (it == _arc_index.end()) {

             std::ostringstream msg;

             msg << "Arc with label not found: " << tokens[label_index];

             throw FormatError(msg.str());

           }

           a = it->second;

         }

         for (int i = 0; i < static_cast<int>(_arc_maps.size()); ++i) {

           _arc_maps[i].second->set(a, tokens[map_index[i]]);

         }

       }

       if (readSuccess()) {

         line.putback(c);

       }

     }

     void readAttributes() {

       std::set<std::string> read_attr;

       char c;

       while (readLine() && line >> c && c != '@') {

         line.putback(c);

         std::string attr, token;

         if (!_reader_bits::readToken(line, attr))

           throw FormatError("Attribute name not found");

         if (!_reader_bits::readToken(line, token))

           throw FormatError("Attribute value not found");

         if (line >> c)

           throw FormatError("Extra character at the end of line");

         {

           std::set<std::string>::iterator it = read_attr.find(attr);

           if (it != read_attr.end()) {

             std::ostringstream msg;

             msg << "Multiple occurence of attribute: " << attr;

             throw FormatError(msg.str());

           }

           read_attr.insert(attr);

         }

         {

           typename Attributes::iterator it = _attributes.lower_bound(attr);

           while (it != _attributes.end() && it->first == attr) {

             it->second->set(token);

             ++it;

           }

         }

       }

       if (readSuccess()) {

         line.putback(c);

       }

       for (typename Attributes::iterator it = _attributes.begin();

            it != _attributes.end(); ++it) {

         if (read_attr.find(it->first) == read_attr.end()) {

           std::ostringstream msg;

           msg << "Attribute not found: " << it->first;

           throw FormatError(msg.str());

         }

       }

     }

   public:

     /// \name Execution of the reader

     /// @{

     /// \brief Start the batch processing

     ///

     /// This function starts the batch processing

     void run() {

       LEMON_ASSERT(_is != 0, "This reader assigned to an other reader");

       bool nodes_done = _skip_nodes;

       bool arcs_done = _skip_arcs;

       bool attributes_done = false;

       line_num = 0;

       readLine();

       skipSection();

       while (readSuccess()) {

         try {

           char c;

           std::string section, caption;

           line >> c;

           _reader_bits::readToken(line, section);

           _reader_bits::readToken(line, caption);

           if (line >> c)

             throw FormatError("Extra character at the end of line");

           if (section == "nodes" && !nodes_done) {

             if (_nodes_caption.empty() || _nodes_caption == caption) {

               readNodes();

               nodes_done = true;

             }

           } else if ((section == "arcs" || section == "edges") &&

                      !arcs_done) {

             if (_arcs_caption.empty() || _arcs_caption == caption) {

               readArcs();

               arcs_done = true;

             }

           } else if (section == "attributes" && !attributes_done) {

             if (_attributes_caption.empty() || _attributes_caption == caption) {

               readAttributes();

               attributes_done = true;

             }

           } else {

             readLine();

             skipSection();

           }

         } catch (FormatError& error) {

           error.line(line_num);

           error.file(_filename);

           throw;

         }

       }

       if (!nodes_done) {

         throw FormatError("Section @nodes not found");

       }

       if (!arcs_done) {

         throw FormatError("Section @arcs not found");

       }

       if (!attributes_done && !_attributes.empty()) {

         throw FormatError("Section @attributes not found");

       }

     }

     /// @}

   };

   /// \brief Return a \ref DigraphReader class

   ///

   /// This function just returns a \ref DigraphReader class.

   /// \relates DigraphReader

   template <typename Digraph>

   DigraphReader<Digraph> digraphReader(Digraph& digraph, std::istream& is) {

     DigraphReader<Digraph> tmp(digraph, is);

     return tmp;

   }

   /// \brief Return a \ref DigraphReader class

   ///

   /// This function just returns a \ref DigraphReader class.

   /// \relates DigraphReader

   template <typename Digraph>

   DigraphReader<Digraph> digraphReader(Digraph& digraph,

                                        const std::string& fn) {

     DigraphReader<Digraph> tmp(digraph, fn);

     return tmp;

   }

   /// \brief Return a \ref DigraphReader class

   ///

   /// This function just returns a \ref DigraphReader class.

   /// \relates DigraphReader

   template <typename Digraph>

   DigraphReader<Digraph> digraphReader(Digraph& digraph, const char* fn) {

     DigraphReader<Digraph> tmp(digraph, fn);

     return tmp;

   }

   template <typename Graph>

   class GraphReader;

   template <typename Graph>

   GraphReader<Graph> graphReader(Graph& graph, 

                                  std::istream& is = std::cin);

   template <typename Graph>

   GraphReader<Graph> graphReader(Graph& graph, const std::string& fn);

   template <typename Graph>

   GraphReader<Graph> graphReader(Graph& graph, const char *fn);

   /// \ingroup lemon_io

   ///

   /// \brief \ref lgf-format "LGF" reader for undirected graphs

   ///

   /// This utility reads an \ref lgf-format "LGF" file.

   ///

   /// It can be used almost the same way as \c DigraphReader.

   /// The only difference is that this class can handle edges and

   /// edge maps as well as arcs and arc maps.

   ///

   /// The columns in the \c \@edges (or \c \@arcs) section are the

   /// edge maps. However, if there are two maps with the same name

   /// prefixed with \c '+' and \c '-', then these can be read into an

   /// arc map.  Similarly, an attribute can be read into an arc, if

   /// it's value is an edge label prefixed with \c '+' or \c '-'.

   template <typename _Graph>

   class GraphReader {

   public:

     typedef _Graph Graph;

     TEMPLATE_GRAPH_TYPEDEFS(Graph);

   private:

     std::istream* _is;

     bool local_is;

     std::string _filename;

     Graph& _graph;

     std::string _nodes_caption;

     std::string _edges_caption;

     std::string _attributes_caption;

     typedef std::map<std::string, Node> NodeIndex;

     NodeIndex _node_index;

     typedef std::map<std::string, Edge> EdgeIndex;

     EdgeIndex _edge_index;

     typedef std::vector<std::pair<std::string,

       _reader_bits::MapStorageBase<Node>*> > NodeMaps;

     NodeMaps _node_maps;

     typedef std::vector<std::pair<std::string,

       _reader_bits::MapStorageBase<Edge>*> > EdgeMaps;

     EdgeMaps _edge_maps;

     typedef std::multimap<std::string, _reader_bits::ValueStorageBase*>

       Attributes;

     Attributes _attributes;

     bool _use_nodes;

     bool _use_edges;

     bool _skip_nodes;

     bool _skip_edges;

     int line_num;

     std::istringstream line;

   public:

     /// \brief Constructor

     ///

     /// Construct an undirected graph reader, which reads from the given

     /// input stream.

     GraphReader(Graph& graph, std::istream& is = std::cin)

       : _is(&is), local_is(false), _graph(graph),

         _use_nodes(false), _use_edges(false),

         _skip_nodes(false), _skip_edges(false) {}

     /// \brief Constructor

     ///

     /// Construct an undirected graph reader, which reads from the given

     /// file.

     GraphReader(Graph& graph, const std::string& fn)

       : _is(new std::ifstream(fn.c_str())), local_is(true),

         _filename(fn), _graph(graph),

         _use_nodes(false), _use_edges(false),

         _skip_nodes(false), _skip_edges(false) {

       if (!(*_is)) {

         delete _is;

         throw IoError("Cannot open file", fn);

       }

     }

     /// \brief Constructor

     ///

     /// Construct an undirected graph reader, which reads from the given

     /// file.

     GraphReader(Graph& graph, const char* fn)

       : _is(new std::ifstream(fn)), local_is(true),

         _filename(fn), _graph(graph),

         _use_nodes(false), _use_edges(false),

         _skip_nodes(false), _skip_edges(false) {

       if (!(*_is)) {

         delete _is;

         throw IoError("Cannot open file", fn);

       }

     }

     /// \brief Destructor

     ~GraphReader() {

       for (typename NodeMaps::iterator it = _node_maps.begin();

            it != _node_maps.end(); ++it) {

         delete it->second;

       }

       for (typename EdgeMaps::iterator it = _edge_maps.begin();

            it != _edge_maps.end(); ++it) {

         delete it->second;

       }

       for (typename Attributes::iterator it = _attributes.begin();

            it != _attributes.end(); ++it) {

         delete it->second;

       }

       if (local_is) {

         delete _is;

       }

     }

   private:

     template <typename GR>

     friend GraphReader<GR> graphReader(GR& graph, std::istream& is);

     template <typename GR>

     friend GraphReader<GR> graphReader(GR& graph, const std::string& fn); 

     template <typename GR>

     friend GraphReader<GR> graphReader(GR& graph, const char *fn);

     GraphReader(GraphReader& other)

       : _is(other._is), local_is(other.local_is), _graph(other._graph),

         _use_nodes(other._use_nodes), _use_edges(other._use_edges),

         _skip_nodes(other._skip_nodes), _skip_edges(other._skip_edges) {

       other._is = 0;

       other.local_is = false;

       _node_index.swap(other._node_index);

       _edge_index.swap(other._edge_index);

       _node_maps.swap(other._node_maps);

       _edge_maps.swap(other._edge_maps);

       _attributes.swap(other._attributes);

       _nodes_caption = other._nodes_caption;

       _edges_caption = other._edges_caption;

       _attributes_caption = other._attributes_caption;

     }

     GraphReader& operator=(const GraphReader&);

   public:

     /// \name Reading rules

     /// @{

     /// \brief Node map reading rule

     ///

     /// Add a node map reading rule to the reader.

     template <typename Map>

     GraphReader& nodeMap(const std::string& caption, Map& map) {

       checkConcept<concepts::WriteMap<Node, typename Map::Value>, Map>();

       _reader_bits::MapStorageBase<Node>* storage =

         new _reader_bits::MapStorage<Node, Map>(map);

       _node_maps.push_back(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Node map reading rule

     ///

     /// Add a node map reading rule with specialized converter to the

     /// reader.

     template <typename Map, typename Converter>

     GraphReader& nodeMap(const std::string& caption, Map& map,

                            const Converter& converter = Converter()) {

       checkConcept<concepts::WriteMap<Node, typename Map::Value>, Map>();

       _reader_bits::MapStorageBase<Node>* storage =

         new _reader_bits::MapStorage<Node, Map, Converter>(map, converter);

       _node_maps.push_back(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Edge map reading rule

     ///

     /// Add an edge map reading rule to the reader.

     template <typename Map>

     GraphReader& edgeMap(const std::string& caption, Map& map) {

       checkConcept<concepts::WriteMap<Edge, typename Map::Value>, Map>();

       _reader_bits::MapStorageBase<Edge>* storage =

         new _reader_bits::MapStorage<Edge, Map>(map);

       _edge_maps.push_back(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Edge map reading rule

     ///

     /// Add an edge map reading rule with specialized converter to the

     /// reader.

     template <typename Map, typename Converter>

     GraphReader& edgeMap(const std::string& caption, Map& map,

                           const Converter& converter = Converter()) {

       checkConcept<concepts::WriteMap<Edge, typename Map::Value>, Map>();

       _reader_bits::MapStorageBase<Edge>* storage =

         new _reader_bits::MapStorage<Edge, Map, Converter>(map, converter);

       _edge_maps.push_back(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Arc map reading rule

     ///

     /// Add an arc map reading rule to the reader.

     template <typename Map>

     GraphReader& arcMap(const std::string& caption, Map& map) {

       checkConcept<concepts::WriteMap<Arc, typename Map::Value>, Map>();

       _reader_bits::MapStorageBase<Edge>* forward_storage =

         new _reader_bits::GraphArcMapStorage<Graph, true, Map>(_graph, map);

       _edge_maps.push_back(std::make_pair('+' + caption, forward_storage));

       _reader_bits::MapStorageBase<Edge>* backward_storage =

         new _reader_bits::GraphArcMapStorage<Graph, false, Map>(_graph, map);

       _edge_maps.push_back(std::make_pair('-' + caption, backward_storage));

       return *this;

     }

     /// \brief Arc map reading rule

     ///

     /// Add an arc map reading rule with specialized converter to the

     /// reader.

     template <typename Map, typename Converter>

     GraphReader& arcMap(const std::string& caption, Map& map,

                           const Converter& converter = Converter()) {

       checkConcept<concepts::WriteMap<Arc, typename Map::Value>, Map>();

       _reader_bits::MapStorageBase<Edge>* forward_storage =

         new _reader_bits::GraphArcMapStorage<Graph, true, Map, Converter>

         (_graph, map, converter);

       _edge_maps.push_back(std::make_pair('+' + caption, forward_storage));

       _reader_bits::MapStorageBase<Edge>* backward_storage =

         new _reader_bits::GraphArcMapStorage<Graph, false, Map, Converter>

         (_graph, map, converter);

       _edge_maps.push_back(std::make_pair('-' + caption, backward_storage));

       return *this;

     }

     /// \brief Attribute reading rule

     ///

     /// Add an attribute reading rule to the reader.

     template <typename Value>

     GraphReader& attribute(const std::string& caption, Value& value) {

       _reader_bits::ValueStorageBase* storage =

         new _reader_bits::ValueStorage<Value>(value);

       _attributes.insert(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Attribute reading rule

     ///

     /// Add an attribute reading rule with specialized converter to the

     /// reader.

     template <typename Value, typename Converter>

     GraphReader& attribute(const std::string& caption, Value& value,

                              const Converter& converter = Converter()) {

       _reader_bits::ValueStorageBase* storage =

         new _reader_bits::ValueStorage<Value, Converter>(value, converter);

       _attributes.insert(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Node reading rule

     ///

     /// Add a node reading rule to reader.

     GraphReader& node(const std::string& caption, Node& node) {

       typedef _reader_bits::MapLookUpConverter<Node> Converter;

       Converter converter(_node_index);

       _reader_bits::ValueStorageBase* storage =

         new _reader_bits::ValueStorage<Node, Converter>(node, converter);

       _attributes.insert(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Edge reading rule

     ///

     /// Add an edge reading rule to reader.

     GraphReader& edge(const std::string& caption, Edge& edge) {

       typedef _reader_bits::MapLookUpConverter<Edge> Converter;

       Converter converter(_edge_index);

       _reader_bits::ValueStorageBase* storage =

         new _reader_bits::ValueStorage<Edge, Converter>(edge, converter);

       _attributes.insert(std::make_pair(caption, storage));

       return *this;

     }

     /// \brief Arc reading rule

     ///

     /// Add an arc reading rule to reader.

     GraphReader& arc(const std::string& caption, Arc& arc) {

       typedef _reader_bits::GraphArcLookUpConverter<Graph> Converter;

       Converter converter(_graph, _edge_index);

       _reader_bits::ValueStorageBase* storage =

         new _reader_bits::ValueStorage<Arc, Converter>(arc, converter);

       _attributes.insert(std::make_pair(caption, storage));

       return *this;

     }

     /// @}

     /// \name Select section by name

     /// @{

     /// \brief Set \c \@nodes section to be read

     ///

     /// Set \c \@nodes section to be read.

     GraphReader& nodes(const std::string& caption) {

       _nodes_caption = caption;

       return *this;

     }

     /// \brief Set \c \@edges section to be read

     ///

     /// Set \c \@edges section to be read.

     GraphReader& edges(const std::string& caption) {

       _edges_caption = caption;

       return *this;

     }

     /// \brief Set \c \@attributes section to be read

     ///

     /// Set \c \@attributes section to be read.

     GraphReader& attributes(const std::string& caption) {

       _attributes_caption = caption;

       return *this;

     }

     /// @}

     /// \name Using previously constructed node or edge set

     /// @{

     /// \brief Use previously constructed node set

     ///

     /// Use previously constructed node set, and specify the node

     /// label map.

     template <typename Map>

     GraphReader& useNodes(const Map& map) {

       checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();

       LEMON_ASSERT(!_use_nodes, "Multiple usage of useNodes() member");

       _use_nodes = true;

       _writer_bits::DefaultConverter<typename Map::Value> converter;

       for (NodeIt n(_graph); n != INVALID; ++n) {

         _node_index.insert(std::make_pair(converter(map[n]), n));

       }

       return *this;

     }

     /// \brief Use previously constructed node set

     ///

     /// Use previously constructed node set, and specify the node

     /// label map and a functor which converts the label map values to

     /// \c std::string.

     template <typename Map, typename Converter>

     GraphReader& useNodes(const Map& map,

                             const Converter& converter = Converter()) {

       checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();

       LEMON_ASSERT(!_use_nodes, "Multiple usage of useNodes() member");

       _use_nodes = true;

       for (NodeIt n(_graph); n != INVALID; ++n) {

         _node_index.insert(std::make_pair(converter(map[n]), n));

       }

       return *this;

     }

     /// \brief Use previously constructed edge set

     ///

     /// Use previously constructed edge set, and specify the edge

     /// label map.

     template <typename Map>

     GraphReader& useEdges(const Map& map) {

       checkConcept<concepts::ReadMap<Edge, typename Map::Value>, Map>();

       LEMON_ASSERT(!_use_edges, "Multiple usage of useEdges() member");

       _use_edges = true;

       _writer_bits::DefaultConverter<typename Map::Value> converter;

       for (EdgeIt a(_graph); a != INVALID; ++a) {

         _edge_index.insert(std::make_pair(converter(map[a]), a));

       }

       return *this;

     }

     /// \brief Use previously constructed edge set

     ///

     /// Use previously constructed edge set, and specify the edge

     /// label map and a functor which converts the label map values to

     /// \c std::string.

     template <typename Map, typename Converter>

     GraphReader& useEdges(const Map& map,

                             const Converter& converter = Converter()) {

       checkConcept<concepts::ReadMap<Edge, typename Map::Value>, Map>();

       LEMON_ASSERT(!_use_edges, "Multiple usage of useEdges() member");

       _use_edges = true;

       for (EdgeIt a(_graph); a != INVALID; ++a) {

         _edge_index.insert(std::make_pair(converter(map[a]), a));

       }

       return *this;

     }

     /// \brief Skip the reading of node section

     ///

     /// Omit the reading of the node section. This implies that each node

     /// map reading rule will be abandoned, and the nodes of the graph

     /// will not be constructed, which usually cause that the edge set

     /// could not be read due to lack of node name

     /// could not be read due to lack of node name resolving.

     /// Therefore \c skipEdges() function should also be used, or

     /// \c useNodes() should be used to specify the label of the nodes.

     GraphReader& skipNodes() {

       LEMON_ASSERT(!_skip_nodes, "Skip nodes already set");

       _skip_nodes = true;

       return *this;

     }

     /// \brief Skip the reading of edge section

     ///

     /// Omit the reading of the edge section. This implies that each edge

     /// map reading rule will be abandoned, and the edges of the graph

     /// will not be constructed.

     GraphReader& skipEdges() {

       LEMON_ASSERT(!_skip_edges, "Skip edges already set");

       _skip_edges = true;

       return *this;

     }

     /// @}

   private:

     bool readLine() {

       std::string str;

       while(++line_num, std::getline(*_is, str)) {

         line.clear(); line.str(str);

         char c;

         if (line >> std::ws >> c && c != '#') {

           line.putback(c);

           return true;

         }

       }

       return false;

     }

     bool readSuccess() {

       return static_cast<bool>(*_is);

     }

     void skipSection() {

       char c;

       while (readSuccess() && line >> c && c != '@') {

         readLine();

       }

       line.putback(c);

     }

     void readNodes() {

       std::vector<int> map_index(_node_maps.size());

       int map_num, label_index;

       char c;

       if (!readLine() || !(line >> c) || c == '@') {

         if (readSuccess() && line) line.putback(c);

         if (!_node_maps.empty())

           throw FormatError("Cannot find map names");

         return;

       }

       line.putback(c);

       {

         std::map<std::string, int> maps;

         std::string map;

         int index = 0;

         while (_reader_bits::readToken(line, map)) {

           if (maps.find(map) != maps.end()) {

             std::ostringstream msg;

             msg << "Multiple occurence of node map: " << map;

             throw FormatError(msg.str());

           }

           maps.insert(std::make_pair(map, index));

           ++index;

         }

         for (int i = 0; i < static_cast<int>(_node_maps.size()); ++i) {

           std::map<std::string, int>::iterator jt =

             maps.find(_node_maps[i].first);

           if (jt == maps.end()) {

             std::ostringstream msg;

             msg << "Map not found: " << _node_maps[i].first;

             throw FormatError(msg.str());

           }

           map_index[i] = jt->second;

         }

         {

           std::map<std::string, int>::iterator jt = maps.find("label");

           if (jt != maps.end()) {

             label_index = jt->second;

           } else {

             label_index = -1;

           }

         }

         map_num = maps.size();

       }

       while (readLine() && line >> c && c != '@') {

         line.putback(c);

         std::vector<std::string> tokens(map_num);

         for (int i = 0; i < map_num; ++i) {

           if (!_reader_bits::readToken(line, tokens[i])) {

             std::ostringstream msg;

             msg << "Column not found (" << i + 1 << ")";

             throw FormatError(msg.str());

           }

         }

         if (line >> std::ws >> c)

           throw FormatError("Extra character at the end of line");

         Node n;

         if (!_use_nodes) {

           n = _graph.addNode();

           if (label_index != -1)

             _node_index.insert(std::make_pair(tokens[label_index], n));

         } else {

           if (label_index == -1)

             throw FormatError("Label map not found");

           typename std::map<std::string, Node>::iterator it =

             _node_index.find(tokens[label_index]);

           if (it == _node_index.end()) {

             std::ostringstream msg;

             msg << "Node with label not found: " << tokens[label_index];

             throw FormatError(msg.str());

           }

           n = it->second;

         }

         for (int i = 0; i < static_cast<int>(_node_maps.size()); ++i) {

           _node_maps[i].second->set(n, tokens[map_index[i]]);

         }

       }

       if (readSuccess()) {

         line.putback(c);

       }

     }

     void readEdges() {

       std::vector<int> map_index(_edge_maps.size());

       int map_num, label_index;

       char c;

       if (!readLine() || !(line >> c) || c == '@') {

         if (readSuccess() && line) line.putback(c);

         if (!_edge_maps.empty())

           throw FormatError("Cannot find map names");

         return;

       }

       line.putback(c);

       {

         std::map<std::string, int> maps;

         std::string map;

         int index = 0;

         while (_reader_bits::readToken(line, map)) {

           if (maps.find(map) != maps.end()) {

             std::ostringstream msg;

             msg << "Multiple occurence of edge map: " << map;

             throw FormatError(msg.str());

           }

           maps.insert(std::make_pair(map, index));

           ++index;

         }

         for (int i = 0; i < static_cast<int>(_edge_maps.size()); ++i) {

           std::map<std::string, int>::iterator jt =

             maps.find(_edge_maps[i].first);

           if (jt == maps.end()) {

             std::ostringstream msg;

             msg << "Map not found: " << _edge_maps[i].first;

             throw FormatError(msg.str());

           }

           map_index[i] = jt->second;

         }

         {

           std::map<std::string, int>::iterator jt = maps.find("label");

           if (jt != maps.end()) {

             label_index = jt->second;

           } else {

             label_index = -1;

           }

         }

         map_num = maps.size();

       }

       while (readLine() && line >> c && c != '@') {

         line.putback(c);

         std::string source_token;

         std::string target_token;

         if (!_reader_bits::readToken(line, source_token))

           throw FormatError("Node u not found");

         if (!_reader_bits::readToken(line, target_token))

           throw FormatError("Node v not found");

         std::vector<std::string> tokens(map_num);

         for (int i = 0; i < map_num; ++i) {

           if (!_reader_bits::readToken(line, tokens[i])) {

             std::ostringstream msg;

             msg << "Column not found (" << i + 1 << ")";

             throw FormatError(msg.str());

           }

         }

         if (line >> std::ws >> c)

           throw FormatError("Extra character at the end of line");

         Edge e;

         if (!_use_edges) {

           typename NodeIndex::iterator it;

           it = _node_index.find(source_token);

           if (it == _node_index.end()) {

             std::ostringstream msg;

             msg << "Item not found: " << source_token;

             throw FormatError(msg.str());

           }

           Node source = it->second;

           it = _node_index.find(target_token);

           if (it == _node_index.end()) {

             std::ostringstream msg;

             msg << "Item not found: " << target_token;

             throw FormatError(msg.str());

           }

           Node target = it->second;

           e = _graph.addEdge(source, target);

           if (label_index != -1)

             _edge_index.insert(std::make_pair(tokens[label_index], e));

         } else {

           if (label_index == -1)

             throw FormatError("Label map not found");

           typename std::map<std::string, Edge>::iterator it =

             _edge_index.find(tokens[label_index]);

           if (it == _edge_index.end()) {

             std::ostringstream msg;

             msg << "Edge with label not found: " << tokens[label_index];

             throw FormatError(msg.str());

           }

           e = it->second;

         }

         for (int i = 0; i < static_cast<int>(_edge_maps.size()); ++i) {

           _edge_maps[i].second->set(e, tokens[map_index[i]]);

         }

       }

       if (readSuccess()) {

         line.putback(c);

       }

     }

     void readAttributes() {

       std::set<std::string> read_attr;

       char c;

       while (readLine() && line >> c && c != '@') {

         line.putback(c);

         std::string attr, token;

         if (!_reader_bits::readToken(line, attr))

           throw FormatError("Attribute name not found");

         if (!_reader_bits::readToken(line, token))

           throw FormatError("Attribute value not found");

         if (line >> c)

           throw FormatError("Extra character at the end of line");

         {

           std::set<std::string>::iterator it = read_attr.find(attr);

           if (it != read_attr.end()) {

             std::ostringstream msg;

             msg << "Multiple occurence of attribute: " << attr;

             throw FormatError(msg.str());

           }

           read_attr.insert(attr);

         }

         {

           typename Attributes::iterator it = _attributes.lower_bound(attr);

           while (it != _attributes.end() && it->first == attr) {

             it->second->set(token);

             ++it;

           }

         }

       }

       if (readSuccess()) {

         line.putback(c);

       }

       for (typename Attributes::iterator it = _attributes.begin();

            it != _attributes.end(); ++it) {

         if (read_attr.find(it->first) == read_attr.end()) {

           std::ostringstream msg;

           msg << "Attribute not found: " << it->first;

           throw FormatError(msg.str());

         }

       }

     }

   public:

     /// \name Execution of the reader

     /// @{

     /// \brief Start the batch processing

     ///

     /// This function starts the batch processing

     void run() {

       LEMON_ASSERT(_is != 0, "This reader assigned to an other reader");

       bool nodes_done = _skip_nodes;

       bool edges_done = _skip_edges;

       bool attributes_done = false;

       line_num = 0;

       readLine();

       skipSection();

       while (readSuccess()) {

         try {

           char c;

           std::string section, caption;

           line >> c;

           _reader_bits::readToken(line, section);

           _reader_bits::readToken(line, caption);

           if (line >> c)

             throw FormatError("Extra character at the end of line");

           if (section == "nodes" && !nodes_done) {

             if (_nodes_caption.empty() || _nodes_caption == caption) {

               readNodes();

               nodes_done = true;

             }

           } else if ((section == "edges" || section == "arcs") &&

                      !edges_done) {

             if (_edges_caption.empty() || _edges_caption == caption) {

               readEdges();

               edges_done = true;

             }

           } else if (section == "attributes" && !attributes_done) {

             if (_attributes_caption.empty() || _attributes_caption == caption) {

               readAttributes();

               attributes_done = true;

             }

           } else {

             readLine();

             skipSection();

           }

         } catch (FormatError& error) {

           error.line(line_num);

           error.file(_filename);

           throw;

         }

       }

       if (!nodes_done) {

         throw FormatError("Section @nodes not found");

       }

       if (!edges_done) {

         throw FormatError("Section @edges not found");

       }

       if (!attributes_done && !_attributes.empty()) {

         throw FormatError("Section @attributes not found");

       }

     }

     /// @}

   };

   /// \brief Return a \ref GraphReader class

   ///

   /// This function just returns a \ref GraphReader class.

   /// \relates GraphReader

   template <typename Graph>

   GraphReader<Graph> graphReader(Graph& graph, std::istream& is) {

     GraphReader<Graph> tmp(graph, is);

     return tmp;

   }

   /// \brief Return a \ref GraphReader class

   ///

   /// This function just returns a \ref GraphReader class.

   /// \relates GraphReader

   template <typename Graph>

   GraphReader<Graph> graphReader(Graph& graph, const std::string& fn) {

     GraphReader<Graph> tmp(graph, fn);

     return tmp;

   }

   /// \brief Return a \ref GraphReader class

   ///

   /// This function just returns a \ref GraphReader class.

   /// \relates GraphReader

   template <typename Graph>

   GraphReader<Graph> graphReader(Graph& graph, const char* fn) {

     GraphReader<Graph> tmp(graph, fn);

     return tmp;

   }

   class SectionReader;

   SectionReader sectionReader(std::istream& is);

   SectionReader sectionReader(const std::string& fn);

   SectionReader sectionReader(const char* fn);

   /// \ingroup lemon_io

   ///

   /// \brief Section reader class

   ///

   /// In the \ref lgf-format "LGF" file extra sections can be placed,

   /// which contain any data in arbitrary format. Such sections can be

   /// read with this class. A reading rule can be added to the class

   /// with two different functions. With the \c sectionLines() function a

   /// functor can process the section line-by-line, while with the \c

   /// sectionStream() member the section can be read from an input

   /// stream.

   class SectionReader {

   private:

     std::istream* _is;

     bool local_is;

     std::string _filename;

     typedef std::map<std::string, _reader_bits::Section*> Sections;

     Sections _sections;

     int line_num;

     std::istringstream line;

   public:

     /// \brief Constructor

     ///

     /// Construct a section reader, which reads from the given input

     /// stream.

     SectionReader(std::istream& is)

       : _is(&is), local_is(false) {}

     /// \brief Constructor

     ///

     /// Construct a section reader, which reads from the given file.

     SectionReader(const std::string& fn)

       : _is(new std::ifstream(fn.c_str())), local_is(true),

         _filename(fn) {

       if (!(*_is)) {

         delete _is;

         throw IoError("Cannot open file", fn);

       }

     }

     /// \brief Constructor

     ///

     /// Construct a section reader, which reads from the given file.

     SectionReader(const char* fn)

       : _is(new std::ifstream(fn)), local_is(true),

         _filename(fn) {

       if (!(*_is)) {

         delete _is;

         throw IoError("Cannot open file", fn);

       }

     }

     /// \brief Destructor

     ~SectionReader() {

       for (Sections::iterator it = _sections.begin();

            it != _sections.end(); ++it) {

         delete it->second;

       }

       if (local_is) {

         delete _is;

       }

     }

   private:

     friend SectionReader sectionReader(std::istream& is);

     friend SectionReader sectionReader(const std::string& fn);

     friend SectionReader sectionReader(const char* fn);

     SectionReader(SectionReader& other)

       : _is(other._is), local_is(other.local_is) {

       other._is = 0;

       other.local_is = false;

       _sections.swap(other._sections);

     }

     SectionReader& operator=(const SectionReader&);

   public:

     /// \name Section readers

     /// @{

     /// \brief Add a section processor with line oriented reading

     ///

     /// The first parameter is the type descriptor of the section, the

     /// second is a functor, which takes just one \c std::string

     /// parameter. At the reading process, each line of the section

     /// will be given to the functor object. However, the empty lines

     /// and the comment lines are filtered out, and the leading

     /// whitespaces are trimmed from each processed string.

     ///

     /// For example let's see a section, which contain several

     /// integers, which should be inserted into a vector.

     ///\code

     ///  @numbers

     ///  12 45 23

     ///  4

     ///  23 6

     ///\endcode

     ///

     /// The functor is implemented as a struct:

     ///\code

     ///  struct NumberSection {

     ///    std::vector<int>& _data;

     ///    NumberSection(std::vector<int>& data) : _data(data) {}

     ///    void operator()(const std::string& line) {

     ///      std::istringstream ls(line);

     ///      int value;

     ///      while (ls >> value) _data.push_back(value);

     ///    }

     ///  };

     ///

     ///  // ...

     ///

     ///  reader.sectionLines("numbers", NumberSection(vec));

     ///\endcode

     template <typename Functor>

     SectionReader& sectionLines(const std::string& type, Functor functor) {

       LEMON_ASSERT(!type.empty(), "Type is empty.");

       LEMON_ASSERT(_sections.find(type) == _sections.end(),

                    "Multiple reading of section.");

       _sections.insert(std::make_pair(type,

         new _reader_bits::LineSection<Functor>(functor)));

       return *this;

     }

     /// \brief Add a section processor with stream oriented reading

     ///

     /// The first parameter is the type of the section, the second is

     /// a functor, which takes an \c std::istream& and an \c int&

     /// parameter, the latter regard to the line number of stream. The

     /// functor can read the input while the section go on, and the

     /// line number should be modified accordingly.

     template <typename Functor>

     SectionReader& sectionStream(const std::string& type, Functor functor) {

       LEMON_ASSERT(!type.empty(), "Type is empty.");

       LEMON_ASSERT(_sections.find(type) == _sections.end(),

                    "Multiple reading of section.");

       _sections.insert(std::make_pair(type,

          new _reader_bits::StreamSection<Functor>(functor)));

       return *this;

     }

     /// @}

   private:

     bool readLine() {

       std::string str;

       while(++line_num, std::getline(*_is, str)) {

         line.clear(); line.str(str);

         char c;

         if (line >> std::ws >> c && c != '#') {

           line.putback(c);

           return true;

         }

       }

       return false;

     }

     bool readSuccess() {

       return static_cast<bool>(*_is);

     }

     void skipSection() {

       char c;

       while (readSuccess() && line >> c && c != '@') {

         readLine();

       }

       line.putback(c);

     }

   public:

     /// \name Execution of the reader

     /// @{

     /// \brief Start the batch processing

     ///

     /// This function starts the batch processing.

     void run() {

       LEMON_ASSERT(_is != 0, "This reader assigned to an other reader");

       std::set<std::string> extra_sections;

       line_num = 0;

       readLine();

       skipSection();

       while (readSuccess()) {

         try {

           char c;

           std::string section, caption;

           line >> c;

           _reader_bits::readToken(line, section);

           _reader_bits::readToken(line, caption);

           if (line >> c)

             throw FormatError("Extra character at the end of line");

           if (extra_sections.find(section) != extra_sections.end()) {

             std::ostringstream msg;

             msg << "Multiple occurence of section: " << section;

             throw FormatError(msg.str());

           }

           Sections::iterator it = _sections.find(section);

           if (it != _sections.end()) {

             extra_sections.insert(section);

             it->second->process(*_is, line_num);

           }

           readLine();

           skipSection();

         } catch (FormatError& error) {

           error.line(line_num);

           error.file(_filename);

           throw;

         }

       }

       for (Sections::iterator it = _sections.begin();

            it != _sections.end(); ++it) {

         if (extra_sections.find(it->first) == extra_sections.end()) {

           std::ostringstream os;

           os << "Cannot find section: " << it->first;

           throw FormatError(os.str());

         }

       }

     }

     /// @}

   };

   /// \brief Return a \ref SectionReader class

   ///

   /// This function just returns a \ref SectionReader class.

   /// \relates SectionReader

   inline SectionReader sectionReader(std::istream& is) {

     SectionReader tmp(is);

     return tmp;

   }

   /// \brief Return a \ref SectionReader class

   ///

   /// This function just returns a \ref SectionReader class.

   /// \relates SectionReader

   inline SectionReader sectionReader(const std::string& fn) {

     SectionReader tmp(fn);

     return tmp;

   }

   /// \brief Return a \ref SectionReader class

   ///

   /// This function just returns a \ref SectionReader class.

   /// \relates SectionReader

   inline SectionReader sectionReader(const char* fn) {

     SectionReader tmp(fn);

     return tmp;

   }

   /// \ingroup lemon_io

   ///

   /// \brief Reader for the contents of the \ref lgf-format "LGF" file

   ///

   /// This class can be used to read the sections, the map names and

   /// the attributes from a file. Usually, the LEMON programs know

   /// that, which type of graph, which maps and which attributes

   /// should be read from a file, but in general tools (like glemon)

   /// the contents of an LGF file should be guessed somehow. This class

   /// reads the graph and stores the appropriate information for

   /// reading the graph.

   ///

   ///\code

   /// LgfContents contents("graph.lgf");

   /// contents.run();

   ///

   /// // Does it contain any node section and arc section?

   /// if (contents.nodeSectionNum() == 0 || contents.arcSectionNum()) {

   ///   std::cerr << "Failure, cannot find graph." << std::endl;

   ///   return -1;

   /// }

   /// std::cout << "The name of the default node section: "

   ///           << contents.nodeSection(0) << std::endl;

   /// std::cout << "The number of the arc maps: "

   ///           << contents.arcMaps(0).size() << std::endl;

   /// std::cout << "The name of second arc map: "

   ///           << contents.arcMaps(0)[1] << std::endl;

   ///\endcode

   class LgfContents {

   private:

     std::istream* _is;

     bool local_is;

     std::vector<std::string> _node_sections;

     std::vector<std::string> _edge_sections;

     std::vector<std::string> _attribute_sections;

     std::vector<std::string> _extra_sections;

     std::vector<bool> _arc_sections;

     std::vector<std::vector<std::string> > _node_maps;

     std::vector<std::vector<std::string> > _edge_maps;

     std::vector<std::vector<std::string> > _attributes;

     int line_num;

     std::istringstream line;

   public:

     /// \brief Constructor

     ///

     /// Construct an \e LGF contents reader, which reads from the given

     /// input stream.

     LgfContents(std::istream& is)

       : _is(&is), local_is(false) {}

     /// \brief Constructor

     ///

     /// Construct an \e LGF contents reader, which reads from the given

     /// file.

     LgfContents(const std::string& fn)

       : _is(new std::ifstream(fn.c_str())), local_is(true) {

       if (!(*_is)) {

         delete _is;

         throw IoError("Cannot open file", fn);

       }

     }

     /// \brief Constructor

     ///

     /// Construct an \e LGF contents reader, which reads from the given

     /// file.

     LgfContents(const char* fn)

       : _is(new std::ifstream(fn)), local_is(true) {

       if (!(*_is)) {

         delete _is;

         throw IoError("Cannot open file", fn);

       }

     }

     /// \brief Destructor

     ~LgfContents() {

       if (local_is) delete _is;

     }

   private:

     LgfContents(const LgfContents&);

     LgfContents& operator=(const LgfContents&);

   public:

     /// \name Node sections

     /// @{

     /// \brief Gives back the number of node sections in the file.

     ///

     /// Gives back the number of node sections in the file.

     int nodeSectionNum() const {

       return _node_sections.size();

     }

     /// \brief Returns the node section name at the given position.

     ///

     /// Returns the node section name at the given position.

     const std::string& nodeSection(int i) const {

       return _node_sections[i];

     }

     /// \brief Gives back the node maps for the given section.

     ///

     /// Gives back the node maps for the given section.

     const std::vector<std::string>& nodeMapNames(int i) const {

       return _node_maps[i];

     }

     /// @}

     /// \name Arc/Edge sections

     /// @{

     /// \brief Gives back the number of arc/edge sections in the file.

     ///

     /// Gives back the number of arc/edge sections in the file.

     /// \note It is synonym of \c edgeSectionNum().

     int arcSectionNum() const {

       return _edge_sections.size();

     }

     /// \brief Returns the arc/edge section name at the given position.

     ///

     /// Returns the arc/edge section name at the given position.

     /// \note It is synonym of \c edgeSection().

     const std::string& arcSection(int i) const {

       return _edge_sections[i];

     }

     /// \brief Gives back the arc/edge maps for the given section.

     ///

     /// Gives back the arc/edge maps for the given section.

     /// \note It is synonym of \c edgeMapNames().

     const std::vector<std::string>& arcMapNames(int i) const {

       return _edge_maps[i];

     }

     /// @}

     /// \name Synonyms

     /// @{

     /// \brief Gives back the number of arc/edge sections in the file.

     ///

     /// Gives back the number of arc/edge sections in the file.

     /// \note It is synonym of \c arcSectionNum().

     int edgeSectionNum() const {

       return _edge_sections.size();

     }

     /// \brief Returns the section name at the given position.

     ///

     /// Returns the section name at the given position.

     /// \note It is synonym of \c arcSection().

     const std::string& edgeSection(int i) const {

       return _edge_sections[i];

     }

     /// \brief Gives back the edge maps for the given section.

     ///

     /// Gives back the edge maps for the given section.

     /// \note It is synonym of \c arcMapNames().

     const std::vector<std::string>& edgeMapNames(int i) const {

       return _edge_maps[i];

     }

     /// @}

     /// \name Attribute sections

     /// @{

     /// \brief Gives back the number of attribute sections in the file.

     ///

     /// Gives back the number of attribute sections in the file.

     int attributeSectionNum() const {

       return _attribute_sections.size();

     }

     /// \brief Returns the attribute section name at the given position.

     ///

     /// Returns the attribute section name at the given position.

     const std::string& attributeSectionNames(int i) const {

       return _attribute_sections[i];

     }

     /// \brief Gives back the attributes for the given section.

     ///

     /// Gives back the attributes for the given section.

     const std::vector<std::string>& attributes(int i) const {

       return _attributes[i];

     }

     /// @}

     /// \name Extra sections

     /// @{

     /// \brief Gives back the number of extra sections in the file.

     ///

     /// Gives back the number of extra sections in the file.

     int extraSectionNum() const {

       return _extra_sections.size();

     }

     /// \brief Returns the extra section type at the given position.

     ///

     /// Returns the section type at the given position.

     const std::string& extraSection(int i) const {

       return _extra_sections[i];

     }

     /// @}

   private:

     bool readLine() {

       std::string str;

       while(++line_num, std::getline(*_is, str)) {

         line.clear(); line.str(str);

         char c;

         if (line >> std::ws >> c && c != '#') {

           line.putback(c);

           return true;

         }

       }

       return false;

     }

     bool readSuccess() {

       return static_cast<bool>(*_is);

     }

     void skipSection() {

       char c;

       while (readSuccess() && line >> c && c != '@') {

         readLine();

       }

       line.putback(c);

     }

     void readMaps(std::vector<std::string>& maps) {

       char c;

       if (!readLine() || !(line >> c) || c == '@') {

         if (readSuccess() && line) line.putback(c);

         return;

       }

       line.putback(c);

       std::string map;

       while (_reader_bits::readToken(line, map)) {

         maps.push_back(map);

       }

     }

     void readAttributes(std::vector<std::string>& attrs) {

       readLine();

       char c;

       while (readSuccess() && line >> c && c != '@') {

         line.putback(c);

         std::string attr;

         _reader_bits::readToken(line, attr);

         attrs.push_back(attr);

         readLine();

       }

       line.putback(c);

     }

   public:

     /// \name Execution of the contents reader

     /// @{

     /// \brief Starts the reading

     ///

     /// This function starts the reading.

     void run() {

       readLine();

       skipSection();

       while (readSuccess()) {

         char c;

         line >> c;

         std::string section, caption;

         _reader_bits::readToken(line, section);

         _reader_bits::readToken(line, caption);

         if (section == "nodes") {

           _node_sections.push_back(caption);

           _node_maps.push_back(std::vector<std::string>());

           readMaps(_node_maps.back());

           readLine(); skipSection();

         } else if (section == "arcs" || section == "edges") {

           _edge_sections.push_back(caption);

           _arc_sections.push_back(section == "arcs");

           _edge_maps.push_back(std::vector<std::string>());

           readMaps(_edge_maps.back());