0
2
0
... | ... |
@@ -29,110 +29,111 @@ |
29 | 29 |
///\brief HypercubeGraph class. |
30 | 30 |
|
31 | 31 |
namespace lemon { |
32 | 32 |
|
33 | 33 |
class HypercubeGraphBase { |
34 | 34 |
|
35 | 35 |
public: |
36 | 36 |
|
37 | 37 |
typedef HypercubeGraphBase Graph; |
38 | 38 |
|
39 | 39 |
class Node; |
40 | 40 |
class Edge; |
41 | 41 |
class Arc; |
42 | 42 |
|
43 | 43 |
public: |
44 | 44 |
|
45 | 45 |
HypercubeGraphBase() {} |
46 | 46 |
|
47 | 47 |
protected: |
48 | 48 |
|
49 | 49 |
void construct(int dim) { |
50 | 50 |
LEMON_ASSERT(dim >= 1, "The number of dimensions must be at least 1."); |
51 | 51 |
_dim = dim; |
52 | 52 |
_node_num = 1 << dim; |
53 |
_edge_num = dim * (1 << dim-1); |
|
53 |
_edge_num = dim * (1 << (dim-1)); |
|
54 | 54 |
} |
55 | 55 |
|
56 | 56 |
public: |
57 | 57 |
|
58 | 58 |
typedef True NodeNumTag; |
59 | 59 |
typedef True EdgeNumTag; |
60 | 60 |
typedef True ArcNumTag; |
61 | 61 |
|
62 | 62 |
int nodeNum() const { return _node_num; } |
63 | 63 |
int edgeNum() const { return _edge_num; } |
64 | 64 |
int arcNum() const { return 2 * _edge_num; } |
65 | 65 |
|
66 | 66 |
int maxNodeId() const { return _node_num - 1; } |
67 | 67 |
int maxEdgeId() const { return _edge_num - 1; } |
68 | 68 |
int maxArcId() const { return 2 * _edge_num - 1; } |
69 | 69 |
|
70 | 70 |
static Node nodeFromId(int id) { return Node(id); } |
71 | 71 |
static Edge edgeFromId(int id) { return Edge(id); } |
72 | 72 |
static Arc arcFromId(int id) { return Arc(id); } |
73 | 73 |
|
74 | 74 |
static int id(Node node) { return node._id; } |
75 | 75 |
static int id(Edge edge) { return edge._id; } |
76 | 76 |
static int id(Arc arc) { return arc._id; } |
77 | 77 |
|
78 | 78 |
Node u(Edge edge) const { |
79 |
int base = edge._id & ((1 << _dim-1) - 1); |
|
80 |
int k = edge._id >> _dim-1; |
|
81 |
|
|
79 |
int base = edge._id & ((1 << (_dim-1)) - 1); |
|
80 |
int k = edge._id >> (_dim-1); |
|
81 |
return ((base >> k) << (k+1)) | (base & ((1 << k) - 1)); |
|
82 | 82 |
} |
83 | 83 |
|
84 | 84 |
Node v(Edge edge) const { |
85 |
int base = edge._id & ((1 << _dim-1) - 1); |
|
86 |
int k = edge._id >> _dim-1; |
|
87 |
|
|
85 |
int base = edge._id & ((1 << (_dim-1)) - 1); |
|
86 |
int k = edge._id >> (_dim-1); |
|
87 |
return ((base >> k) << (k+1)) | (base & ((1 << k) - 1)) | (1 << k); |
|
88 | 88 |
} |
89 | 89 |
|
90 | 90 |
Node source(Arc arc) const { |
91 | 91 |
return (arc._id & 1) == 1 ? u(arc) : v(arc); |
92 | 92 |
} |
93 | 93 |
|
94 | 94 |
Node target(Arc arc) const { |
95 | 95 |
return (arc._id & 1) == 1 ? v(arc) : u(arc); |
96 | 96 |
} |
97 | 97 |
|
98 | 98 |
typedef True FindEdgeTag; |
99 | 99 |
typedef True FindArcTag; |
100 | 100 |
|
101 | 101 |
Edge findEdge(Node u, Node v, Edge prev = INVALID) const { |
102 | 102 |
if (prev != INVALID) return INVALID; |
103 | 103 |
int d = u._id ^ v._id; |
104 | 104 |
int k = 0; |
105 | 105 |
if (d == 0) return INVALID; |
106 | 106 |
for ( ; (d & 1) == 0; d >>= 1) ++k; |
107 | 107 |
if (d >> 1 != 0) return INVALID; |
108 |
return (k << _dim-1) | ((u._id >> k+1) << k) | |
|
108 |
return (k << (_dim-1)) | ((u._id >> (k+1)) << k) | |
|
109 |
(u._id & ((1 << k) - 1)); |
|
109 | 110 |
} |
110 | 111 |
|
111 | 112 |
Arc findArc(Node u, Node v, Arc prev = INVALID) const { |
112 | 113 |
Edge edge = findEdge(u, v, prev); |
113 | 114 |
if (edge == INVALID) return INVALID; |
114 |
int k = edge._id >> _dim-1; |
|
115 |
int k = edge._id >> (_dim-1); |
|
115 | 116 |
return ((u._id >> k) & 1) == 1 ? edge._id << 1 : (edge._id << 1) | 1; |
116 | 117 |
} |
117 | 118 |
|
118 | 119 |
class Node { |
119 | 120 |
friend class HypercubeGraphBase; |
120 | 121 |
|
121 | 122 |
protected: |
122 | 123 |
int _id; |
123 | 124 |
Node(int id) : _id(id) {} |
124 | 125 |
public: |
125 | 126 |
Node() {} |
126 | 127 |
Node (Invalid) : _id(-1) {} |
127 | 128 |
bool operator==(const Node node) const {return _id == node._id;} |
128 | 129 |
bool operator!=(const Node node) const {return _id != node._id;} |
129 | 130 |
bool operator<(const Node node) const {return _id < node._id;} |
130 | 131 |
}; |
131 | 132 |
|
132 | 133 |
class Edge { |
133 | 134 |
friend class HypercubeGraphBase; |
134 | 135 |
friend class Arc; |
135 | 136 |
|
136 | 137 |
protected: |
137 | 138 |
int _id; |
138 | 139 |
|
... | ... |
@@ -173,109 +174,109 @@ |
173 | 174 |
|
174 | 175 |
void first(Edge& edge) const { |
175 | 176 |
edge._id = _edge_num - 1; |
176 | 177 |
} |
177 | 178 |
|
178 | 179 |
static void next(Edge& edge) { |
179 | 180 |
--edge._id; |
180 | 181 |
} |
181 | 182 |
|
182 | 183 |
void first(Arc& arc) const { |
183 | 184 |
arc._id = 2 * _edge_num - 1; |
184 | 185 |
} |
185 | 186 |
|
186 | 187 |
static void next(Arc& arc) { |
187 | 188 |
--arc._id; |
188 | 189 |
} |
189 | 190 |
|
190 | 191 |
void firstInc(Edge& edge, bool& dir, const Node& node) const { |
191 | 192 |
edge._id = node._id >> 1; |
192 | 193 |
dir = (node._id & 1) == 0; |
193 | 194 |
} |
194 | 195 |
|
195 | 196 |
void nextInc(Edge& edge, bool& dir) const { |
196 | 197 |
Node n = dir ? u(edge) : v(edge); |
197 |
int k = (edge._id >> _dim-1) + 1; |
|
198 |
int k = (edge._id >> (_dim-1)) + 1; |
|
198 | 199 |
if (k < _dim) { |
199 |
edge._id = (k << _dim-1) | |
|
200 |
((n._id >> k+1) << k) | (n._id & ((1 << k) - 1)); |
|
200 |
edge._id = (k << (_dim-1)) | |
|
201 |
((n._id >> (k+1)) << k) | (n._id & ((1 << k) - 1)); |
|
201 | 202 |
dir = ((n._id >> k) & 1) == 0; |
202 | 203 |
} else { |
203 | 204 |
edge._id = -1; |
204 | 205 |
dir = true; |
205 | 206 |
} |
206 | 207 |
} |
207 | 208 |
|
208 | 209 |
void firstOut(Arc& arc, const Node& node) const { |
209 | 210 |
arc._id = ((node._id >> 1) << 1) | (~node._id & 1); |
210 | 211 |
} |
211 | 212 |
|
212 | 213 |
void nextOut(Arc& arc) const { |
213 | 214 |
Node n = (arc._id & 1) == 1 ? u(arc) : v(arc); |
214 | 215 |
int k = (arc._id >> _dim) + 1; |
215 | 216 |
if (k < _dim) { |
216 |
arc._id = (k << _dim-1) | |
|
217 |
((n._id >> k+1) << k) | (n._id & ((1 << k) - 1)); |
|
217 |
arc._id = (k << (_dim-1)) | |
|
218 |
((n._id >> (k+1)) << k) | (n._id & ((1 << k) - 1)); |
|
218 | 219 |
arc._id = (arc._id << 1) | (~(n._id >> k) & 1); |
219 | 220 |
} else { |
220 | 221 |
arc._id = -1; |
221 | 222 |
} |
222 | 223 |
} |
223 | 224 |
|
224 | 225 |
void firstIn(Arc& arc, const Node& node) const { |
225 | 226 |
arc._id = ((node._id >> 1) << 1) | (node._id & 1); |
226 | 227 |
} |
227 | 228 |
|
228 | 229 |
void nextIn(Arc& arc) const { |
229 | 230 |
Node n = (arc._id & 1) == 1 ? v(arc) : u(arc); |
230 | 231 |
int k = (arc._id >> _dim) + 1; |
231 | 232 |
if (k < _dim) { |
232 |
arc._id = (k << _dim-1) | |
|
233 |
((n._id >> k+1) << k) | (n._id & ((1 << k) - 1)); |
|
233 |
arc._id = (k << (_dim-1)) | |
|
234 |
((n._id >> (k+1)) << k) | (n._id & ((1 << k) - 1)); |
|
234 | 235 |
arc._id = (arc._id << 1) | ((n._id >> k) & 1); |
235 | 236 |
} else { |
236 | 237 |
arc._id = -1; |
237 | 238 |
} |
238 | 239 |
} |
239 | 240 |
|
240 | 241 |
static bool direction(Arc arc) { |
241 | 242 |
return (arc._id & 1) == 1; |
242 | 243 |
} |
243 | 244 |
|
244 | 245 |
static Arc direct(Edge edge, bool dir) { |
245 | 246 |
return Arc((edge._id << 1) | (dir ? 1 : 0)); |
246 | 247 |
} |
247 | 248 |
|
248 | 249 |
int dimension() const { |
249 | 250 |
return _dim; |
250 | 251 |
} |
251 | 252 |
|
252 | 253 |
bool projection(Node node, int n) const { |
253 | 254 |
return static_cast<bool>(node._id & (1 << n)); |
254 | 255 |
} |
255 | 256 |
|
256 | 257 |
int dimension(Edge edge) const { |
257 |
return edge._id >> _dim-1; |
|
258 |
return edge._id >> (_dim-1); |
|
258 | 259 |
} |
259 | 260 |
|
260 | 261 |
int dimension(Arc arc) const { |
261 | 262 |
return arc._id >> _dim; |
262 | 263 |
} |
263 | 264 |
|
264 | 265 |
int index(Node node) const { |
265 | 266 |
return node._id; |
266 | 267 |
} |
267 | 268 |
|
268 | 269 |
Node operator()(int ix) const { |
269 | 270 |
return Node(ix); |
270 | 271 |
} |
271 | 272 |
|
272 | 273 |
private: |
273 | 274 |
int _dim; |
274 | 275 |
int _node_num, _edge_num; |
275 | 276 |
}; |
276 | 277 |
|
277 | 278 |
|
278 | 279 |
typedef GraphExtender<HypercubeGraphBase> ExtendedHypercubeGraphBase; |
279 | 280 |
|
280 | 281 |
/// \ingroup graphs |
281 | 282 |
/// |
... | ... |
@@ -300,80 +300,80 @@ |
300 | 300 |
checkGraphOutArcList(G, n, nb); |
301 | 301 |
checkGraphInArcList(G, n, nb); |
302 | 302 |
checkGraphIncEdgeList(G, n, nb); |
303 | 303 |
} |
304 | 304 |
|
305 | 305 |
checkArcDirections(G); |
306 | 306 |
|
307 | 307 |
checkGraphConArcList(G, 2 * (width * (height - 1) + (width - 1) * height)); |
308 | 308 |
checkGraphConEdgeList(G, width * (height - 1) + (width - 1) * height); |
309 | 309 |
|
310 | 310 |
checkNodeIds(G); |
311 | 311 |
checkArcIds(G); |
312 | 312 |
checkEdgeIds(G); |
313 | 313 |
checkGraphNodeMap(G); |
314 | 314 |
checkGraphArcMap(G); |
315 | 315 |
checkGraphEdgeMap(G); |
316 | 316 |
|
317 | 317 |
} |
318 | 318 |
|
319 | 319 |
void checkHypercubeGraph(int dim) { |
320 | 320 |
GRAPH_TYPEDEFS(HypercubeGraph); |
321 | 321 |
|
322 | 322 |
HypercubeGraph G(dim); |
323 | 323 |
checkGraphNodeList(G, 1 << dim); |
324 |
checkGraphEdgeList(G, dim * (1 << dim-1)); |
|
324 |
checkGraphEdgeList(G, dim * (1 << (dim-1))); |
|
325 | 325 |
checkGraphArcList(G, dim * (1 << dim)); |
326 | 326 |
|
327 | 327 |
Node n = G.nodeFromId(dim); |
328 | 328 |
|
329 | 329 |
for (NodeIt n(G); n != INVALID; ++n) { |
330 | 330 |
checkGraphIncEdgeList(G, n, dim); |
331 | 331 |
for (IncEdgeIt e(G, n); e != INVALID; ++e) { |
332 | 332 |
check( (G.u(e) == n && |
333 |
G.id(G.v(e)) == G.id(n) ^ (1 << G.dimension(e))) || |
|
333 |
G.id(G.v(e)) == (G.id(n) ^ (1 << G.dimension(e)))) || |
|
334 | 334 |
(G.v(e) == n && |
335 |
G.id(G.u(e)) == G.id(n) ^ (1 << G.dimension(e))), |
|
335 |
G.id(G.u(e)) == (G.id(n) ^ (1 << G.dimension(e)))), |
|
336 | 336 |
"Wrong edge or wrong dimension"); |
337 | 337 |
} |
338 | 338 |
|
339 | 339 |
checkGraphOutArcList(G, n, dim); |
340 | 340 |
for (OutArcIt a(G, n); a != INVALID; ++a) { |
341 | 341 |
check(G.source(a) == n && |
342 |
G.id(G.target(a)) == G.id(n) ^ (1 << G.dimension(a)), |
|
342 |
G.id(G.target(a)) == (G.id(n) ^ (1 << G.dimension(a))), |
|
343 | 343 |
"Wrong arc or wrong dimension"); |
344 | 344 |
} |
345 | 345 |
|
346 | 346 |
checkGraphInArcList(G, n, dim); |
347 | 347 |
for (InArcIt a(G, n); a != INVALID; ++a) { |
348 | 348 |
check(G.target(a) == n && |
349 |
G.id(G.source(a)) == G.id(n) ^ (1 << G.dimension(a)), |
|
349 |
G.id(G.source(a)) == (G.id(n) ^ (1 << G.dimension(a))), |
|
350 | 350 |
"Wrong arc or wrong dimension"); |
351 | 351 |
} |
352 | 352 |
} |
353 | 353 |
|
354 | 354 |
checkGraphConArcList(G, (1 << dim) * dim); |
355 |
checkGraphConEdgeList(G, dim * (1 << dim-1)); |
|
355 |
checkGraphConEdgeList(G, dim * (1 << (dim-1))); |
|
356 | 356 |
|
357 | 357 |
checkArcDirections(G); |
358 | 358 |
|
359 | 359 |
checkNodeIds(G); |
360 | 360 |
checkArcIds(G); |
361 | 361 |
checkEdgeIds(G); |
362 | 362 |
checkGraphNodeMap(G); |
363 | 363 |
checkGraphArcMap(G); |
364 | 364 |
checkGraphEdgeMap(G); |
365 | 365 |
} |
366 | 366 |
|
367 | 367 |
void checkGraphs() { |
368 | 368 |
{ // Checking ListGraph |
369 | 369 |
checkGraph<ListGraph>(); |
370 | 370 |
checkGraphValidityErase<ListGraph>(); |
371 | 371 |
} |
372 | 372 |
{ // Checking SmartGraph |
373 | 373 |
checkGraph<SmartGraph>(); |
374 | 374 |
checkGraphValidity<SmartGraph>(); |
375 | 375 |
} |
376 | 376 |
{ // Checking FullGraph |
377 | 377 |
checkFullGraph(7); |
378 | 378 |
checkFullGraph(8); |
379 | 379 |
} |
0 comments (0 inline)