[1182] | 1 | #include <iostream> |
---|
| 2 | #include <vector> |
---|
| 3 | |
---|
| 4 | #include <lemon/maps.h> |
---|
| 5 | #include <lemon/kruskal.h> |
---|
| 6 | #include <lemon/list_graph.h> |
---|
| 7 | |
---|
| 8 | |
---|
| 9 | using namespace std; |
---|
| 10 | using namespace lemon; |
---|
| 11 | |
---|
| 12 | |
---|
| 13 | int main() { |
---|
| 14 | |
---|
| 15 | typedef ListGraph::Node Node; |
---|
| 16 | typedef ListGraph::Edge Edge; |
---|
| 17 | typedef ListGraph::NodeIt NodeIt; |
---|
| 18 | typedef ListGraph::EdgeIt EdgeIt; |
---|
| 19 | |
---|
| 20 | ListGraph G; |
---|
| 21 | |
---|
| 22 | Node s=G.addNode(); |
---|
| 23 | Node v1=G.addNode(); |
---|
| 24 | Node v2=G.addNode(); |
---|
| 25 | Node v3=G.addNode(); |
---|
| 26 | Node v4=G.addNode(); |
---|
| 27 | Node t=G.addNode(); |
---|
| 28 | |
---|
| 29 | Edge e1 = G.addEdge(s, v1); |
---|
| 30 | Edge e2 = G.addEdge(s, v2); |
---|
| 31 | Edge e3 = G.addEdge(v1, v2); |
---|
| 32 | Edge e4 = G.addEdge(v2, v1); |
---|
| 33 | Edge e5 = G.addEdge(v1, v3); |
---|
| 34 | Edge e6 = G.addEdge(v3, v2); |
---|
| 35 | Edge e7 = G.addEdge(v2, v4); |
---|
| 36 | Edge e8 = G.addEdge(v4, v3); |
---|
| 37 | Edge e9 = G.addEdge(v3, t); |
---|
| 38 | Edge e10 = G.addEdge(v4, t); |
---|
| 39 | |
---|
| 40 | typedef ListGraph::EdgeMap<int> ECostMap; |
---|
| 41 | typedef ListGraph::EdgeMap<bool> EBoolMap; |
---|
| 42 | |
---|
| 43 | ECostMap edge_cost_map(G, 2); |
---|
| 44 | EBoolMap tree_map(G); |
---|
| 45 | |
---|
| 46 | |
---|
| 47 | //Test with const map. |
---|
| 48 | std::cout << "The weight of the minimum spanning tree is " << kruskalEdgeMap(G, ConstMap<ListGraph::Edge,int>(2), tree_map)<<std::endl; |
---|
| 49 | |
---|
| 50 | /* |
---|
| 51 | ==10, |
---|
| 52 | "Total cost should be 10"); |
---|
| 53 | //Test with a edge map (filled with uniform costs). |
---|
| 54 | check(kruskalEdgeMap(G, edge_cost_map, tree_map)==10, |
---|
| 55 | "Total cost should be 10"); |
---|
| 56 | |
---|
| 57 | edge_cost_map.set(e1, -10); |
---|
| 58 | edge_cost_map.set(e2, -9); |
---|
| 59 | edge_cost_map.set(e3, -8); |
---|
| 60 | edge_cost_map.set(e4, -7); |
---|
| 61 | edge_cost_map.set(e5, -6); |
---|
| 62 | edge_cost_map.set(e6, -5); |
---|
| 63 | edge_cost_map.set(e7, -4); |
---|
| 64 | edge_cost_map.set(e8, -3); |
---|
| 65 | edge_cost_map.set(e9, -2); |
---|
| 66 | edge_cost_map.set(e10, -1); |
---|
| 67 | |
---|
| 68 | vector<Edge> tree_edge_vec; |
---|
| 69 | |
---|
| 70 | //Test with a edge map and inserter. |
---|
| 71 | check(kruskalEdgeMap_IteratorOut(G, edge_cost_map, |
---|
| 72 | back_inserter(tree_edge_vec)) |
---|
| 73 | ==-31, |
---|
| 74 | "Total cost should be -31."); |
---|
| 75 | |
---|
| 76 | tree_edge_vec.clear(); |
---|
| 77 | |
---|
| 78 | //The above test could also be coded like this: |
---|
| 79 | check(kruskal(G, |
---|
| 80 | makeKruskalMapInput(G, edge_cost_map), |
---|
| 81 | makeKruskalSequenceOutput(back_inserter(tree_edge_vec))) |
---|
| 82 | ==-31, |
---|
| 83 | "Total cost should be -31."); |
---|
| 84 | |
---|
| 85 | check(tree_edge_vec.size()==5,"The tree should have 5 edges."); |
---|
| 86 | |
---|
| 87 | check(tree_edge_vec[0]==e1 && |
---|
| 88 | tree_edge_vec[1]==e2 && |
---|
| 89 | tree_edge_vec[2]==e5 && |
---|
| 90 | tree_edge_vec[3]==e7 && |
---|
| 91 | tree_edge_vec[4]==e9, |
---|
| 92 | "Wrong tree."); |
---|
| 93 | */ |
---|
| 94 | return 0; |
---|
| 95 | } |
---|