Example
#include <iostream>
#include <seqan/graph_algorithms.h>
using namespace seqan;
int main()
{
typedef Graph<Directed<> > TGraph;
typedef VertexDescriptor<TGraph>::Type TVertexDescriptor;
typedef EdgeDescriptor<TGraph>::Type TEdgeDescriptor;
typedef Size<TGraph>::Type TSize;
// Create graph with 5 directed edges (3,0), (1,2), ...
TSize numEdges = 5;
TVertexDescriptor edges[] = {3, 0, 1, 2, 2, 1, 1, 3, 3, 2};
TGraph g;
addEdges(g, edges, numEdges);
// Print graph to stdout.
std::cout << g << "\n";
// Compute transitive closure.
String<bool> closure;
transitiveClosure(closure, g);
// Print result to stdout.
TSize len = static_cast<TSize>(std::sqrt((double) length(closure)));
for (TSize row = 0; row < len; ++row)
{
for (TSize col = 0; col < len; ++col)
std::cout << getValue(closure, row * len + col) << ",";
std::cout << std::endl;
}
return 0;
}
Adjacency list:
0 ->
1 -> 3,2,
2 -> 1,
3 -> 2,0,
Edge list:
Source: 1,Target: 3 (Id: 3)
Source: 1,Target: 2 (Id: 1)
Source: 2,Target: 1 (Id: 2)
Source: 3,Target: 2 (Id: 4)
Source: 3,Target: 0 (Id: 0)
1,0,0,0,
1,1,1,1,
1,1,1,1,
1,1,1,1,