alpar@906
|
1 |
/* -*- C++ -*-
|
alpar@921
|
2 |
* src/test/min_cost_flow_test.cc - Part of LEMON, a generic C++ optimization library
|
alpar@906
|
3 |
*
|
alpar@906
|
4 |
* Copyright (C) 2004 Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
|
alpar@906
|
5 |
* (Egervary Combinatorial Optimization Research Group, EGRES).
|
alpar@906
|
6 |
*
|
alpar@906
|
7 |
* Permission to use, modify and distribute this software is granted
|
alpar@906
|
8 |
* provided that this copyright notice appears in all copies. For
|
alpar@906
|
9 |
* precise terms see the accompanying LICENSE file.
|
alpar@906
|
10 |
*
|
alpar@906
|
11 |
* This software is provided "AS IS" with no warranty of any kind,
|
alpar@906
|
12 |
* express or implied, and with no claim as to its suitability for any
|
alpar@906
|
13 |
* purpose.
|
alpar@906
|
14 |
*
|
alpar@906
|
15 |
*/
|
alpar@906
|
16 |
|
alpar@899
|
17 |
#include <iostream>
|
alpar@899
|
18 |
#include "test_tools.h"
|
alpar@921
|
19 |
#include <lemon/list_graph.h>
|
alpar@921
|
20 |
#include <lemon/min_cost_flow.h>
|
alpar@899
|
21 |
//#include <path.h>
|
alpar@899
|
22 |
//#include <maps.h>
|
alpar@899
|
23 |
|
alpar@899
|
24 |
using namespace std;
|
alpar@921
|
25 |
using namespace lemon;
|
alpar@899
|
26 |
|
alpar@899
|
27 |
|
alpar@899
|
28 |
|
alpar@899
|
29 |
bool passed = true;
|
alpar@899
|
30 |
/*
|
alpar@899
|
31 |
void check(bool rc, char *msg="") {
|
alpar@899
|
32 |
passed = passed && rc;
|
alpar@899
|
33 |
if(!rc) {
|
alpar@899
|
34 |
std::cerr << "Test failed! ("<< msg << ")" << std::endl; \
|
alpar@899
|
35 |
|
alpar@899
|
36 |
|
alpar@899
|
37 |
}
|
alpar@899
|
38 |
}
|
alpar@899
|
39 |
*/
|
alpar@899
|
40 |
|
alpar@899
|
41 |
|
alpar@899
|
42 |
int main()
|
alpar@899
|
43 |
{
|
alpar@899
|
44 |
|
alpar@899
|
45 |
typedef ListGraph::Node Node;
|
alpar@899
|
46 |
typedef ListGraph::Edge Edge;
|
alpar@899
|
47 |
|
alpar@899
|
48 |
ListGraph graph;
|
alpar@899
|
49 |
|
alpar@899
|
50 |
//Ahuja könyv példája
|
alpar@899
|
51 |
|
alpar@899
|
52 |
Node s=graph.addNode();
|
alpar@899
|
53 |
Node v1=graph.addNode();
|
alpar@899
|
54 |
Node v2=graph.addNode();
|
alpar@899
|
55 |
Node v3=graph.addNode();
|
alpar@899
|
56 |
Node v4=graph.addNode();
|
alpar@899
|
57 |
Node v5=graph.addNode();
|
alpar@899
|
58 |
Node t=graph.addNode();
|
alpar@899
|
59 |
|
alpar@899
|
60 |
Edge s_v1=graph.addEdge(s, v1);
|
alpar@899
|
61 |
Edge v1_v2=graph.addEdge(v1, v2);
|
alpar@899
|
62 |
Edge s_v3=graph.addEdge(s, v3);
|
alpar@899
|
63 |
Edge v2_v4=graph.addEdge(v2, v4);
|
alpar@899
|
64 |
Edge v2_v5=graph.addEdge(v2, v5);
|
alpar@899
|
65 |
Edge v3_v5=graph.addEdge(v3, v5);
|
alpar@899
|
66 |
Edge v4_t=graph.addEdge(v4, t);
|
alpar@899
|
67 |
Edge v5_t=graph.addEdge(v5, t);
|
alpar@899
|
68 |
|
alpar@899
|
69 |
|
alpar@899
|
70 |
ListGraph::EdgeMap<int> length(graph);
|
alpar@899
|
71 |
|
alpar@899
|
72 |
length.set(s_v1, 6);
|
alpar@899
|
73 |
length.set(v1_v2, 4);
|
alpar@899
|
74 |
length.set(s_v3, 10);
|
alpar@899
|
75 |
length.set(v2_v4, 5);
|
alpar@899
|
76 |
length.set(v2_v5, 1);
|
alpar@899
|
77 |
length.set(v3_v5, 4);
|
alpar@899
|
78 |
length.set(v4_t, 8);
|
alpar@899
|
79 |
length.set(v5_t, 8);
|
alpar@899
|
80 |
|
alpar@899
|
81 |
ListGraph::EdgeMap<int> capacity(graph);
|
alpar@899
|
82 |
|
alpar@899
|
83 |
capacity.set(s_v1, 2);
|
alpar@899
|
84 |
capacity.set(v1_v2, 2);
|
alpar@899
|
85 |
capacity.set(s_v3, 1);
|
alpar@899
|
86 |
capacity.set(v2_v4, 1);
|
alpar@899
|
87 |
capacity.set(v2_v5, 1);
|
alpar@899
|
88 |
capacity.set(v3_v5, 1);
|
alpar@899
|
89 |
capacity.set(v4_t, 1);
|
alpar@899
|
90 |
capacity.set(v5_t, 2);
|
alpar@899
|
91 |
|
alpar@899
|
92 |
// ConstMap<Edge, int> const1map(1);
|
alpar@899
|
93 |
std::cout << "Mincostflows algorithm test..." << std::endl;
|
alpar@899
|
94 |
|
alpar@899
|
95 |
MinCostFlow< ListGraph, ListGraph::EdgeMap<int>, ListGraph::EdgeMap<int> >
|
alpar@899
|
96 |
surb_test(graph, length, capacity);
|
alpar@899
|
97 |
|
alpar@899
|
98 |
int k=1;
|
alpar@899
|
99 |
|
alpar@899
|
100 |
check( surb_test.run(s,t,k) == 1 && surb_test.totalLength() == 19,"One path, total length should be 19");
|
alpar@899
|
101 |
|
alpar@899
|
102 |
check(surb_test.checkComplementarySlackness(), "Is the primal-dual solution pair really optimal?");
|
alpar@899
|
103 |
|
alpar@899
|
104 |
k=2;
|
alpar@899
|
105 |
|
alpar@899
|
106 |
check( surb_test.run(s,t,k) == 2 && surb_test.totalLength() == 41,"Two paths, total length should be 41");
|
alpar@899
|
107 |
|
alpar@899
|
108 |
check(surb_test.checkComplementarySlackness(), "Is the primal-dual solution pair really optimal?");
|
alpar@899
|
109 |
|
alpar@899
|
110 |
|
alpar@899
|
111 |
k=4;
|
alpar@899
|
112 |
|
alpar@899
|
113 |
check( surb_test.run(s,t,k) == 3 && surb_test.totalLength() == 64,"Three paths, total length should be 64");
|
alpar@899
|
114 |
|
alpar@899
|
115 |
check(surb_test.checkComplementarySlackness(), "Is the primal-dual solution pair really optimal?");
|
alpar@899
|
116 |
|
alpar@899
|
117 |
|
alpar@899
|
118 |
cout << (passed ? "All tests passed." : "Some of the tests failed!!!")
|
alpar@899
|
119 |
<< endl;
|
alpar@899
|
120 |
|
alpar@899
|
121 |
return passed ? 0 : 1;
|
alpar@899
|
122 |
|
alpar@899
|
123 |
}
|