LEMON/LEMON-main Changeset - r289:d91884dcd572

Repositories » LEMON » LEMON-main

Summary
Changelog
Files
Switch To
- loading...
Options
- Lightweight changelog
- Search
Pull Requests

Changeset - r289:d91884dcd572

« 288:47b3a3

293:47fbc8 »
292:e7af73 »

r289:d91884dcd572

Mon, 29 Sep 2008 12:34:08

[Not Reviewed]

0 comments (0 inline)

kpeter (Peter Kovacs)
kpeter@inf.elte.hu

Use DEBUG instead of ASSERT in graph extenders (ticket #17)

0 1 0

default

1 file changed with 8 insertions and 8 deletions:

lemon/bits/base_extender.h

↑ Collapse diff ↑

lemon/bits/base_extender.h

Show inline comments

/* -*- 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.
 *
 */

#ifndef LEMON_BITS_BASE_EXTENDER_H
#define LEMON_BITS_BASE_EXTENDER_H

#include <lemon/core.h>
#include <lemon/error.h>

#include <lemon/bits/map_extender.h>
#include <lemon/bits/default_map.h>

#include <lemon/concept_check.h>
#include <lemon/concepts/maps.h>

///\ingroup digraphbits
///\file
///\brief Extenders for the digraph types
namespace lemon {

  /// \ingroup digraphbits
  ///
  /// \brief BaseDigraph to BaseGraph extender
  template <typename Base>
  class UndirDigraphExtender : public Base {

  public:

    typedef Base Parent;
    typedef typename Parent::Arc Edge;
    typedef typename Parent::Node Node;

    typedef True UndirectedTag;

    class Arc : public Edge {
      friend class UndirDigraphExtender;

    protected:
      bool forward;

      Arc(const Edge &ue, bool _forward) :
        Edge(ue), forward(_forward) {}

    public:
      Arc() {}

      // Invalid arc constructor
      Arc(Invalid i) : Edge(i), forward(true) {}

      bool operator==(const Arc &that) const {
        return forward==that.forward && Edge(*this)==Edge(that);
      }
      bool operator!=(const Arc &that) const {
        return forward!=that.forward || Edge(*this)!=Edge(that);
      }
      bool operator<(const Arc &that) const {
        return forward<that.forward ||
          (!(that.forward<forward) && Edge(*this)<Edge(that));
      }
    };

    /// First node of the edge
    Node u(const Edge &e) const {
      return Parent::source(e);
    }

    /// Source of the given arc
    Node source(const Arc &e) const {
      return e.forward ? Parent::source(e) : Parent::target(e);
    }

    /// Second node of the edge
    Node v(const Edge &e) const {
      return Parent::target(e);
    }

    /// Target of the given arc
    Node target(const Arc &e) const {
      return e.forward ? Parent::target(e) : Parent::source(e);
    }

    /// \brief Directed arc from an edge.
    ///
    /// Returns a directed arc corresponding to the specified edge.
    /// If the given bool is true, the first node of the given edge and
    /// the source node of the returned arc are the same.
    static Arc direct(const Edge &e, bool d) {
      return Arc(e, d);
    }

    /// Returns whether the given directed arc has the same orientation
    /// as the corresponding edge.
    static bool direction(const Arc &a) { return a.forward; }

    using Parent::first;
    using Parent::next;

    void first(Arc &e) const {
      Parent::first(e);
      e.forward=true;
    }

    void next(Arc &e) const {
      if( e.forward ) {
        e.forward = false;
      }
      else {
        Parent::next(e);
        e.forward = true;
      }
    }

    void firstOut(Arc &e, const Node &n) const {
      Parent::firstIn(e,n);
      if( Edge(e) != INVALID ) {
        e.forward = false;
      }
      else {
        Parent::firstOut(e,n);
        e.forward = true;
      }
    }
    void nextOut(Arc &e) const {
      if( ! e.forward ) {
        Node n = Parent::target(e);
        Parent::nextIn(e);
        if( Edge(e) == INVALID ) {
          Parent::firstOut(e, n);
          e.forward = true;
        }
      }
      else {
        Parent::nextOut(e);
      }
    }

    void firstIn(Arc &e, const Node &n) const {
      Parent::firstOut(e,n);
      if( Edge(e) != INVALID ) {
        e.forward = false;
      }
      else {
        Parent::firstIn(e,n);
        e.forward = true;
      }
    }
    void nextIn(Arc &e) const {
      if( ! e.forward ) {
        Node n = Parent::source(e);
        Parent::nextOut(e);
        if( Edge(e) == INVALID ) {
          Parent::firstIn(e, n);
          e.forward = true;
        }
      }
      else {
        Parent::nextIn(e);
      }
    }

    void firstInc(Edge &e, bool &d, const Node &n) const {
      d = true;
      Parent::firstOut(e, n);
      if (e != INVALID) return;
      d = false;
      Parent::firstIn(e, n);
    }

    void nextInc(Edge &e, bool &d) const {
      if (d) {
        Node s = Parent::source(e);
        Parent::nextOut(e);
        if (e != INVALID) return;
        d = false;
        Parent::firstIn(e, s);
      } else {
        Parent::nextIn(e);
      }
    }

    Node nodeFromId(int ix) const {
      return Parent::nodeFromId(ix);
    }

    Arc arcFromId(int ix) const {
      return direct(Parent::arcFromId(ix >> 1), bool(ix & 1));
    }

    Edge edgeFromId(int ix) const {
      return Parent::arcFromId(ix);
    }

    int id(const Node &n) const {
      return Parent::id(n);
    }

    int id(const Edge &e) const {
      return Parent::id(e);
    }

    int id(const Arc &e) const {
      return 2 * Parent::id(e) + int(e.forward);
    }

    int maxNodeId() const {
      return Parent::maxNodeId();
    }

    int maxArcId() const {
      return 2 * Parent::maxArcId() + 1;
    }

    int maxEdgeId() const {
      return Parent::maxArcId();
    }

    int arcNum() const {
      return 2 * Parent::arcNum();
    }

    int edgeNum() const {
      return Parent::arcNum();
    }

    Arc findArc(Node s, Node t, Arc p = INVALID) const {
      if (p == INVALID) {
        Edge arc = Parent::findArc(s, t);
        if (arc != INVALID) return direct(arc, true);
        arc = Parent::findArc(t, s);
        if (arc != INVALID) return direct(arc, false);
      } else if (direction(p)) {
        Edge arc = Parent::findArc(s, t, p);
        if (arc != INVALID) return direct(arc, true);
        arc = Parent::findArc(t, s);
        if (arc != INVALID) return direct(arc, false);
      } else {
        Edge arc = Parent::findArc(t, s, p);
        if (arc != INVALID) return direct(arc, false);
      }
      return INVALID;
    }

    Edge findEdge(Node s, Node t, Edge p = INVALID) const {
      if (s != t) {
        if (p == INVALID) {
          Edge arc = Parent::findArc(s, t);
          if (arc != INVALID) return arc;
          arc = Parent::findArc(t, s);
          if (arc != INVALID) return arc;
        } else if (Parent::s(p) == s) {
          Edge arc = Parent::findArc(s, t, p);
          if (arc != INVALID) return arc;
          arc = Parent::findArc(t, s);
          if (arc != INVALID) return arc;
        } else {
          Edge arc = Parent::findArc(t, s, p);
          if (arc != INVALID) return arc;
        }
      } else {
        return Parent::findArc(s, t, p);
      }
      return INVALID;
    }
  };

  template <typename Base>
  class BidirBpGraphExtender : public Base {
  public:
    typedef Base Parent;
    typedef BidirBpGraphExtender Digraph;

    typedef typename Parent::Node Node;
    typedef typename Parent::Edge Edge;


    using Parent::first;
    using Parent::next;

    using Parent::id;

    class Red : public Node {
      friend class BidirBpGraphExtender;
    public:
      Red() {}
      Red(const Node& node) : Node(node) {
        LEMON_ASSERT(Parent::red(node) || node == INVALID,
                     typename Parent::NodeSetError());
        LEMON_DEBUG(Parent::red(node) || node == INVALID,
                    typename Parent::NodeSetError());
      }
      Red& operator=(const Node& node) {
        LEMON_ASSERT(Parent::red(node) || node == INVALID,
                     typename Parent::NodeSetError());
        LEMON_DEBUG(Parent::red(node) || node == INVALID,
                    typename Parent::NodeSetError());
        Node::operator=(node);
        return *this;
      }
      Red(Invalid) : Node(INVALID) {}
      Red& operator=(Invalid) {
        Node::operator=(INVALID);
        return *this;
      }
    };

    void first(Red& node) const {
      Parent::firstRed(static_cast<Node&>(node));
    }
    void next(Red& node) const {
      Parent::nextRed(static_cast<Node&>(node));
    }

    int id(const Red& node) const {
      return Parent::redId(node);
    }

    class Blue : public Node {
      friend class BidirBpGraphExtender;
    public:
      Blue() {}
      Blue(const Node& node) : Node(node) {
        LEMON_ASSERT(Parent::blue(node) || node == INVALID,
                     typename Parent::NodeSetError());
        LEMON_DEBUG(Parent::blue(node) || node == INVALID,
                    typename Parent::NodeSetError());
      }
      Blue& operator=(const Node& node) {
        LEMON_ASSERT(Parent::blue(node) || node == INVALID,
                     typename Parent::NodeSetError());
        LEMON_DEBUG(Parent::blue(node) || node == INVALID,
                    typename Parent::NodeSetError());
        Node::operator=(node);
        return *this;
      }
      Blue(Invalid) : Node(INVALID) {}
      Blue& operator=(Invalid) {
        Node::operator=(INVALID);
        return *this;
      }
    };

    void first(Blue& node) const {
      Parent::firstBlue(static_cast<Node&>(node));
    }
    void next(Blue& node) const {
      Parent::nextBlue(static_cast<Node&>(node));
    }

    int id(const Blue& node) const {
      return Parent::redId(node);
    }

    Node source(const Edge& arc) const {
      return red(arc);
    }
    Node target(const Edge& arc) const {
      return blue(arc);
    }

    void firstInc(Edge& arc, bool& dir, const Node& node) const {
      if (Parent::red(node)) {
        Parent::firstFromRed(arc, node);
        dir = true;
      } else {
        Parent::firstFromBlue(arc, node);
        dir = static_cast<Edge&>(arc) == INVALID;
      }
    }
    void nextInc(Edge& arc, bool& dir) const {
      if (dir) {
        Parent::nextFromRed(arc);
      } else {
        Parent::nextFromBlue(arc);
        if (arc == INVALID) dir = true;
      }
    }

    class Arc : public Edge {
      friend class BidirBpGraphExtender;
    protected:
      bool forward;

      Arc(const Edge& arc, bool _forward)
        : Edge(arc), forward(_forward) {}

    public:
      Arc() {}
      Arc (Invalid) : Edge(INVALID), forward(true) {}
      bool operator==(const Arc& i) const {
        return Edge::operator==(i) && forward == i.forward;
      }
      bool operator!=(const Arc& i) const {
        return Edge::operator!=(i) || forward != i.forward;
      }
      bool operator<(const Arc& i) const {
        return Edge::operator<(i) ||
          (!(i.forward<forward) && Edge(*this)<Edge(i));
      }
    };

    void first(Arc& arc) const {
      Parent::first(static_cast<Edge&>(arc));
      arc.forward = true;
    }

    void next(Arc& arc) const {
      if (!arc.forward) {
        Parent::next(static_cast<Edge&>(arc));
      }
      arc.forward = !arc.forward;
    }

    void firstOut(Arc& arc, const Node& node) const {
      if (Parent::red(node)) {
        Parent::firstFromRed(arc, node);
        arc.forward = true;
      } else {
        Parent::firstFromBlue(arc, node);
        arc.forward = static_cast<Edge&>(arc) == INVALID;
      }
    }
    void nextOut(Arc& arc) const {
      if (arc.forward) {
        Parent::nextFromRed(arc);
      } else {
        Parent::nextFromBlue(arc);
        arc.forward = static_cast<Edge&>(arc) == INVALID;
      }
    }

    void firstIn(Arc& arc, const Node& node) const {
      if (Parent::blue(node)) {
        Parent::firstFromBlue(arc, node);
        arc.forward = true;
      } else {
        Parent::firstFromRed(arc, node);
        arc.forward = static_cast<Edge&>(arc) == INVALID;
      }
    }
    void nextIn(Arc& arc) const {
      if (arc.forward) {
        Parent::nextFromBlue(arc);
      } else {
        Parent::nextFromRed(arc);
        arc.forward = static_cast<Edge&>(arc) == INVALID;
      }
    }

    Node source(const Arc& arc) const {
      return arc.forward ? Parent::red(arc) : Parent::blue(arc);
    }
    Node target(const Arc& arc) const {
      return arc.forward ? Parent::blue(arc) : Parent::red(arc);
    }

    int id(const Arc& arc) const {
      return (Parent::id(static_cast<const Edge&>(arc)) << 1) +
        (arc.forward ? 0 : 1);
    }
    Arc arcFromId(int ix) const {
      return Arc(Parent::fromEdgeId(ix >> 1), (ix & 1) == 0);
    }
    int maxArcId() const {
      return (Parent::maxEdgeId() << 1) + 1;
    }

    bool direction(const Arc& arc) const {
      return arc.forward;
    }

    Arc direct(const Edge& arc, bool dir) const {
      return Arc(arc, dir);
    }

    int arcNum() const {
      return 2 * Parent::edgeNum();
    }

    int edgeNum() const {
      return Parent::edgeNum();
    }


  };
}

#endif

0 comments (0 inline)

RhodeCode

Login to your account

1	1	/* -- mode: C++; indent-tabs-mode: nil; --
2	2	*
3	3	* This file is a part of LEMON, a generic C++ optimization library.
4	4	*
5	5	* Copyright (C) 2003-2008
6	6	* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
7	7	* (Egervary Research Group on Combinatorial Optimization, EGRES).
8	8	*
9	9	* Permission to use, modify and distribute this software is granted
10	10	* provided that this copyright notice appears in all copies. For
11	11	* precise terms see the accompanying LICENSE file.
12	12	*
13	13	* This software is provided "AS IS" with no warranty of any kind,
14	14	* express or implied, and with no claim as to its suitability for any
15	15	* purpose.
16	16	*
17	17	*/
18	18
19	19	#ifndef LEMON_BITS_BASE_EXTENDER_H
20	20	#define LEMON_BITS_BASE_EXTENDER_H
21	21
22	22	#include <lemon/core.h>
23	23	#include <lemon/error.h>
24	24
25	25	#include <lemon/bits/map_extender.h>
26	26	#include <lemon/bits/default_map.h>
27	27
28	28	#include <lemon/concept_check.h>
29	29	#include <lemon/concepts/maps.h>
30	30
31	31	///\ingroup digraphbits
32	32	///\file
33	33	///\brief Extenders for the digraph types
34	34	namespace lemon {
35	35
36	36	/// \ingroup digraphbits
37	37	///
38	38	/// \brief BaseDigraph to BaseGraph extender
39	39	template <typename Base>
40	40	class UndirDigraphExtender : public Base {
41	41
42	42	public:
43	43
44	44	typedef Base Parent;
45	45	typedef typename Parent::Arc Edge;
46	46	typedef typename Parent::Node Node;
47	47
48	48	typedef True UndirectedTag;
49	49
50	50	class Arc : public Edge {
51	51	friend class UndirDigraphExtender;
52	52
53	53	protected:
54	54	bool forward;
55	55
56	56	Arc(const Edge &ue, bool _forward) :
57	57	Edge(ue), forward(_forward) {}
58	58
59	59	public:
60	60	Arc() {}
61	61
62	62	// Invalid arc constructor
63	63	Arc(Invalid i) : Edge(i), forward(true) {}
64	64
65	65	bool operator==(const Arc &that) const {
66	66	return forward==that.forward && Edge(*this)==Edge(that);
67	67	}
68	68	bool operator!=(const Arc &that) const {
69	69	return forward!=that.forward \|\| Edge(*this)!=Edge(that);
70	70	}
71	71	bool operator<(const Arc &that) const {
72	72	return forward<that.forward \|\|
73	73	(!(that.forward<forward) && Edge(*this)<Edge(that));
74	74	}
75	75	};
76	76
77	77	/// First node of the edge
78	78	Node u(const Edge &e) const {
79	79	return Parent::source(e);
80	80	}
81	81
82	82	/// Source of the given arc
83	83	Node source(const Arc &e) const {
84	84	return e.forward ? Parent::source(e) : Parent::target(e);
85	85	}
86	86
87	87	/// Second node of the edge
88	88	Node v(const Edge &e) const {
89	89	return Parent::target(e);
90	90	}
91	91
92	92	/// Target of the given arc
93	93	Node target(const Arc &e) const {
94	94	return e.forward ? Parent::target(e) : Parent::source(e);
95	95	}
96	96
97	97	/// \brief Directed arc from an edge.
98	98	///
99	99	/// Returns a directed arc corresponding to the specified edge.
100	100	/// If the given bool is true, the first node of the given edge and
101	101	/// the source node of the returned arc are the same.
102	102	static Arc direct(const Edge &e, bool d) {
103	103	return Arc(e, d);
104	104	}
105	105
106	106	/// Returns whether the given directed arc has the same orientation
107	107	/// as the corresponding edge.
108	108	static bool direction(const Arc &a) { return a.forward; }
109	109
110	110	using Parent::first;
111	111	using Parent::next;
112	112
113	113	void first(Arc &e) const {
114	114	Parent::first(e);
115	115	e.forward=true;
116	116	}
117	117
118	118	void next(Arc &e) const {
119	119	if( e.forward ) {
120	120	e.forward = false;
121	121	}
122	122	else {
123	123	Parent::next(e);
124	124	e.forward = true;
125	125	}
126	126	}
127	127
128	128	void firstOut(Arc &e, const Node &n) const {
129	129	Parent::firstIn(e,n);
130	130	if( Edge(e) != INVALID ) {
131	131	e.forward = false;
132	132	}
133	133	else {
134	134	Parent::firstOut(e,n);
135	135	e.forward = true;
136	136	}
137	137	}
138	138	void nextOut(Arc &e) const {
139	139	if( ! e.forward ) {
140	140	Node n = Parent::target(e);
141	141	Parent::nextIn(e);
142	142	if( Edge(e) == INVALID ) {
143	143	Parent::firstOut(e, n);
144	144	e.forward = true;
145	145	}
146	146	}
147	147	else {
148	148	Parent::nextOut(e);
149	149	}
150	150	}
151	151
152	152	void firstIn(Arc &e, const Node &n) const {
153	153	Parent::firstOut(e,n);
154	154	if( Edge(e) != INVALID ) {
155	155	e.forward = false;
156	156	}
157	157	else {
158	158	Parent::firstIn(e,n);
159	159	e.forward = true;
160	160	}
161	161	}
162	162	void nextIn(Arc &e) const {
163	163	if( ! e.forward ) {
164	164	Node n = Parent::source(e);
165	165	Parent::nextOut(e);
166	166	if( Edge(e) == INVALID ) {
167	167	Parent::firstIn(e, n);
168	168	e.forward = true;
169	169	}
170	170	}
171	171	else {
172	172	Parent::nextIn(e);
173	173	}
174	174	}
175	175
176	176	void firstInc(Edge &e, bool &d, const Node &n) const {
177	177	d = true;
178	178	Parent::firstOut(e, n);
179	179	if (e != INVALID) return;
180	180	d = false;
181	181	Parent::firstIn(e, n);
182	182	}
183	183
184	184	void nextInc(Edge &e, bool &d) const {
185	185	if (d) {
186	186	Node s = Parent::source(e);
187	187	Parent::nextOut(e);
188	188	if (e != INVALID) return;
189	189	d = false;
190	190	Parent::firstIn(e, s);
191	191	} else {
192	192	Parent::nextIn(e);
193	193	}
194	194	}
195	195
196	196	Node nodeFromId(int ix) const {
197	197	return Parent::nodeFromId(ix);
198	198	}
199	199
200	200	Arc arcFromId(int ix) const {
201	201	return direct(Parent::arcFromId(ix >> 1), bool(ix & 1));
202	202	}
203	203
204	204	Edge edgeFromId(int ix) const {
205	205	return Parent::arcFromId(ix);
206	206	}
207	207
208	208	int id(const Node &n) const {
209	209	return Parent::id(n);
210	210	}
211	211
212	212	int id(const Edge &e) const {
213	213	return Parent::id(e);
214	214	}
215	215
216	216	int id(const Arc &e) const {
217	217	return 2 * Parent::id(e) + int(e.forward);
218	218	}
219	219
220	220	int maxNodeId() const {
221	221	return Parent::maxNodeId();
222	222	}
223	223
224	224	int maxArcId() const {
225	225	return 2 * Parent::maxArcId() + 1;
226	226	}
227	227
228	228	int maxEdgeId() const {
229	229	return Parent::maxArcId();
230	230	}
231	231
232	232	int arcNum() const {
233	233	return 2 * Parent::arcNum();
234	234	}
235	235
236	236	int edgeNum() const {
237	237	return Parent::arcNum();
238	238	}
239	239
240	240	Arc findArc(Node s, Node t, Arc p = INVALID) const {
241	241	if (p == INVALID) {
242	242	Edge arc = Parent::findArc(s, t);
243	243	if (arc != INVALID) return direct(arc, true);
244	244	arc = Parent::findArc(t, s);
245	245	if (arc != INVALID) return direct(arc, false);
246	246	} else if (direction(p)) {
247	247	Edge arc = Parent::findArc(s, t, p);
248	248	if (arc != INVALID) return direct(arc, true);
249	249	arc = Parent::findArc(t, s);
250	250	if (arc != INVALID) return direct(arc, false);
251	251	} else {
252	252	Edge arc = Parent::findArc(t, s, p);
253	253	if (arc != INVALID) return direct(arc, false);
254	254	}
255	255	return INVALID;
256	256	}
257	257
258	258	Edge findEdge(Node s, Node t, Edge p = INVALID) const {
259	259	if (s != t) {
260	260	if (p == INVALID) {
261	261	Edge arc = Parent::findArc(s, t);
262	262	if (arc != INVALID) return arc;
263	263	arc = Parent::findArc(t, s);
264	264	if (arc != INVALID) return arc;
265	265	} else if (Parent::s(p) == s) {
266	266	Edge arc = Parent::findArc(s, t, p);
267	267	if (arc != INVALID) return arc;
268	268	arc = Parent::findArc(t, s);
269	269	if (arc != INVALID) return arc;
270	270	} else {
271	271	Edge arc = Parent::findArc(t, s, p);
272	272	if (arc != INVALID) return arc;
273	273	}
274	274	} else {
275	275	return Parent::findArc(s, t, p);
276	276	}
277	277	return INVALID;
278	278	}
279	279	};
280	280
281	281	template <typename Base>
282	282	class BidirBpGraphExtender : public Base {
283	283	public:
284	284	typedef Base Parent;
285	285	typedef BidirBpGraphExtender Digraph;
286	286
287	287	typedef typename Parent::Node Node;
288	288	typedef typename Parent::Edge Edge;
289	289
290	290
291	291	using Parent::first;
292	292	using Parent::next;
293	293
294	294	using Parent::id;
295	295
296	296	class Red : public Node {
297	297	friend class BidirBpGraphExtender;
298	298	public:
299	299	Red() {}
300	300	Red(const Node& node) : Node(node) {
301		LEMON_ASSERT(Parent::red(node) \|\| node == INVALID,
302		typename Parent::NodeSetError());
	301	LEMON_DEBUG(Parent::red(node) \|\| node == INVALID,
	302	typename Parent::NodeSetError());
303	303	}
304	304	Red& operator=(const Node& node) {
305		LEMON_ASSERT(Parent::red(node) \|\| node == INVALID,
306		typename Parent::NodeSetError());
	305	LEMON_DEBUG(Parent::red(node) \|\| node == INVALID,
	306	typename Parent::NodeSetError());
307	307	Node::operator=(node);
308	308	return *this;
309	309	}
310	310	Red(Invalid) : Node(INVALID) {}
311	311	Red& operator=(Invalid) {
312	312	Node::operator=(INVALID);
313	313	return *this;
314	314	}
315	315	};
316	316
317	317	void first(Red& node) const {
318	318	Parent::firstRed(static_cast<Node&>(node));
319	319	}
320	320	void next(Red& node) const {
321	321	Parent::nextRed(static_cast<Node&>(node));
322	322	}
323	323
324	324	int id(const Red& node) const {
325	325	return Parent::redId(node);
326	326	}
327	327
328	328	class Blue : public Node {
329	329	friend class BidirBpGraphExtender;
330	330	public:
331	331	Blue() {}
332	332	Blue(const Node& node) : Node(node) {
333		LEMON_ASSERT(Parent::blue(node) \|\| node == INVALID,
334		typename Parent::NodeSetError());
	333	LEMON_DEBUG(Parent::blue(node) \|\| node == INVALID,
	334	typename Parent::NodeSetError());
335	335	}
336	336	Blue& operator=(const Node& node) {
337		LEMON_ASSERT(Parent::blue(node) \|\| node == INVALID,
338		typename Parent::NodeSetError());
	337	LEMON_DEBUG(Parent::blue(node) \|\| node == INVALID,
	338	typename Parent::NodeSetError());
339	339	Node::operator=(node);
340	340	return *this;
341	341	}
342	342	Blue(Invalid) : Node(INVALID) {}
343	343	Blue& operator=(Invalid) {
344	344	Node::operator=(INVALID);
345	345	return *this;
346	346	}
347	347	};
348	348
349	349	void first(Blue& node) const {
350	350	Parent::firstBlue(static_cast<Node&>(node));
351	351	}
352	352	void next(Blue& node) const {
353	353	Parent::nextBlue(static_cast<Node&>(node));
354	354	}
355	355
356	356	int id(const Blue& node) const {
357	357	return Parent::redId(node);
358	358	}
359	359
360	360	Node source(const Edge& arc) const {
361	361	return red(arc);
362	362	}
363	363	Node target(const Edge& arc) const {
364	364	return blue(arc);
365	365	}
366	366
367	367	void firstInc(Edge& arc, bool& dir, const Node& node) const {
368	368	if (Parent::red(node)) {
369	369	Parent::firstFromRed(arc, node);
370	370	dir = true;
371	371	} else {
372	372	Parent::firstFromBlue(arc, node);
373	373	dir = static_cast<Edge&>(arc) == INVALID;
374	374	}
375	375	}
376	376	void nextInc(Edge& arc, bool& dir) const {
377	377	if (dir) {
378	378	Parent::nextFromRed(arc);
379	379	} else {
380	380	Parent::nextFromBlue(arc);
381	381	if (arc == INVALID) dir = true;
382	382	}
383	383	}
384	384
385	385	class Arc : public Edge {
386	386	friend class BidirBpGraphExtender;
387	387	protected:
388	388	bool forward;
389	389
390	390	Arc(const Edge& arc, bool _forward)
391	391	: Edge(arc), forward(_forward) {}
392	392
393	393	public:
394	394	Arc() {}
395	395	Arc (Invalid) : Edge(INVALID), forward(true) {}
396	396	bool operator==(const Arc& i) const {
397	397	return Edge::operator==(i) && forward == i.forward;
398	398	}
399	399	bool operator!=(const Arc& i) const {
400	400	return Edge::operator!=(i) \|\| forward != i.forward;
401	401	}
402	402	bool operator<(const Arc& i) const {
403	403	return Edge::operator<(i) \|\|
404	404	(!(i.forward<forward) && Edge(*this)<Edge(i));
405	405	}
406	406	};
407	407
408	408	void first(Arc& arc) const {
409	409	Parent::first(static_cast<Edge&>(arc));
410	410	arc.forward = true;
411	411	}
412	412
413	413	void next(Arc& arc) const {
414	414	if (!arc.forward) {
415	415	Parent::next(static_cast<Edge&>(arc));
416	416	}
417	417	arc.forward = !arc.forward;
418	418	}
419	419
420	420	void firstOut(Arc& arc, const Node& node) const {
421	421	if (Parent::red(node)) {
422	422	Parent::firstFromRed(arc, node);
423	423	arc.forward = true;
424	424	} else {
425	425	Parent::firstFromBlue(arc, node);
426	426	arc.forward = static_cast<Edge&>(arc) == INVALID;
427	427	}
428	428	}
429	429	void nextOut(Arc& arc) const {
430	430	if (arc.forward) {
431	431	Parent::nextFromRed(arc);
432	432	} else {
433	433	Parent::nextFromBlue(arc);
434	434	arc.forward = static_cast<Edge&>(arc) == INVALID;
435	435	}
436	436	}
437	437
438	438	void firstIn(Arc& arc, const Node& node) const {
439	439	if (Parent::blue(node)) {
440	440	Parent::firstFromBlue(arc, node);
441	441	arc.forward = true;
442	442	} else {
443	443	Parent::firstFromRed(arc, node);
444	444	arc.forward = static_cast<Edge&>(arc) == INVALID;
445	445	}
446	446	}
447	447	void nextIn(Arc& arc) const {
448	448	if (arc.forward) {
449	449	Parent::nextFromBlue(arc);
450	450	} else {
451	451	Parent::nextFromRed(arc);
452	452	arc.forward = static_cast<Edge&>(arc) == INVALID;
453	453	}
454	454	}
455	455
456	456	Node source(const Arc& arc) const {
457	457	return arc.forward ? Parent::red(arc) : Parent::blue(arc);
458	458	}
459	459	Node target(const Arc& arc) const {
460	460	return arc.forward ? Parent::blue(arc) : Parent::red(arc);
461	461	}
462	462
463	463	int id(const Arc& arc) const {
464	464	return (Parent::id(static_cast<const Edge&>(arc)) << 1) +
465	465	(arc.forward ? 0 : 1);
466	466	}
467	467	Arc arcFromId(int ix) const {
468	468	return Arc(Parent::fromEdgeId(ix >> 1), (ix & 1) == 0);
469	469	}
470	470	int maxArcId() const {
471	471	return (Parent::maxEdgeId() << 1) + 1;
472	472	}
473	473
474	474	bool direction(const Arc& arc) const {
475	475	return arc.forward;
476	476	}
477	477
478	478	Arc direct(const Edge& arc, bool dir) const {
479	479	return Arc(arc, dir);
480	480	}
481	481
482	482	int arcNum() const {
483	483	return 2 * Parent::edgeNum();
484	484	}
485	485
486	486	int edgeNum() const {
487	487	return Parent::edgeNum();
488	488	}
489	489
490	490
491	491	};
492	492	}
493	493
494	494	#endif