0
2
0
... | ... |
@@ -184,32 +184,73 @@ |
184 | 184 |
int target = nodeRef[digraph.target(*it)].id; |
185 | 185 |
arcRef[*it] = Arc(arc_index); |
186 | 186 |
arc_source[arc_index] = source; |
187 | 187 |
arc_target[arc_index] = target; |
188 | 188 |
arc_next_in[arc_index] = node_first_in[target]; |
189 | 189 |
node_first_in[target] = arc_index; |
190 | 190 |
arc_next_out[arc_index] = arc_index + 1; |
191 | 191 |
++arc_index; |
192 | 192 |
} |
193 | 193 |
arc_next_out[arc_index - 1] = -1; |
194 | 194 |
} else { |
195 | 195 |
node_first_out[source] = arc_index; |
196 | 196 |
} |
197 | 197 |
} |
198 | 198 |
node_first_out[node_num] = arc_num; |
199 | 199 |
} |
200 |
|
|
201 |
template <typename ArcListIterator> |
|
202 |
void build(int n, ArcListIterator first, ArcListIterator last) { |
|
203 |
built = true; |
|
204 |
|
|
205 |
node_num = n; |
|
206 |
arc_num = std::distance(first, last); |
|
207 |
|
|
208 |
node_first_out = new int[node_num + 1]; |
|
209 |
node_first_in = new int[node_num]; |
|
210 |
|
|
211 |
arc_source = new int[arc_num]; |
|
212 |
arc_target = new int[arc_num]; |
|
213 |
arc_next_out = new int[arc_num]; |
|
214 |
arc_next_in = new int[arc_num]; |
|
215 |
|
|
216 |
for (int i = 0; i != node_num; ++i) { |
|
217 |
node_first_in[i] = -1; |
|
218 |
} |
|
219 |
|
|
220 |
int arc_index = 0; |
|
221 |
for (int i = 0; i != node_num; ++i) { |
|
222 |
node_first_out[i] = arc_index; |
|
223 |
for ( ; first != last && (*first).first == i; ++first) { |
|
224 |
int j = (*first).second; |
|
225 |
LEMON_ASSERT(j >= 0 && j < node_num, |
|
226 |
"Wrong arc list for StaticDigraph::build()"); |
|
227 |
arc_source[arc_index] = i; |
|
228 |
arc_target[arc_index] = j; |
|
229 |
arc_next_in[arc_index] = node_first_in[j]; |
|
230 |
node_first_in[j] = arc_index; |
|
231 |
arc_next_out[arc_index] = arc_index + 1; |
|
232 |
++arc_index; |
|
233 |
} |
|
234 |
if (arc_index > node_first_out[i]) |
|
235 |
arc_next_out[arc_index - 1] = -1; |
|
236 |
} |
|
237 |
LEMON_ASSERT(first == last, |
|
238 |
"Wrong arc list for StaticDigraph::build()"); |
|
239 |
node_first_out[node_num] = arc_num; |
|
240 |
} |
|
200 | 241 |
|
201 | 242 |
protected: |
202 | 243 |
|
203 | 244 |
void fastFirstOut(Arc& e, const Node& n) const { |
204 | 245 |
e.id = node_first_out[n.id]; |
205 | 246 |
} |
206 | 247 |
|
207 | 248 |
static void fastNextOut(Arc& e) { |
208 | 249 |
++e.id; |
209 | 250 |
} |
210 | 251 |
void fastLastOut(Arc& e, const Node& n) const { |
211 | 252 |
e.id = node_first_out[n.id + 1]; |
212 | 253 |
} |
213 | 254 |
|
214 | 255 |
protected: |
215 | 256 |
bool built; |
... | ... |
@@ -315,32 +356,70 @@ |
315 | 356 |
/// concept. |
316 | 357 |
/// \param arcRef The arc references will be copied into this map. |
317 | 358 |
/// Its key type must be \c Digraph::Arc and its value type must be |
318 | 359 |
/// \c StaticDigraph::Arc. |
319 | 360 |
/// It must conform to the \ref concepts::WriteMap "WriteMap" concept. |
320 | 361 |
/// |
321 | 362 |
/// \note If you do not need the arc references, then you could use |
322 | 363 |
/// \ref NullMap for the last parameter. However the node references |
323 | 364 |
/// are required by the function itself, thus they must be readable |
324 | 365 |
/// from the map. |
325 | 366 |
template <typename Digraph, typename NodeRefMap, typename ArcRefMap> |
326 | 367 |
void build(const Digraph& digraph, NodeRefMap& nodeRef, ArcRefMap& arcRef) { |
327 | 368 |
if (built) Parent::clear(); |
328 | 369 |
Parent::build(digraph, nodeRef, arcRef); |
329 | 370 |
} |
330 | 371 |
|
372 |
/// \brief Build the digraph from an arc list. |
|
373 |
/// |
|
374 |
/// This function builds the digraph from the given arc list. |
|
375 |
/// It can be called more than once, but in such case, the whole |
|
376 |
/// structure and all maps will be cleared and rebuilt. |
|
377 |
/// |
|
378 |
/// The list of the arcs must be given in the range <tt>[begin, end)</tt> |
|
379 |
/// specified by STL compatible itartors whose \c value_type must be |
|
380 |
/// <tt>std::pair<int,int></tt>. |
|
381 |
/// Each arc must be specified by a pair of integer indices |
|
382 |
/// from the range <tt>[0..n-1]</tt>. <i>The pairs must be in a |
|
383 |
/// non-decreasing order with respect to their first values.</i> |
|
384 |
/// If the k-th pair in the list is <tt>(i,j)</tt>, then |
|
385 |
/// <tt>arc(k-1)</tt> will connect <tt>node(i)</tt> to <tt>node(j)</tt>. |
|
386 |
/// |
|
387 |
/// \param n The number of nodes. |
|
388 |
/// \param begin An iterator pointing to the beginning of the arc list. |
|
389 |
/// \param end An iterator pointing to the end of the arc list. |
|
390 |
/// |
|
391 |
/// For example, a simple digraph can be constructed like this. |
|
392 |
/// \code |
|
393 |
/// std::vector<std::pair<int,int> > arcs; |
|
394 |
/// arcs.push_back(std::make_pair(0,1)); |
|
395 |
/// arcs.push_back(std::make_pair(0,2)); |
|
396 |
/// arcs.push_back(std::make_pair(1,3)); |
|
397 |
/// arcs.push_back(std::make_pair(1,2)); |
|
398 |
/// arcs.push_back(std::make_pair(3,0)); |
|
399 |
/// StaticDigraph gr; |
|
400 |
/// gr.build(4, arcs.begin(), arcs.end()); |
|
401 |
/// \endcode |
|
402 |
template <typename ArcListIterator> |
|
403 |
void build(int n, ArcListIterator begin, ArcListIterator end) { |
|
404 |
if (built) Parent::clear(); |
|
405 |
StaticDigraphBase::build(n, begin, end); |
|
406 |
notifier(Node()).build(); |
|
407 |
notifier(Arc()).build(); |
|
408 |
} |
|
409 |
|
|
331 | 410 |
/// \brief Clear the digraph. |
332 | 411 |
/// |
333 | 412 |
/// This function erases all nodes and arcs from the digraph. |
334 | 413 |
void clear() { |
335 | 414 |
Parent::clear(); |
336 | 415 |
} |
337 | 416 |
|
338 | 417 |
protected: |
339 | 418 |
|
340 | 419 |
using Parent::fastFirstOut; |
341 | 420 |
using Parent::fastNextOut; |
342 | 421 |
using Parent::fastLastOut; |
343 | 422 |
|
344 | 423 |
public: |
345 | 424 |
|
346 | 425 |
class OutArcIt : public Arc { |
... | ... |
@@ -428,32 +428,64 @@ |
428 | 428 |
|
429 | 429 |
digraphCopy(g, G).nodeRef(nref).run(); |
430 | 430 |
|
431 | 431 |
checkGraphNodeList(G, 3); |
432 | 432 |
checkGraphArcList(G, 4); |
433 | 433 |
|
434 | 434 |
checkGraphOutArcList(G, nref[n1], 1); |
435 | 435 |
checkGraphOutArcList(G, nref[n2], 3); |
436 | 436 |
checkGraphOutArcList(G, nref[n3], 0); |
437 | 437 |
|
438 | 438 |
checkGraphInArcList(G, nref[n1], 1); |
439 | 439 |
checkGraphInArcList(G, nref[n2], 1); |
440 | 440 |
checkGraphInArcList(G, nref[n3], 2); |
441 | 441 |
|
442 | 442 |
checkGraphConArcList(G, 4); |
443 | 443 |
|
444 |
std::vector<std::pair<int,int> > arcs; |
|
445 |
arcs.push_back(std::make_pair(0,1)); |
|
446 |
arcs.push_back(std::make_pair(0,2)); |
|
447 |
arcs.push_back(std::make_pair(1,3)); |
|
448 |
arcs.push_back(std::make_pair(1,2)); |
|
449 |
arcs.push_back(std::make_pair(3,0)); |
|
450 |
arcs.push_back(std::make_pair(3,3)); |
|
451 |
arcs.push_back(std::make_pair(4,2)); |
|
452 |
arcs.push_back(std::make_pair(4,3)); |
|
453 |
arcs.push_back(std::make_pair(4,1)); |
|
454 |
|
|
455 |
G.build(6, arcs.begin(), arcs.end()); |
|
456 |
|
|
457 |
checkGraphNodeList(G, 6); |
|
458 |
checkGraphArcList(G, 9); |
|
459 |
|
|
460 |
checkGraphOutArcList(G, G.node(0), 2); |
|
461 |
checkGraphOutArcList(G, G.node(1), 2); |
|
462 |
checkGraphOutArcList(G, G.node(2), 0); |
|
463 |
checkGraphOutArcList(G, G.node(3), 2); |
|
464 |
checkGraphOutArcList(G, G.node(4), 3); |
|
465 |
checkGraphOutArcList(G, G.node(5), 0); |
|
466 |
|
|
467 |
checkGraphInArcList(G, G.node(0), 1); |
|
468 |
checkGraphInArcList(G, G.node(1), 2); |
|
469 |
checkGraphInArcList(G, G.node(2), 3); |
|
470 |
checkGraphInArcList(G, G.node(3), 3); |
|
471 |
checkGraphInArcList(G, G.node(4), 0); |
|
472 |
checkGraphInArcList(G, G.node(5), 0); |
|
473 |
|
|
474 |
checkGraphConArcList(G, 9); |
|
475 |
|
|
444 | 476 |
checkNodeIds(G); |
445 | 477 |
checkArcIds(G); |
446 | 478 |
checkGraphNodeMap(G); |
447 | 479 |
checkGraphArcMap(G); |
448 | 480 |
|
449 | 481 |
int n = G.nodeNum(); |
450 | 482 |
int m = G.arcNum(); |
451 | 483 |
check(G.index(G.node(n-1)) == n-1, "Wrong index."); |
452 | 484 |
check(G.index(G.arc(m-1)) == m-1, "Wrong index."); |
453 | 485 |
} |
454 | 486 |
|
455 | 487 |
void checkFullDigraph(int num) { |
456 | 488 |
typedef FullDigraph Digraph; |
457 | 489 |
DIGRAPH_TYPEDEFS(Digraph); |
458 | 490 |
Digraph G(num); |
459 | 491 |
|
0 comments (0 inline)