COIN-OR::LEMON - Graph Library

Changeset 57:c1acf0018c0a in lemon for lemon


Ignore:
Timestamp:
01/20/08 20:43:48 (16 years ago)
Author:
Balazs Dezso <deba@…>
Branch:
default
Phase:
public
Message:

Port ListDigraph? and ListGraph? from svn -r 3433
Details:

  • port Digraph and Graph concepts
  • port ListDigraph? and ListGraph?
  • port Basic graph constructing tools
  • port Digraph and Graph tests
Location:
lemon
Files:
11 added
2 edited

Legend:

Unmodified
Added
Removed
  • lemon/Makefile.am

    r32 r57  
    2222
    2323bits_HEADERS += \
     24        lemon/bits/alteration_notifier.h \
     25        lemon/bits/array_map.h \
     26        lemon/bits/base_extender.h \
     27        lemon/bits/default_map.h \
    2428        lemon/bits/invalid.h \
    25         lemon/bits/utility.h
     29        lemon/bits/map_extender.h \
     30        lemon/bits/utility.h \
     31        lemon/bits/vector_map.h
    2632
    2733concept_HEADERS +=
     34        lemon/concept_check.h \
     35        lemon/concepts/digraph.h \
     36        lemon/concepts/graph.h \
     37        lemon/concepts/graph_components.h
  • lemon/list_graph.h

    r39 r57  
    1717 */
    1818
     19#ifndef LEMON_LIST_GRAPH_H
     20#define LEMON_LIST_GRAPH_H
     21
     22///\ingroup graphs
     23///\file
     24///\brief ListDigraph, ListGraph classes.
     25
     26#include <lemon/bits/graph_extender.h>
     27
     28#include <vector>
     29#include <list>
     30
     31namespace lemon {
     32
     33  class ListDigraphBase {
     34
     35  protected:
     36    struct NodeT {
     37      int first_in, first_out;
     38      int prev, next;
     39    };
     40 
     41    struct ArcT {
     42      int target, source;
     43      int prev_in, prev_out;
     44      int next_in, next_out;
     45    };
     46
     47    std::vector<NodeT> nodes;
     48
     49    int first_node;
     50
     51    int first_free_node;
     52
     53    std::vector<ArcT> arcs;
     54
     55    int first_free_arc;
     56   
     57  public:
     58   
     59    typedef ListDigraphBase Digraph;
     60   
     61    class Node {
     62      friend class ListDigraphBase;
     63    protected:
     64
     65      int id;
     66      explicit Node(int pid) { id = pid;}
     67
     68    public:
     69      Node() {}
     70      Node (Invalid) { id = -1; }
     71      bool operator==(const Node& node) const {return id == node.id;}
     72      bool operator!=(const Node& node) const {return id != node.id;}
     73      bool operator<(const Node& node) const {return id < node.id;}
     74    };
     75
     76    class Arc {
     77      friend class ListDigraphBase;
     78    protected:
     79
     80      int id;
     81      explicit Arc(int pid) { id = pid;}
     82
     83    public:
     84      Arc() {}
     85      Arc (Invalid) { id = -1; }
     86      bool operator==(const Arc& arc) const {return id == arc.id;}
     87      bool operator!=(const Arc& arc) const {return id != arc.id;}
     88      bool operator<(const Arc& arc) const {return id < arc.id;}
     89    };
     90
     91
     92
     93    ListDigraphBase()
     94      : nodes(), first_node(-1),
     95        first_free_node(-1), arcs(), first_free_arc(-1) {}
     96
     97   
     98    int maxNodeId() const { return nodes.size()-1; }
     99    int maxArcId() const { return arcs.size()-1; }
     100
     101    Node source(Arc e) const { return Node(arcs[e.id].source); }
     102    Node target(Arc e) const { return Node(arcs[e.id].target); }
     103
     104
     105    void first(Node& node) const {
     106      node.id = first_node;
     107    }
     108
     109    void next(Node& node) const {
     110      node.id = nodes[node.id].next;
     111    }
     112
     113
     114    void first(Arc& e) const {
     115      int n;
     116      for(n = first_node;
     117          n!=-1 && nodes[n].first_in == -1;
     118          n = nodes[n].next);
     119      e.id = (n == -1) ? -1 : nodes[n].first_in;
     120    }
     121
     122    void next(Arc& arc) const {
     123      if (arcs[arc.id].next_in != -1) {
     124        arc.id = arcs[arc.id].next_in;
     125      } else {
     126        int n;
     127        for(n = nodes[arcs[arc.id].target].next;
     128          n!=-1 && nodes[n].first_in == -1;
     129          n = nodes[n].next);
     130        arc.id = (n == -1) ? -1 : nodes[n].first_in;
     131      }     
     132    }
     133
     134    void firstOut(Arc &e, const Node& v) const {
     135      e.id = nodes[v.id].first_out;
     136    }
     137    void nextOut(Arc &e) const {
     138      e.id=arcs[e.id].next_out;
     139    }
     140
     141    void firstIn(Arc &e, const Node& v) const {
     142      e.id = nodes[v.id].first_in;
     143    }
     144    void nextIn(Arc &e) const {
     145      e.id=arcs[e.id].next_in;
     146    }
     147
     148   
     149    static int id(Node v) { return v.id; }
     150    static int id(Arc e) { return e.id; }
     151
     152    static Node nodeFromId(int id) { return Node(id);}
     153    static Arc arcFromId(int id) { return Arc(id);}
     154
     155    Node addNode() {     
     156      int n;
     157     
     158      if(first_free_node==-1) {
     159        n = nodes.size();
     160        nodes.push_back(NodeT());
     161      } else {
     162        n = first_free_node;
     163        first_free_node = nodes[n].next;
     164      }
     165     
     166      nodes[n].next = first_node;
     167      if(first_node != -1) nodes[first_node].prev = n;
     168      first_node = n;
     169      nodes[n].prev = -1;
     170     
     171      nodes[n].first_in = nodes[n].first_out = -1;
     172     
     173      return Node(n);
     174    }
     175   
     176    Arc addArc(Node u, Node v) {
     177      int n;     
     178
     179      if (first_free_arc == -1) {
     180        n = arcs.size();
     181        arcs.push_back(ArcT());
     182      } else {
     183        n = first_free_arc;
     184        first_free_arc = arcs[n].next_in;
     185      }
     186     
     187      arcs[n].source = u.id;
     188      arcs[n].target = v.id;
     189
     190      arcs[n].next_out = nodes[u.id].first_out;
     191      if(nodes[u.id].first_out != -1) {
     192        arcs[nodes[u.id].first_out].prev_out = n;
     193      }
     194     
     195      arcs[n].next_in = nodes[v.id].first_in;
     196      if(nodes[v.id].first_in != -1) {
     197        arcs[nodes[v.id].first_in].prev_in = n;
     198      }
     199     
     200      arcs[n].prev_in = arcs[n].prev_out = -1;
     201       
     202      nodes[u.id].first_out = nodes[v.id].first_in = n;
     203
     204      return Arc(n);
     205    }
     206   
     207    void erase(const Node& node) {
     208      int n = node.id;
     209     
     210      if(nodes[n].next != -1) {
     211        nodes[nodes[n].next].prev = nodes[n].prev;
     212      }
     213     
     214      if(nodes[n].prev != -1) {
     215        nodes[nodes[n].prev].next = nodes[n].next;
     216      } else {
     217        first_node = nodes[n].next;
     218      }
     219     
     220      nodes[n].next = first_free_node;
     221      first_free_node = n;
     222
     223    }
     224   
     225    void erase(const Arc& arc) {
     226      int n = arc.id;
     227     
     228      if(arcs[n].next_in!=-1) {
     229        arcs[arcs[n].next_in].prev_in = arcs[n].prev_in;
     230      }
     231
     232      if(arcs[n].prev_in!=-1) {
     233        arcs[arcs[n].prev_in].next_in = arcs[n].next_in;
     234      } else {
     235        nodes[arcs[n].target].first_in = arcs[n].next_in;
     236      }
     237
     238     
     239      if(arcs[n].next_out!=-1) {
     240        arcs[arcs[n].next_out].prev_out = arcs[n].prev_out;
     241      }
     242
     243      if(arcs[n].prev_out!=-1) {
     244        arcs[arcs[n].prev_out].next_out = arcs[n].next_out;
     245      } else {
     246        nodes[arcs[n].source].first_out = arcs[n].next_out;
     247      }
     248     
     249      arcs[n].next_in = first_free_arc;
     250      first_free_arc = n;     
     251
     252    }
     253
     254    void clear() {
     255      arcs.clear();
     256      nodes.clear();
     257      first_node = first_free_node = first_free_arc = -1;
     258    }
     259
     260  protected:
     261    void changeTarget(Arc e, Node n)
     262    {
     263      if(arcs[e.id].next_in != -1)
     264        arcs[arcs[e.id].next_in].prev_in = arcs[e.id].prev_in;
     265      if(arcs[e.id].prev_in != -1)
     266        arcs[arcs[e.id].prev_in].next_in = arcs[e.id].next_in;
     267      else nodes[arcs[e.id].target].first_in = arcs[e.id].next_in;
     268      if (nodes[n.id].first_in != -1) {
     269        arcs[nodes[n.id].first_in].prev_in = e.id;
     270      }
     271      arcs[e.id].target = n.id;
     272      arcs[e.id].prev_in = -1;
     273      arcs[e.id].next_in = nodes[n.id].first_in;
     274      nodes[n.id].first_in = e.id;
     275    }
     276    void changeSource(Arc e, Node n)
     277    {
     278      if(arcs[e.id].next_out != -1)
     279        arcs[arcs[e.id].next_out].prev_out = arcs[e.id].prev_out;
     280      if(arcs[e.id].prev_out != -1)
     281        arcs[arcs[e.id].prev_out].next_out = arcs[e.id].next_out;
     282      else nodes[arcs[e.id].source].first_out = arcs[e.id].next_out;
     283      if (nodes[n.id].first_out != -1) {
     284        arcs[nodes[n.id].first_out].prev_out = e.id;
     285      }
     286      arcs[e.id].source = n.id;
     287      arcs[e.id].prev_out = -1;
     288      arcs[e.id].next_out = nodes[n.id].first_out;
     289      nodes[n.id].first_out = e.id;
     290    }
     291
     292  };
     293
     294  typedef DigraphExtender<ListDigraphBase> ExtendedListDigraphBase;
     295
     296  /// \addtogroup digraphs
     297  /// @{
     298
     299  ///A list digraph class.
     300
     301  ///This is a simple and fast digraph implementation.
     302  ///
     303  ///It conforms to the \ref concepts::Digraph "Digraph concept" and it
     304  ///also provides several additional useful extra functionalities.
     305  ///The most of the member functions and nested classes are
     306  ///documented only in the concept class.
     307  ///
     308  ///An important extra feature of this digraph implementation is that
     309  ///its maps are real \ref concepts::ReferenceMap "reference map"s.
     310  ///
     311  ///\sa concepts::Digraph.
     312
     313  class ListDigraph : public ExtendedListDigraphBase {
     314  private:
     315    ///ListDigraph is \e not copy constructible. Use DigraphCopy() instead.
     316   
     317    ///ListDigraph is \e not copy constructible. Use DigraphCopy() instead.
     318    ///
     319    ListDigraph(const ListDigraph &) :ExtendedListDigraphBase() {};
     320    ///\brief Assignment of ListDigraph to another one is \e not allowed.
     321    ///Use DigraphCopy() instead.
     322
     323    ///Assignment of ListDigraph to another one is \e not allowed.
     324    ///Use DigraphCopy() instead.
     325    void operator=(const ListDigraph &) {}
     326  public:
     327
     328    typedef ExtendedListDigraphBase Parent;
     329
     330    /// Constructor
     331   
     332    /// Constructor.
     333    ///
     334    ListDigraph() {}
     335
     336    ///Add a new node to the digraph.
     337   
     338    /// \return the new node.
     339    ///
     340    Node addNode() { return Parent::addNode(); }
     341
     342    ///Add a new arc to the digraph.
     343   
     344    ///Add a new arc to the digraph with source node \c s
     345    ///and target node \c t.
     346    ///\return the new arc.
     347    Arc addArc(const Node& s, const Node& t) {
     348      return Parent::addArc(s, t);
     349    }
     350
     351    /// Changes the target of \c e to \c n
     352
     353    /// Changes the target of \c e to \c n
     354    ///
     355    ///\note The <tt>ArcIt</tt>s and <tt>OutArcIt</tt>s referencing
     356    ///the changed arc remain valid. However <tt>InArcIt</tt>s are
     357    ///invalidated.
     358    ///\warning This functionality cannot be used together with the Snapshot
     359    ///feature.
     360    void changeTarget(Arc e, Node n) {
     361      Parent::changeTarget(e,n);
     362    }
     363    /// Changes the source of \c e to \c n
     364
     365    /// Changes the source of \c e to \c n
     366    ///
     367    ///\note The <tt>ArcIt</tt>s and <tt>InArcIt</tt>s referencing
     368    ///the changed arc remain valid. However <tt>OutArcIt</tt>s are
     369    ///invalidated.
     370    ///\warning This functionality cannot be used together with the Snapshot
     371    ///feature.
     372    void changeSource(Arc e, Node n) {
     373      Parent::changeSource(e,n);
     374    }
     375
     376    /// Invert the direction of an arc.
     377
     378    ///\note The <tt>ArcIt</tt>s referencing the changed arc remain
     379    ///valid. However <tt>OutArcIt</tt>s and <tt>InArcIt</tt>s are
     380    ///invalidated.
     381    ///\warning This functionality cannot be used together with the Snapshot
     382    ///feature.
     383    void reverseArc(Arc e) {
     384      Node t=target(e);
     385      changeTarget(e,source(e));
     386      changeSource(e,t);
     387    }
     388
     389    /// Using this it is possible to avoid the superfluous memory
     390    /// allocation: if you know that the digraph you want to build will
     391    /// be very large (e.g. it will contain millions of nodes and/or arcs)
     392    /// then it is worth reserving space for this amount before starting
     393    /// to build the digraph.
     394    /// \sa reserveArc
     395    void reserveNode(int n) { nodes.reserve(n); };
     396
     397    /// \brief Using this it is possible to avoid the superfluous memory
     398    /// allocation.
     399
     400    /// Using this it is possible to avoid the superfluous memory
     401    /// allocation: if you know that the digraph you want to build will
     402    /// be very large (e.g. it will contain millions of nodes and/or arcs)
     403    /// then it is worth reserving space for this amount before starting
     404    /// to build the digraph.
     405    /// \sa reserveNode
     406    void reserveArc(int m) { arcs.reserve(m); };
     407
     408    ///Contract two nodes.
     409
     410    ///This function contracts two nodes.
     411    ///
     412    ///Node \p b will be removed but instead of deleting
     413    ///incident arcs, they will be joined to \p a.
     414    ///The last parameter \p r controls whether to remove loops. \c true
     415    ///means that loops will be removed.
     416    ///
     417    ///\note The <tt>ArcIt</tt>s
     418    ///referencing a moved arc remain
     419    ///valid. However <tt>InArcIt</tt>s and <tt>OutArcIt</tt>s
     420    ///may be invalidated.
     421    ///\warning This functionality cannot be used together with the Snapshot
     422    ///feature.
     423    void contract(Node a, Node b, bool r = true)
     424    {
     425      for(OutArcIt e(*this,b);e!=INVALID;) {
     426        OutArcIt f=e;
     427        ++f;
     428        if(r && target(e)==a) erase(e);
     429        else changeSource(e,a);
     430        e=f;
     431      }
     432      for(InArcIt e(*this,b);e!=INVALID;) {
     433        InArcIt f=e;
     434        ++f;
     435        if(r && source(e)==a) erase(e);
     436        else changeTarget(e,a);
     437        e=f;
     438      }
     439      erase(b);
     440    }
     441
     442    ///Split a node.
     443
     444    ///This function splits a node. First a new node is added to the digraph,
     445    ///then the source of each outgoing arc of \c n is moved to this new node.
     446    ///If \c connect is \c true (this is the default value), then a new arc
     447    ///from \c n to the newly created node is also added.
     448    ///\return The newly created node.
     449    ///
     450    ///\note The <tt>ArcIt</tt>s referencing a moved arc remain
     451    ///valid. However <tt>InArcIt</tt>s and <tt>OutArcIt</tt>s may
     452    ///be invalidated. 
     453    ///
     454    ///\warning This functionality cannot be used together with the
     455    ///Snapshot feature.  \todo It could be implemented in a bit
     456    ///faster way.
     457    Node split(Node n, bool connect = true) {
     458      Node b = addNode();
     459      for(OutArcIt e(*this,n);e!=INVALID;) {
     460        OutArcIt f=e;
     461        ++f;
     462        changeSource(e,b);
     463        e=f;
     464      }
     465      if (connect) addArc(n,b);
     466      return b;
     467    }
     468     
     469    ///Split an arc.
     470
     471    ///This function splits an arc. First a new node \c b is added to
     472    ///the digraph, then the original arc is re-targeted to \c
     473    ///b. Finally an arc from \c b to the original target is added.
     474    ///\return The newly created node. 
     475    ///\warning This functionality
     476    ///cannot be used together with the Snapshot feature.
     477    Node split(Arc e) {
     478      Node b = addNode();
     479      addArc(b,target(e));
     480      changeTarget(e,b);
     481      return b;
     482    }
     483     
     484    /// \brief Class to make a snapshot of the digraph and restore
     485    /// to it later.
     486    ///
     487    /// Class to make a snapshot of the digraph and to restore it
     488    /// later.
     489    ///
     490    /// The newly added nodes and arcs can be removed using the
     491    /// restore() function.
     492    ///
     493    /// \warning Arc and node deletions cannot be restored. This
     494    /// events invalidate the snapshot.
     495    class Snapshot {
     496    protected:
     497
     498      typedef Parent::NodeNotifier NodeNotifier;
     499
     500      class NodeObserverProxy : public NodeNotifier::ObserverBase {
     501      public:
     502
     503        NodeObserverProxy(Snapshot& _snapshot)
     504          : snapshot(_snapshot) {}
     505
     506        using NodeNotifier::ObserverBase::attach;
     507        using NodeNotifier::ObserverBase::detach;
     508        using NodeNotifier::ObserverBase::attached;
     509       
     510      protected:
     511       
     512        virtual void add(const Node& node) {
     513          snapshot.addNode(node);
     514        }
     515        virtual void add(const std::vector<Node>& nodes) {
     516          for (int i = nodes.size() - 1; i >= 0; ++i) {
     517            snapshot.addNode(nodes[i]);
     518          }
     519        }
     520        virtual void erase(const Node& node) {
     521          snapshot.eraseNode(node);
     522        }
     523        virtual void erase(const std::vector<Node>& nodes) {
     524          for (int i = 0; i < int(nodes.size()); ++i) {
     525            snapshot.eraseNode(nodes[i]);
     526          }
     527        }
     528        virtual void build() {
     529          Node node;
     530          std::vector<Node> nodes;
     531          for (notifier()->first(node); node != INVALID;
     532               notifier()->next(node)) {
     533            nodes.push_back(node);
     534          }
     535          for (int i = nodes.size() - 1; i >= 0; --i) {
     536            snapshot.addNode(nodes[i]);
     537          }
     538        }
     539        virtual void clear() {
     540          Node node;
     541          for (notifier()->first(node); node != INVALID;
     542               notifier()->next(node)) {
     543            snapshot.eraseNode(node);
     544          }
     545        }
     546
     547        Snapshot& snapshot;
     548      };
     549
     550      class ArcObserverProxy : public ArcNotifier::ObserverBase {
     551      public:
     552
     553        ArcObserverProxy(Snapshot& _snapshot)
     554          : snapshot(_snapshot) {}
     555
     556        using ArcNotifier::ObserverBase::attach;
     557        using ArcNotifier::ObserverBase::detach;
     558        using ArcNotifier::ObserverBase::attached;
     559       
     560      protected:
     561
     562        virtual void add(const Arc& arc) {
     563          snapshot.addArc(arc);
     564        }
     565        virtual void add(const std::vector<Arc>& arcs) {
     566          for (int i = arcs.size() - 1; i >= 0; ++i) {
     567            snapshot.addArc(arcs[i]);
     568          }
     569        }
     570        virtual void erase(const Arc& arc) {
     571          snapshot.eraseArc(arc);
     572        }
     573        virtual void erase(const std::vector<Arc>& arcs) {
     574          for (int i = 0; i < int(arcs.size()); ++i) {
     575            snapshot.eraseArc(arcs[i]);
     576          }
     577        }
     578        virtual void build() {
     579          Arc arc;
     580          std::vector<Arc> arcs;
     581          for (notifier()->first(arc); arc != INVALID;
     582               notifier()->next(arc)) {
     583            arcs.push_back(arc);
     584          }
     585          for (int i = arcs.size() - 1; i >= 0; --i) {
     586            snapshot.addArc(arcs[i]);
     587          }
     588        }
     589        virtual void clear() {
     590          Arc arc;
     591          for (notifier()->first(arc); arc != INVALID;
     592               notifier()->next(arc)) {
     593            snapshot.eraseArc(arc);
     594          }
     595        }
     596
     597        Snapshot& snapshot;
     598      };
     599     
     600      ListDigraph *digraph;
     601
     602      NodeObserverProxy node_observer_proxy;
     603      ArcObserverProxy arc_observer_proxy;
     604
     605      std::list<Node> added_nodes;
     606      std::list<Arc> added_arcs;
     607
     608
     609      void addNode(const Node& node) {
     610        added_nodes.push_front(node);       
     611      }
     612      void eraseNode(const Node& node) {
     613        std::list<Node>::iterator it =
     614          std::find(added_nodes.begin(), added_nodes.end(), node);
     615        if (it == added_nodes.end()) {
     616          clear();
     617          arc_observer_proxy.detach();
     618          throw NodeNotifier::ImmediateDetach();
     619        } else {
     620          added_nodes.erase(it);
     621        }
     622      }
     623
     624      void addArc(const Arc& arc) {
     625        added_arcs.push_front(arc);       
     626      }
     627      void eraseArc(const Arc& arc) {
     628        std::list<Arc>::iterator it =
     629          std::find(added_arcs.begin(), added_arcs.end(), arc);
     630        if (it == added_arcs.end()) {
     631          clear();
     632          node_observer_proxy.detach();
     633          throw ArcNotifier::ImmediateDetach();
     634        } else {
     635          added_arcs.erase(it);
     636        }       
     637      }
     638
     639      void attach(ListDigraph &_digraph) {
     640        digraph = &_digraph;
     641        node_observer_proxy.attach(digraph->notifier(Node()));
     642        arc_observer_proxy.attach(digraph->notifier(Arc()));
     643      }
     644           
     645      void detach() {
     646        node_observer_proxy.detach();
     647        arc_observer_proxy.detach();
     648      }
     649
     650      bool attached() const {
     651        return node_observer_proxy.attached();
     652      }
     653
     654      void clear() {
     655        added_nodes.clear();
     656        added_arcs.clear();       
     657      }
     658
     659    public:
     660
     661      /// \brief Default constructor.
     662      ///
     663      /// Default constructor.
     664      /// To actually make a snapshot you must call save().
     665      Snapshot()
     666        : digraph(0), node_observer_proxy(*this),
     667          arc_observer_proxy(*this) {}
     668     
     669      /// \brief Constructor that immediately makes a snapshot.
     670      ///     
     671      /// This constructor immediately makes a snapshot of the digraph.
     672      /// \param _digraph The digraph we make a snapshot of.
     673      Snapshot(ListDigraph &_digraph)
     674        : node_observer_proxy(*this),
     675          arc_observer_proxy(*this) {
     676        attach(_digraph);
     677      }
     678     
     679      /// \brief Make a snapshot.
     680      ///
     681      /// Make a snapshot of the digraph.
     682      ///
     683      /// This function can be called more than once. In case of a repeated
     684      /// call, the previous snapshot gets lost.
     685      /// \param _digraph The digraph we make the snapshot of.
     686      void save(ListDigraph &_digraph) {
     687        if (attached()) {
     688          detach();
     689          clear();
     690        }
     691        attach(_digraph);
     692      }
     693     
     694      /// \brief Undo the changes until the last snapshot.
     695      //
     696      /// Undo the changes until the last snapshot created by save().
     697      void restore() {
     698        detach();
     699        for(std::list<Arc>::iterator it = added_arcs.begin();
     700            it != added_arcs.end(); ++it) {
     701          digraph->erase(*it);
     702        }
     703        for(std::list<Node>::iterator it = added_nodes.begin();
     704            it != added_nodes.end(); ++it) {
     705          digraph->erase(*it);
     706        }
     707        clear();
     708      }
     709
     710      /// \brief Gives back true when the snapshot is valid.
     711      ///
     712      /// Gives back true when the snapshot is valid.
     713      bool valid() const {
     714        return attached();
     715      }
     716    };
     717   
     718  };
     719
     720  ///@}
     721
     722  class ListGraphBase {
     723
     724  protected:
     725
     726    struct NodeT {
     727      int first_out;
     728      int prev, next;
     729    };
     730 
     731    struct ArcT {
     732      int target;
     733      int prev_out, next_out;
     734    };
     735
     736    std::vector<NodeT> nodes;
     737
     738    int first_node;
     739
     740    int first_free_node;
     741
     742    std::vector<ArcT> arcs;
     743
     744    int first_free_arc;
     745   
     746  public:
     747   
     748    typedef ListGraphBase Digraph;
     749
     750    class Node;
     751    class Arc;
     752    class Edge;
     753   
     754    class Node {
     755      friend class ListGraphBase;
     756    protected:
     757
     758      int id;
     759      explicit Node(int pid) { id = pid;}
     760
     761    public:
     762      Node() {}
     763      Node (Invalid) { id = -1; }
     764      bool operator==(const Node& node) const {return id == node.id;}
     765      bool operator!=(const Node& node) const {return id != node.id;}
     766      bool operator<(const Node& node) const {return id < node.id;}
     767    };
     768
     769    class Edge {
     770      friend class ListGraphBase;
     771    protected:
     772
     773      int id;
     774      explicit Edge(int pid) { id = pid;}
     775
     776    public:
     777      Edge() {}
     778      Edge (Invalid) { id = -1; }
     779      bool operator==(const Edge& arc) const {return id == arc.id;}
     780      bool operator!=(const Edge& arc) const {return id != arc.id;}
     781      bool operator<(const Edge& arc) const {return id < arc.id;}
     782    };
     783
     784    class Arc {
     785      friend class ListGraphBase;
     786    protected:
     787
     788      int id;
     789      explicit Arc(int pid) { id = pid;}
     790
     791    public:
     792      operator Edge() const { return edgeFromId(id / 2); }
     793
     794      Arc() {}
     795      Arc (Invalid) { id = -1; }
     796      bool operator==(const Arc& arc) const {return id == arc.id;}
     797      bool operator!=(const Arc& arc) const {return id != arc.id;}
     798      bool operator<(const Arc& arc) const {return id < arc.id;}
     799    };
     800
     801
     802
     803    ListGraphBase()
     804      : nodes(), first_node(-1),
     805        first_free_node(-1), arcs(), first_free_arc(-1) {}
     806
     807   
     808    int maxNodeId() const { return nodes.size()-1; }
     809    int maxEdgeId() const { return arcs.size() / 2 - 1; }
     810    int maxArcId() const { return arcs.size()-1; }
     811
     812    Node source(Arc e) const { return Node(arcs[e.id ^ 1].target); }
     813    Node target(Arc e) const { return Node(arcs[e.id].target); }
     814
     815    Node u(Edge e) const { return Node(arcs[2 * e.id].target); }
     816    Node v(Edge e) const { return Node(arcs[2 * e.id + 1].target); }
     817
     818    static bool direction(Arc e) {
     819      return (e.id & 1) == 1;
     820    }
     821
     822    static Arc direct(Edge e, bool d) {
     823      return Arc(e.id * 2 + (d ? 1 : 0));
     824    }
     825
     826    void first(Node& node) const {
     827      node.id = first_node;
     828    }
     829
     830    void next(Node& node) const {
     831      node.id = nodes[node.id].next;
     832    }
     833
     834    void first(Arc& e) const {
     835      int n = first_node;
     836      while (n != -1 && nodes[n].first_out == -1) {
     837        n = nodes[n].next;
     838      }
     839      e.id = (n == -1) ? -1 : nodes[n].first_out;
     840    }
     841
     842    void next(Arc& e) const {
     843      if (arcs[e.id].next_out != -1) {
     844        e.id = arcs[e.id].next_out;
     845      } else {
     846        int n = nodes[arcs[e.id ^ 1].target].next;
     847        while(n != -1 && nodes[n].first_out == -1) {
     848          n = nodes[n].next;
     849        }
     850        e.id = (n == -1) ? -1 : nodes[n].first_out;
     851      }     
     852    }
     853
     854    void first(Edge& e) const {
     855      int n = first_node;
     856      while (n != -1) {
     857        e.id = nodes[n].first_out;
     858        while ((e.id & 1) != 1) {
     859          e.id = arcs[e.id].next_out;
     860        }
     861        if (e.id != -1) {
     862          e.id /= 2;
     863          return;
     864        }
     865        n = nodes[n].next;
     866      }
     867      e.id = -1;
     868    }
     869
     870    void next(Edge& e) const {
     871      int n = arcs[e.id * 2].target;
     872      e.id = arcs[(e.id * 2) | 1].next_out;
     873      while ((e.id & 1) != 1) {
     874        e.id = arcs[e.id].next_out;
     875      }
     876      if (e.id != -1) {
     877        e.id /= 2;
     878        return;
     879      }
     880      n = nodes[n].next;
     881      while (n != -1) {
     882        e.id = nodes[n].first_out;
     883        while ((e.id & 1) != 1) {
     884          e.id = arcs[e.id].next_out;
     885        }
     886        if (e.id != -1) {
     887          e.id /= 2;
     888          return;
     889        }
     890        n = nodes[n].next;
     891      }
     892      e.id = -1;
     893    }
     894
     895    void firstOut(Arc &e, const Node& v) const {
     896      e.id = nodes[v.id].first_out;
     897    }
     898    void nextOut(Arc &e) const {
     899      e.id = arcs[e.id].next_out;
     900    }
     901
     902    void firstIn(Arc &e, const Node& v) const {
     903      e.id = ((nodes[v.id].first_out) ^ 1);
     904      if (e.id == -2) e.id = -1;
     905    }
     906    void nextIn(Arc &e) const {
     907      e.id = ((arcs[e.id ^ 1].next_out) ^ 1);
     908      if (e.id == -2) e.id = -1;
     909    }
     910
     911    void firstInc(Edge &e, bool& d, const Node& v) const {
     912      int de = nodes[v.id].first_out;
     913      if (de != -1 ) {
     914        e.id = de / 2;
     915        d = ((de & 1) == 1);
     916      } else {
     917        e.id = -1;
     918        d = true;
     919      }
     920    }
     921    void nextInc(Edge &e, bool& d) const {
     922      int de = (arcs[(e.id * 2) | (d ? 1 : 0)].next_out);
     923      if (de != -1 ) {
     924        e.id = de / 2;
     925        d = ((de & 1) == 1);
     926      } else {
     927        e.id = -1;
     928        d = true;
     929      }
     930    }
     931   
     932    static int id(Node v) { return v.id; }
     933    static int id(Arc e) { return e.id; }
     934    static int id(Edge e) { return e.id; }
     935
     936    static Node nodeFromId(int id) { return Node(id);}
     937    static Arc arcFromId(int id) { return Arc(id);}
     938    static Edge edgeFromId(int id) { return Edge(id);}
     939
     940    Node addNode() {     
     941      int n;
     942     
     943      if(first_free_node==-1) {
     944        n = nodes.size();
     945        nodes.push_back(NodeT());
     946      } else {
     947        n = first_free_node;
     948        first_free_node = nodes[n].next;
     949      }
     950     
     951      nodes[n].next = first_node;
     952      if (first_node != -1) nodes[first_node].prev = n;
     953      first_node = n;
     954      nodes[n].prev = -1;
     955     
     956      nodes[n].first_out = -1;
     957     
     958      return Node(n);
     959    }
     960   
     961    Edge addEdge(Node u, Node v) {
     962      int n;     
     963
     964      if (first_free_arc == -1) {
     965        n = arcs.size();
     966        arcs.push_back(ArcT());
     967        arcs.push_back(ArcT());
     968      } else {
     969        n = first_free_arc;
     970        first_free_arc = arcs[n].next_out;
     971      }
     972     
     973      arcs[n].target = u.id;
     974      arcs[n | 1].target = v.id;
     975
     976      arcs[n].next_out = nodes[v.id].first_out;
     977      if (nodes[v.id].first_out != -1) {
     978        arcs[nodes[v.id].first_out].prev_out = n;
     979      }     
     980      arcs[n].prev_out = -1;
     981      nodes[v.id].first_out = n;
     982     
     983      arcs[n | 1].next_out = nodes[u.id].first_out;
     984      if (nodes[u.id].first_out != -1) {
     985        arcs[nodes[u.id].first_out].prev_out = (n | 1);
     986      }
     987      arcs[n | 1].prev_out = -1;     
     988      nodes[u.id].first_out = (n | 1);
     989
     990      return Edge(n / 2);
     991    }
     992   
     993    void erase(const Node& node) {
     994      int n = node.id;
     995     
     996      if(nodes[n].next != -1) {
     997        nodes[nodes[n].next].prev = nodes[n].prev;
     998      }
     999     
     1000      if(nodes[n].prev != -1) {
     1001        nodes[nodes[n].prev].next = nodes[n].next;
     1002      } else {
     1003        first_node = nodes[n].next;
     1004      }
     1005     
     1006      nodes[n].next = first_free_node;
     1007      first_free_node = n;
     1008
     1009    }
     1010   
     1011    void erase(const Edge& arc) {
     1012      int n = arc.id * 2;
     1013     
     1014      if (arcs[n].next_out != -1) {
     1015        arcs[arcs[n].next_out].prev_out = arcs[n].prev_out;
     1016      }
     1017
     1018      if (arcs[n].prev_out != -1) {
     1019        arcs[arcs[n].prev_out].next_out = arcs[n].next_out;
     1020      } else {
     1021        nodes[arcs[n | 1].target].first_out = arcs[n].next_out;
     1022      }
     1023
     1024      if (arcs[n | 1].next_out != -1) {
     1025        arcs[arcs[n | 1].next_out].prev_out = arcs[n | 1].prev_out;
     1026      }
     1027
     1028      if (arcs[n | 1].prev_out != -1) {
     1029        arcs[arcs[n | 1].prev_out].next_out = arcs[n | 1].next_out;
     1030      } else {
     1031        nodes[arcs[n].target].first_out = arcs[n | 1].next_out;
     1032      }
     1033     
     1034      arcs[n].next_out = first_free_arc;
     1035      first_free_arc = n;     
     1036
     1037    }
     1038
     1039    void clear() {
     1040      arcs.clear();
     1041      nodes.clear();
     1042      first_node = first_free_node = first_free_arc = -1;
     1043    }
     1044
     1045  protected:
     1046
     1047    void changeTarget(Edge e, Node n) {
     1048      if(arcs[2 * e.id].next_out != -1) {
     1049        arcs[arcs[2 * e.id].next_out].prev_out = arcs[2 * e.id].prev_out;
     1050      }
     1051      if(arcs[2 * e.id].prev_out != -1) {
     1052        arcs[arcs[2 * e.id].prev_out].next_out =
     1053          arcs[2 * e.id].next_out;
     1054      } else {
     1055        nodes[arcs[(2 * e.id) | 1].target].first_out =
     1056          arcs[2 * e.id].next_out;
     1057      }
     1058
     1059      if (nodes[n.id].first_out != -1) {
     1060        arcs[nodes[n.id].first_out].prev_out = 2 * e.id;
     1061      }
     1062      arcs[(2 * e.id) | 1].target = n.id;
     1063      arcs[2 * e.id].prev_out = -1;
     1064      arcs[2 * e.id].next_out = nodes[n.id].first_out;
     1065      nodes[n.id].first_out = 2 * e.id;
     1066    }
     1067
     1068    void changeSource(Edge e, Node n) {
     1069      if(arcs[(2 * e.id) | 1].next_out != -1) {
     1070        arcs[arcs[(2 * e.id) | 1].next_out].prev_out =
     1071          arcs[(2 * e.id) | 1].prev_out;
     1072      }
     1073      if(arcs[(2 * e.id) | 1].prev_out != -1) {
     1074        arcs[arcs[(2 * e.id) | 1].prev_out].next_out =
     1075          arcs[(2 * e.id) | 1].next_out;
     1076      } else {
     1077        nodes[arcs[2 * e.id].target].first_out =
     1078          arcs[(2 * e.id) | 1].next_out;
     1079      }
     1080
     1081      if (nodes[n.id].first_out != -1) {
     1082        arcs[nodes[n.id].first_out].prev_out = ((2 * e.id) | 1);
     1083      }
     1084      arcs[2 * e.id].target = n.id;
     1085      arcs[(2 * e.id) | 1].prev_out = -1;
     1086      arcs[(2 * e.id) | 1].next_out = nodes[n.id].first_out;
     1087      nodes[n.id].first_out = ((2 * e.id) | 1);
     1088    }
     1089
     1090  };
     1091
     1092//   typedef GraphExtender<UndirDigraphExtender<ListDigraphBase> >
     1093//   ExtendedListGraphBase;
     1094
     1095  typedef GraphExtender<ListGraphBase> ExtendedListGraphBase;
     1096
     1097
     1098
     1099  /// \addtogroup digraphs
     1100  /// @{
     1101
     1102  ///An undirected list digraph class.
     1103
     1104  ///This is a simple and fast undirected digraph implementation.
     1105  ///
     1106  ///An important extra feature of this digraph implementation is that
     1107  ///its maps are real \ref concepts::ReferenceMap "reference map"s.
     1108  ///
     1109  ///It conforms to the
     1110  ///\ref concepts::Graph "Graph concept".
     1111  ///
     1112  ///\sa concepts::Graph.
     1113  ///
     1114  class ListGraph : public ExtendedListGraphBase {
     1115  private:
     1116    ///ListGraph is \e not copy constructible. Use GraphCopy() instead.
     1117
     1118    ///ListGraph is \e not copy constructible. Use GraphCopy() instead.
     1119    ///
     1120    ListGraph(const ListGraph &) :ExtendedListGraphBase()  {};
     1121    ///\brief Assignment of ListGraph to another one is \e not allowed.
     1122    ///Use GraphCopy() instead.
     1123
     1124    ///Assignment of ListGraph to another one is \e not allowed.
     1125    ///Use GraphCopy() instead.
     1126    void operator=(const ListGraph &) {}
     1127  public:
     1128    /// Constructor
     1129   
     1130    /// Constructor.
     1131    ///
     1132    ListGraph() {}
     1133
     1134    typedef ExtendedListGraphBase Parent;
     1135
     1136    typedef Parent::OutArcIt IncArcIt;
     1137
     1138    /// \brief Add a new node to the digraph.
     1139    ///
     1140    /// \return the new node.
     1141    ///
     1142    Node addNode() { return Parent::addNode(); }
     1143
     1144    /// \brief Add a new edge to the digraph.
     1145    ///
     1146    /// Add a new arc to the digraph with source node \c s
     1147    /// and target node \c t.
     1148    /// \return the new edge.
     1149    Edge addEdge(const Node& s, const Node& t) {
     1150      return Parent::addEdge(s, t);
     1151    }
     1152    /// \brief Changes the source of \c e to \c n
     1153    ///
     1154    /// Changes the source of \c e to \c n
     1155    ///
     1156    ///\note The <tt>ArcIt</tt>s and <tt>InArcIt</tt>s
     1157    ///referencing the changed arc remain
     1158    ///valid. However <tt>OutArcIt</tt>s are invalidated.
     1159    void changeSource(Edge e, Node n) {
     1160      Parent::changeSource(e,n);
     1161    }   
     1162    /// \brief Changes the target of \c e to \c n
     1163    ///
     1164    /// Changes the target of \c e to \c n
     1165    ///
     1166    /// \note The <tt>ArcIt</tt>s referencing the changed arc remain
     1167    /// valid. However the other iterators may be invalidated.
     1168    void changeTarget(Edge e, Node n) {
     1169      Parent::changeTarget(e,n);
     1170    }
     1171    /// \brief Changes the source of \c e to \c n
     1172    ///
     1173    /// Changes the source of \c e to \c n. It changes the proper
     1174    /// node of the represented edge.
     1175    ///
     1176    ///\note The <tt>ArcIt</tt>s and <tt>InArcIt</tt>s
     1177    ///referencing the changed arc remain
     1178    ///valid. However <tt>OutArcIt</tt>s are invalidated.
     1179    void changeSource(Arc e, Node n) {
     1180      if (Parent::direction(e)) {
     1181        Parent::changeSource(e,n);
     1182      } else {
     1183        Parent::changeTarget(e,n);
     1184      }
     1185    }
     1186    /// \brief Changes the target of \c e to \c n
     1187    ///
     1188    /// Changes the target of \c e to \c n. It changes the proper
     1189    /// node of the represented edge.
     1190    ///
     1191    ///\note The <tt>ArcIt</tt>s and <tt>OutArcIt</tt>s
     1192    ///referencing the changed arc remain
     1193    ///valid. However <tt>InArcIt</tt>s are invalidated.
     1194    void changeTarget(Arc e, Node n) {
     1195      if (Parent::direction(e)) {
     1196        Parent::changeTarget(e,n);
     1197      } else {
     1198        Parent::changeSource(e,n);
     1199      }
     1200    }
     1201    /// \brief Contract two nodes.
     1202    ///
     1203    /// This function contracts two nodes.
     1204    ///
     1205    /// Node \p b will be removed but instead of deleting
     1206    /// its neighboring arcs, they will be joined to \p a.
     1207    /// The last parameter \p r controls whether to remove loops. \c true
     1208    /// means that loops will be removed.
     1209    ///
     1210    /// \note The <tt>ArcIt</tt>s referencing a moved arc remain
     1211    /// valid.
     1212    void contract(Node a, Node b, bool r = true) {
     1213      for(IncArcIt e(*this, b); e!=INVALID;) {
     1214        IncArcIt f = e; ++f;
     1215        if (r && runningNode(e) == a) {
     1216          erase(e);
     1217        } else if (source(e) == b) {
     1218          changeSource(e, a);
     1219        } else {
     1220          changeTarget(e, a);
     1221        }
     1222        e = f;
     1223      }
     1224      erase(b);
     1225    }
     1226
     1227
     1228    /// \brief Class to make a snapshot of the digraph and restore
     1229    /// to it later.
     1230    ///
     1231    /// Class to make a snapshot of the digraph and to restore it
     1232    /// later.
     1233    ///
     1234    /// The newly added nodes and edges can be removed
     1235    /// using the restore() function.
     1236    ///
     1237    /// \warning Arc and node deletions cannot be restored. This
     1238    /// events invalidate the snapshot.
     1239    class Snapshot {
     1240    protected:
     1241
     1242      typedef Parent::NodeNotifier NodeNotifier;
     1243
     1244      class NodeObserverProxy : public NodeNotifier::ObserverBase {
     1245      public:
     1246
     1247        NodeObserverProxy(Snapshot& _snapshot)
     1248          : snapshot(_snapshot) {}
     1249
     1250        using NodeNotifier::ObserverBase::attach;
     1251        using NodeNotifier::ObserverBase::detach;
     1252        using NodeNotifier::ObserverBase::attached;
     1253       
     1254      protected:
     1255       
     1256        virtual void add(const Node& node) {
     1257          snapshot.addNode(node);
     1258        }
     1259        virtual void add(const std::vector<Node>& nodes) {
     1260          for (int i = nodes.size() - 1; i >= 0; ++i) {
     1261            snapshot.addNode(nodes[i]);
     1262          }
     1263        }
     1264        virtual void erase(const Node& node) {
     1265          snapshot.eraseNode(node);
     1266        }
     1267        virtual void erase(const std::vector<Node>& nodes) {
     1268          for (int i = 0; i < int(nodes.size()); ++i) {
     1269            snapshot.eraseNode(nodes[i]);
     1270          }
     1271        }
     1272        virtual void build() {
     1273          Node node;
     1274          std::vector<Node> nodes;
     1275          for (notifier()->first(node); node != INVALID;
     1276               notifier()->next(node)) {
     1277            nodes.push_back(node);
     1278          }
     1279          for (int i = nodes.size() - 1; i >= 0; --i) {
     1280            snapshot.addNode(nodes[i]);
     1281          }
     1282        }
     1283        virtual void clear() {
     1284          Node node;
     1285          for (notifier()->first(node); node != INVALID;
     1286               notifier()->next(node)) {
     1287            snapshot.eraseNode(node);
     1288          }
     1289        }
     1290
     1291        Snapshot& snapshot;
     1292      };
     1293
     1294      class EdgeObserverProxy : public EdgeNotifier::ObserverBase {
     1295      public:
     1296
     1297        EdgeObserverProxy(Snapshot& _snapshot)
     1298          : snapshot(_snapshot) {}
     1299
     1300        using EdgeNotifier::ObserverBase::attach;
     1301        using EdgeNotifier::ObserverBase::detach;
     1302        using EdgeNotifier::ObserverBase::attached;
     1303       
     1304      protected:
     1305
     1306        virtual void add(const Edge& arc) {
     1307          snapshot.addEdge(arc);
     1308        }
     1309        virtual void add(const std::vector<Edge>& arcs) {
     1310          for (int i = arcs.size() - 1; i >= 0; ++i) {
     1311            snapshot.addEdge(arcs[i]);
     1312          }
     1313        }
     1314        virtual void erase(const Edge& arc) {
     1315          snapshot.eraseEdge(arc);
     1316        }
     1317        virtual void erase(const std::vector<Edge>& arcs) {
     1318          for (int i = 0; i < int(arcs.size()); ++i) {
     1319            snapshot.eraseEdge(arcs[i]);
     1320          }
     1321        }
     1322        virtual void build() {
     1323          Edge arc;
     1324          std::vector<Edge> arcs;
     1325          for (notifier()->first(arc); arc != INVALID;
     1326               notifier()->next(arc)) {
     1327            arcs.push_back(arc);
     1328          }
     1329          for (int i = arcs.size() - 1; i >= 0; --i) {
     1330            snapshot.addEdge(arcs[i]);
     1331          }
     1332        }
     1333        virtual void clear() {
     1334          Edge arc;
     1335          for (notifier()->first(arc); arc != INVALID;
     1336               notifier()->next(arc)) {
     1337            snapshot.eraseEdge(arc);
     1338          }
     1339        }
     1340
     1341        Snapshot& snapshot;
     1342      };
     1343     
     1344      ListGraph *digraph;
     1345
     1346      NodeObserverProxy node_observer_proxy;
     1347      EdgeObserverProxy arc_observer_proxy;
     1348
     1349      std::list<Node> added_nodes;
     1350      std::list<Edge> added_arcs;
     1351
     1352
     1353      void addNode(const Node& node) {
     1354        added_nodes.push_front(node);       
     1355      }
     1356      void eraseNode(const Node& node) {
     1357        std::list<Node>::iterator it =
     1358          std::find(added_nodes.begin(), added_nodes.end(), node);
     1359        if (it == added_nodes.end()) {
     1360          clear();
     1361          arc_observer_proxy.detach();
     1362          throw NodeNotifier::ImmediateDetach();
     1363        } else {
     1364          added_nodes.erase(it);
     1365        }
     1366      }
     1367
     1368      void addEdge(const Edge& arc) {
     1369        added_arcs.push_front(arc);       
     1370      }
     1371      void eraseEdge(const Edge& arc) {
     1372        std::list<Edge>::iterator it =
     1373          std::find(added_arcs.begin(), added_arcs.end(), arc);
     1374        if (it == added_arcs.end()) {
     1375          clear();
     1376          node_observer_proxy.detach();
     1377          throw EdgeNotifier::ImmediateDetach();
     1378        } else {
     1379          added_arcs.erase(it);
     1380        }       
     1381      }
     1382
     1383      void attach(ListGraph &_digraph) {
     1384        digraph = &_digraph;
     1385        node_observer_proxy.attach(digraph->notifier(Node()));
     1386        arc_observer_proxy.attach(digraph->notifier(Edge()));
     1387      }
     1388           
     1389      void detach() {
     1390        node_observer_proxy.detach();
     1391        arc_observer_proxy.detach();
     1392      }
     1393
     1394      bool attached() const {
     1395        return node_observer_proxy.attached();
     1396      }
     1397
     1398      void clear() {
     1399        added_nodes.clear();
     1400        added_arcs.clear();       
     1401      }
     1402
     1403    public:
     1404
     1405      /// \brief Default constructor.
     1406      ///
     1407      /// Default constructor.
     1408      /// To actually make a snapshot you must call save().
     1409      Snapshot()
     1410        : digraph(0), node_observer_proxy(*this),
     1411          arc_observer_proxy(*this) {}
     1412     
     1413      /// \brief Constructor that immediately makes a snapshot.
     1414      ///     
     1415      /// This constructor immediately makes a snapshot of the digraph.
     1416      /// \param _digraph The digraph we make a snapshot of.
     1417      Snapshot(ListGraph &_digraph)
     1418        : node_observer_proxy(*this),
     1419          arc_observer_proxy(*this) {
     1420        attach(_digraph);
     1421      }
     1422     
     1423      /// \brief Make a snapshot.
     1424      ///
     1425      /// Make a snapshot of the digraph.
     1426      ///
     1427      /// This function can be called more than once. In case of a repeated
     1428      /// call, the previous snapshot gets lost.
     1429      /// \param _digraph The digraph we make the snapshot of.
     1430      void save(ListGraph &_digraph) {
     1431        if (attached()) {
     1432          detach();
     1433          clear();
     1434        }
     1435        attach(_digraph);
     1436      }
     1437     
     1438      /// \brief Undo the changes until the last snapshot.
     1439      //
     1440      /// Undo the changes until the last snapshot created by save().
     1441      void restore() {
     1442        detach();
     1443        for(std::list<Edge>::iterator it = added_arcs.begin();
     1444            it != added_arcs.end(); ++it) {
     1445          digraph->erase(*it);
     1446        }
     1447        for(std::list<Node>::iterator it = added_nodes.begin();
     1448            it != added_nodes.end(); ++it) {
     1449          digraph->erase(*it);
     1450        }
     1451        clear();
     1452      }
     1453
     1454      /// \brief Gives back true when the snapshot is valid.
     1455      ///
     1456      /// Gives back true when the snapshot is valid.
     1457      bool valid() const {
     1458        return attached();
     1459      }
     1460    };
     1461  };
     1462 
     1463  /// @} 
     1464} //namespace lemon
     1465 
     1466
     1467#endif
Note: See TracChangeset for help on using the changeset viewer.