COIN-OR::LEMON - Graph Library

source: lemon-0.x/src/work/proba.cc @ 50:e125f12784e2

Last change on this file since 50:e125f12784e2 was 30:10a3f2e0928c, checked in by jacint, 20 years ago

is_valid changed to valid

File size: 6.5 KB
Line 
1#include <iostream>
2#include <vector>
3#include <string>
4
5#include <marci_list_graph.hh>
6#include <marci_graph_traits.hh>
7#include <marci_property_vector.hh>
8#include <preflow_push_hl.hh>
9#include <preflow_push_max_flow.hh>
10#include <reverse_bfs.hh>
11//#include <dijkstra.hh>
12
13using namespace marci;
14
15
16int main (int, char*[])
17{
18  typedef graph_traits<list_graph>::node_iterator node_iterator;
19  typedef graph_traits<list_graph>::edge_iterator edge_iterator;
20  typedef graph_traits<list_graph>::each_node_iterator each_node_iterator;
21  typedef graph_traits<list_graph>::each_edge_iterator each_edge_iterator;
22  typedef graph_traits<list_graph>::out_edge_iterator out_edge_iterator;
23  typedef graph_traits<list_graph>::in_edge_iterator in_edge_iterator;
24  typedef graph_traits<list_graph>::sym_edge_iterator sym_edge_iterator;
25
26  list_graph flow_test;
27 
28  /*  //Ahuja könyv példája, maxflowvalue=13
29  node_iterator s=flow_test.add_node();
30  node_iterator v1=flow_test.add_node();
31  node_iterator v2=flow_test.add_node();
32  node_iterator v3=flow_test.add_node();
33  node_iterator v4=flow_test.add_node();
34  node_iterator v5=flow_test.add_node();
35  node_iterator t=flow_test.add_node();
36 
37  node_property_vector<list_graph, std::string> node_name(flow_test);
38  node_name.put(s, "s");
39  node_name.put(v1, "v1");
40  node_name.put(v2, "v2");
41  node_name.put(v3, "v3");
42  node_name.put(v4, "v4");
43  node_name.put(v5, "v5");
44  node_name.put(t, "t");
45
46  edge_iterator s_v1=flow_test.add_edge(s, v1);
47  edge_iterator s_v2=flow_test.add_edge(s, v2);
48  edge_iterator s_v3=flow_test.add_edge(s, v3);
49  edge_iterator v2_v4=flow_test.add_edge(v2, v4);
50  edge_iterator v2_v5=flow_test.add_edge(v2, v5);
51  edge_iterator v3_v5=flow_test.add_edge(v3, v5);
52  edge_iterator v4_t=flow_test.add_edge(v4, t);
53  edge_iterator v5_t=flow_test.add_edge(v5, t);
54  edge_iterator v2_s=flow_test.add_edge(v2, s);
55 
56  edge_property_vector<list_graph, int> cap(flow_test); 
57  cap.put(s_v1, 0);
58  cap.put(s_v2, 10);
59  cap.put(s_v3, 10);
60  cap.put(v2_v4, 5);
61  cap.put(v2_v5, 8);
62  cap.put(v3_v5, 5);
63  cap.put(v4_t, 8);
64  cap.put(v5_t, 8);
65  cap.put(v2_s, 0);
66*/
67
68 
69  //Marci példája, maxflowvalue=23
70     node_iterator s=flow_test.add_node();
71  node_iterator v1=flow_test.add_node();
72  node_iterator v2=flow_test.add_node();
73  node_iterator v3=flow_test.add_node();
74  node_iterator v4=flow_test.add_node();
75  node_iterator t=flow_test.add_node();
76  node_iterator w=flow_test.add_node();
77
78 
79  node_property_vector<list_graph, std::string> node_name(flow_test);
80  node_name.put(s, "s");
81  node_name.put(v1, "v1");
82  node_name.put(v2, "v2");
83  node_name.put(v3, "v3");
84  node_name.put(v4, "v4");
85  node_name.put(t, "t");
86  node_name.put(w, "w");
87
88  edge_iterator s_v1=flow_test.add_edge(s, v1);
89  edge_iterator s_v2=flow_test.add_edge(s, v2);
90  edge_iterator v1_v2=flow_test.add_edge(v1, v2);
91  edge_iterator v2_v1=flow_test.add_edge(v2, v1);
92  edge_iterator v1_v3=flow_test.add_edge(v1, v3);
93  edge_iterator v3_v2=flow_test.add_edge(v3, v2);
94  edge_iterator v2_v4=flow_test.add_edge(v2, v4);
95  edge_iterator v4_v3=flow_test.add_edge(v4, v3);
96  edge_iterator v3_t=flow_test.add_edge(v3, t);
97  edge_iterator v4_t=flow_test.add_edge(v4, t);
98  edge_iterator v3_v3=flow_test.add_edge(v3, v3);
99  edge_iterator s_w=flow_test.add_edge(s, w);
100  edge_iterator v2_s=flow_test.add_edge(v2, s);
101 
102
103
104  edge_property_vector<list_graph, int> cap(flow_test); 
105  cap.put(s_v1, 16);
106  cap.put(s_v2, 13);
107  cap.put(v1_v2, 10);
108  cap.put(v2_v1, 4);
109  cap.put(v1_v3, 12);
110  cap.put(v3_v2, 9);
111  cap.put(v2_v4, 14);
112  cap.put(v4_v3, 7);
113  cap.put(v3_t, 20);
114  cap.put(v4_t, 4);
115  cap.put(v3_v3, 4);
116  cap.put(s_w, 4);
117  cap.put(v2_s, 0);
118
119
120  //pelda 3, maxflowvalue=4
121  /*    node_iterator s=flow_test.add_node();
122  node_iterator v1=flow_test.add_node();
123  node_iterator v2=flow_test.add_node();
124  node_iterator t=flow_test.add_node();
125  node_iterator w=flow_test.add_node();
126 
127  node_property_vector<list_graph, std::string> node_name(flow_test);
128  node_name.put(s, "s");
129  node_name.put(v1, "v1");
130  node_name.put(v2, "v2");
131  node_name.put(t, "t");
132  node_name.put(w, "w");
133
134  edge_iterator s_v1=flow_test.add_edge(s, v1);
135  edge_iterator v1_v2=flow_test.add_edge(v1, v2);
136  edge_iterator v2_t=flow_test.add_edge(v2, t);
137  edge_iterator v1_v1=flow_test.add_edge(v1, v1);
138  edge_iterator s_w=flow_test.add_edge(s, w);
139
140
141  edge_property_vector<list_graph, int> cap(flow_test);
142   
143  cap.put(s_v1, 16);
144  cap.put(v1_v2, 10);
145  cap.put(v2_t, 4);
146  cap.put(v1_v1, 3);
147  cap.put(s_w, 5);
148 
149  */
150
151
152
153  std::cout << "Testing reverse_bfs..." << std::endl;
154 
155  reverse_bfs<list_graph> bfs_test(flow_test, t);
156
157  bfs_test.run();
158
159  for (each_node_iterator w=flow_test.first_node(); w.valid(); ++w) {
160    std::cout <<"The distance of " << w << " is " << bfs_test.dist(w) <<std::endl;
161    }
162
163
164
165
166
167  std::cout << "Testing preflow_push_hl..." << std::endl;
168 
169  preflow_push_hl<list_graph, int> preflow_push_test(flow_test, s, t, cap);
170
171  preflow_push_test.run();
172
173  std::cout << "Maximum flow value is: " << preflow_push_test.maxflow() << "."<<std::endl;
174
175  std::cout<< "The flow on edge s-v1 is "<< preflow_push_test.flowonedge(s_v1) << "."<<std::endl;
176
177  edge_property_vector<list_graph, int> flow=preflow_push_test.allflow(); 
178  for (each_edge_iterator e=flow_test.first_edge(); e.valid(); ++e) {
179    std::cout <<"Flow on edge " << flow_test.tail(e) <<"-" << flow_test.head(e)<< " is " <<flow.get(e) <<std::endl;
180    }
181
182  std::cout << "A minimum cut: " <<std::endl; 
183  node_property_vector<list_graph, bool> mincut=preflow_push_test.mincut();
184
185  for (each_node_iterator v=flow_test.first_node(); v.valid(); ++v) {
186      if (mincut.get(v)) std::cout <<node_name.get(v)<< " ";
187    }
188 
189  std::cout<<"\n\n"<<std::endl;
190
191
192
193
194  std::cout << "Testing preflow_push_max_flow..." << std::endl;
195 
196  preflow_push_max_flow<list_graph, int> max_flow_test(flow_test, s, t, cap);
197
198  max_flow_test.run();
199
200  std::cout << "Maximum flow value is: " << max_flow_test.maxflow() << "."<< std::endl;
201
202  std::cout << "A minimum cut: " <<std::endl; 
203  node_property_vector<list_graph, bool> mincut2=max_flow_test.mincut();
204
205  for (each_node_iterator v=flow_test.first_node(); v.valid(); ++v) {
206    if (mincut2.get(v)) std::cout <<node_name.get(v)<< " ";
207  }
208 
209  std::cout << std::endl;
210
211
212
213//    std::cout << "Testing dijkstra..." << std::endl;
214 
215//    dijkstra<list_graph, int> dijkstra_test(flow_test, s, cap);
216
217//    dijkstra_test.run();
218
219//    for (each_node_iterator w=flow_test.first_node(); w.valid(); ++w) {
220//      std::cout <<"The distance of " << w << " is " << dijkstra_test.dist(w) <<std::endl;
221//      }
222
223
224
225
226  return 0;
227}
228
229
230
231
232
233
234
235
236
Note: See TracBrowser for help on using the repository browser.