| [209] | 1 | /* -*- mode: C++; indent-tabs-mode: nil; -*- | 
|---|
| [103] | 2 | * | 
|---|
| [209] | 3 | * This file is a part of LEMON, a generic C++ optimization library. | 
|---|
| [103] | 4 | * | 
|---|
| [440] | 5 | * Copyright (C) 2003-2009 | 
|---|
| [103] | 6 | * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport | 
|---|
|  | 7 | * (Egervary Research Group on Combinatorial Optimization, EGRES). | 
|---|
|  | 8 | * | 
|---|
|  | 9 | * Permission to use, modify and distribute this software is granted | 
|---|
|  | 10 | * provided that this copyright notice appears in all copies. For | 
|---|
|  | 11 | * precise terms see the accompanying LICENSE file. | 
|---|
|  | 12 | * | 
|---|
|  | 13 | * This software is provided "AS IS" with no warranty of any kind, | 
|---|
|  | 14 | * express or implied, and with no claim as to its suitability for any | 
|---|
|  | 15 | * purpose. | 
|---|
|  | 16 | * | 
|---|
|  | 17 | */ | 
|---|
|  | 18 |  | 
|---|
|  | 19 | #include <iostream> | 
|---|
|  | 20 | #include <vector> | 
|---|
|  | 21 |  | 
|---|
|  | 22 | #include "test_tools.h" | 
|---|
|  | 23 | #include <lemon/maps.h> | 
|---|
|  | 24 | #include <lemon/kruskal.h> | 
|---|
|  | 25 | #include <lemon/list_graph.h> | 
|---|
|  | 26 |  | 
|---|
|  | 27 | #include <lemon/concepts/maps.h> | 
|---|
|  | 28 | #include <lemon/concepts/digraph.h> | 
|---|
|  | 29 | #include <lemon/concepts/graph.h> | 
|---|
|  | 30 |  | 
|---|
|  | 31 | using namespace std; | 
|---|
|  | 32 | using namespace lemon; | 
|---|
|  | 33 |  | 
|---|
|  | 34 | void checkCompileKruskal() | 
|---|
|  | 35 | { | 
|---|
|  | 36 | concepts::WriteMap<concepts::Digraph::Arc,bool> w; | 
|---|
|  | 37 | concepts::WriteMap<concepts::Graph::Edge,bool> uw; | 
|---|
|  | 38 |  | 
|---|
|  | 39 | concepts::ReadMap<concepts::Digraph::Arc,int> r; | 
|---|
|  | 40 | concepts::ReadMap<concepts::Graph::Edge,int> ur; | 
|---|
|  | 41 |  | 
|---|
|  | 42 | concepts::Digraph g; | 
|---|
|  | 43 | concepts::Graph ug; | 
|---|
|  | 44 |  | 
|---|
|  | 45 | kruskal(g, r, w); | 
|---|
|  | 46 | kruskal(ug, ur, uw); | 
|---|
|  | 47 |  | 
|---|
|  | 48 | std::vector<std::pair<concepts::Digraph::Arc, int> > rs; | 
|---|
|  | 49 | std::vector<std::pair<concepts::Graph::Edge, int> > urs; | 
|---|
|  | 50 |  | 
|---|
|  | 51 | kruskal(g, rs, w); | 
|---|
|  | 52 | kruskal(ug, urs, uw); | 
|---|
|  | 53 |  | 
|---|
|  | 54 | std::vector<concepts::Digraph::Arc> ws; | 
|---|
|  | 55 | std::vector<concepts::Graph::Edge> uws; | 
|---|
|  | 56 |  | 
|---|
|  | 57 | kruskal(g, r, ws.begin()); | 
|---|
|  | 58 | kruskal(ug, ur, uws.begin()); | 
|---|
|  | 59 | } | 
|---|
|  | 60 |  | 
|---|
|  | 61 | int main() { | 
|---|
|  | 62 |  | 
|---|
|  | 63 | typedef ListGraph::Node Node; | 
|---|
|  | 64 | typedef ListGraph::Edge Edge; | 
|---|
|  | 65 | typedef ListGraph::NodeIt NodeIt; | 
|---|
|  | 66 | typedef ListGraph::ArcIt ArcIt; | 
|---|
|  | 67 |  | 
|---|
|  | 68 | ListGraph G; | 
|---|
|  | 69 |  | 
|---|
|  | 70 | Node s=G.addNode(); | 
|---|
|  | 71 | Node v1=G.addNode(); | 
|---|
|  | 72 | Node v2=G.addNode(); | 
|---|
|  | 73 | Node v3=G.addNode(); | 
|---|
|  | 74 | Node v4=G.addNode(); | 
|---|
|  | 75 | Node t=G.addNode(); | 
|---|
| [209] | 76 |  | 
|---|
| [103] | 77 | Edge e1 = G.addEdge(s, v1); | 
|---|
|  | 78 | Edge e2 = G.addEdge(s, v2); | 
|---|
|  | 79 | Edge e3 = G.addEdge(v1, v2); | 
|---|
|  | 80 | Edge e4 = G.addEdge(v2, v1); | 
|---|
|  | 81 | Edge e5 = G.addEdge(v1, v3); | 
|---|
|  | 82 | Edge e6 = G.addEdge(v3, v2); | 
|---|
|  | 83 | Edge e7 = G.addEdge(v2, v4); | 
|---|
|  | 84 | Edge e8 = G.addEdge(v4, v3); | 
|---|
|  | 85 | Edge e9 = G.addEdge(v3, t); | 
|---|
|  | 86 | Edge e10 = G.addEdge(v4, t); | 
|---|
|  | 87 |  | 
|---|
|  | 88 | typedef ListGraph::EdgeMap<int> ECostMap; | 
|---|
|  | 89 | typedef ListGraph::EdgeMap<bool> EBoolMap; | 
|---|
|  | 90 |  | 
|---|
|  | 91 | ECostMap edge_cost_map(G, 2); | 
|---|
|  | 92 | EBoolMap tree_map(G); | 
|---|
| [209] | 93 |  | 
|---|
| [103] | 94 |  | 
|---|
|  | 95 | //Test with const map. | 
|---|
|  | 96 | check(kruskal(G, ConstMap<ListGraph::Edge,int>(2), tree_map)==10, | 
|---|
| [209] | 97 | "Total cost should be 10"); | 
|---|
| [171] | 98 | //Test with an edge map (filled with uniform costs). | 
|---|
| [103] | 99 | check(kruskal(G, edge_cost_map, tree_map)==10, | 
|---|
| [209] | 100 | "Total cost should be 10"); | 
|---|
| [103] | 101 |  | 
|---|
| [581] | 102 | edge_cost_map[e1] = -10; | 
|---|
|  | 103 | edge_cost_map[e2] = -9; | 
|---|
|  | 104 | edge_cost_map[e3] = -8; | 
|---|
|  | 105 | edge_cost_map[e4] = -7; | 
|---|
|  | 106 | edge_cost_map[e5] = -6; | 
|---|
|  | 107 | edge_cost_map[e6] = -5; | 
|---|
|  | 108 | edge_cost_map[e7] = -4; | 
|---|
|  | 109 | edge_cost_map[e8] = -3; | 
|---|
|  | 110 | edge_cost_map[e9] = -2; | 
|---|
|  | 111 | edge_cost_map[e10] = -1; | 
|---|
| [103] | 112 |  | 
|---|
|  | 113 | vector<Edge> tree_edge_vec(5); | 
|---|
|  | 114 |  | 
|---|
|  | 115 | //Test with a edge map and inserter. | 
|---|
|  | 116 | check(kruskal(G, edge_cost_map, | 
|---|
| [209] | 117 | tree_edge_vec.begin()) | 
|---|
|  | 118 | ==-31, | 
|---|
|  | 119 | "Total cost should be -31."); | 
|---|
|  | 120 |  | 
|---|
| [103] | 121 | tree_edge_vec.clear(); | 
|---|
|  | 122 |  | 
|---|
|  | 123 | check(kruskal(G, edge_cost_map, | 
|---|
| [209] | 124 | back_inserter(tree_edge_vec)) | 
|---|
|  | 125 | ==-31, | 
|---|
|  | 126 | "Total cost should be -31."); | 
|---|
|  | 127 |  | 
|---|
| [103] | 128 | //   tree_edge_vec.clear(); | 
|---|
| [209] | 129 |  | 
|---|
| [103] | 130 | //   //The above test could also be coded like this: | 
|---|
|  | 131 | //   check(kruskal(G, | 
|---|
| [209] | 132 | //                 makeKruskalMapInput(G, edge_cost_map), | 
|---|
|  | 133 | //                 makeKruskalSequenceOutput(back_inserter(tree_edge_vec))) | 
|---|
|  | 134 | //         ==-31, | 
|---|
|  | 135 | //         "Total cost should be -31."); | 
|---|
| [103] | 136 |  | 
|---|
|  | 137 | check(tree_edge_vec.size()==5,"The tree should have 5 edges."); | 
|---|
|  | 138 |  | 
|---|
|  | 139 | check(tree_edge_vec[0]==e1 && | 
|---|
| [209] | 140 | tree_edge_vec[1]==e2 && | 
|---|
|  | 141 | tree_edge_vec[2]==e5 && | 
|---|
|  | 142 | tree_edge_vec[3]==e7 && | 
|---|
|  | 143 | tree_edge_vec[4]==e9, | 
|---|
|  | 144 | "Wrong tree."); | 
|---|
| [103] | 145 |  | 
|---|
|  | 146 | return 0; | 
|---|
|  | 147 | } | 
|---|