diff -r 0286c68fc680 -r a72cac00e274 src/work/proba.cc --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/src/work/proba.cc Tue Jan 20 21:27:45 2004 +0000 @@ -0,0 +1,222 @@ +#include +#include +#include + +#include +#include +#include +#include +#include +#include + +using namespace marci; + + +int main (int, char*[]) +{ + typedef graph_traits::node_iterator node_iterator; + typedef graph_traits::edge_iterator edge_iterator; + typedef graph_traits::each_node_iterator each_node_iterator; + typedef graph_traits::each_edge_iterator each_edge_iterator; + typedef graph_traits::out_edge_iterator out_edge_iterator; + typedef graph_traits::in_edge_iterator in_edge_iterator; + typedef graph_traits::sym_edge_iterator sym_edge_iterator; + + list_graph flow_test; + + //Ahuja könyv példája, maxflowvalue=13 + node_iterator s=flow_test.add_node(); + node_iterator v1=flow_test.add_node(); + node_iterator v2=flow_test.add_node(); + node_iterator v3=flow_test.add_node(); + node_iterator v4=flow_test.add_node(); + node_iterator v5=flow_test.add_node(); + node_iterator t=flow_test.add_node(); + + node_property_vector node_name(flow_test); + node_name.put(s, "s"); + node_name.put(v1, "v1"); + node_name.put(v2, "v2"); + node_name.put(v3, "v3"); + node_name.put(v4, "v4"); + node_name.put(v5, "v5"); + node_name.put(t, "t"); + + edge_iterator s_v1=flow_test.add_edge(s, v1); + edge_iterator s_v2=flow_test.add_edge(s, v2); + edge_iterator s_v3=flow_test.add_edge(s, v3); + + edge_iterator v2_v4=flow_test.add_edge(v2, v4); + edge_iterator v2_v5=flow_test.add_edge(v2, v5); + + edge_iterator v3_v5=flow_test.add_edge(v3, v5); + + edge_iterator v4_t=flow_test.add_edge(v4, t); + edge_iterator v5_t=flow_test.add_edge(v5, t); + + edge_property_vector cap(flow_test); + cap.put(s_v1, 0); + cap.put(s_v2, 10); + cap.put(s_v3, 10); + cap.put(v2_v4, 5); + cap.put(v2_v5, 8); + cap.put(v3_v5, 5); + cap.put(v4_t, 8); + cap.put(v5_t, 8); + + + + //Marci példája, maxflowvalue=23 + /* node_iterator s=flow_test.add_node(); + node_iterator v1=flow_test.add_node(); + node_iterator v2=flow_test.add_node(); + node_iterator v3=flow_test.add_node(); + node_iterator v4=flow_test.add_node(); + node_iterator t=flow_test.add_node(); + node_iterator w=flow_test.add_node(); + + + node_property_vector node_name(flow_test); + node_name.put(s, "s"); + node_name.put(v1, "v1"); + node_name.put(v2, "v2"); + node_name.put(v3, "v3"); + node_name.put(v4, "v4"); + node_name.put(t, "t"); + node_name.put(w, "w"); + + edge_iterator s_v1=flow_test.add_edge(s, v1); + edge_iterator s_v2=flow_test.add_edge(s, v2); + edge_iterator v1_v2=flow_test.add_edge(v1, v2); + edge_iterator v2_v1=flow_test.add_edge(v2, v1); + edge_iterator v1_v3=flow_test.add_edge(v1, v3); + edge_iterator v3_v2=flow_test.add_edge(v3, v2); + edge_iterator v2_v4=flow_test.add_edge(v2, v4); + edge_iterator v4_v3=flow_test.add_edge(v4, v3); + edge_iterator v3_t=flow_test.add_edge(v3, t); + edge_iterator v4_t=flow_test.add_edge(v4, t); + edge_iterator v3_v3=flow_test.add_edge(v3, v3); + edge_iterator s_w=flow_test.add_edge(s, w); + + + + edge_property_vector cap(flow_test); + cap.put(s_v1, 16); + cap.put(s_v2, 13); + cap.put(v1_v2, 10); + cap.put(v2_v1, 4); + cap.put(v1_v3, 12); + cap.put(v3_v2, 9); + cap.put(v2_v4, 14); + cap.put(v4_v3, 7); + cap.put(v3_t, 20); + cap.put(v4_t, 4); + cap.put(v3_v3, 4); + cap.put(s_w, 4); + */ + + + //pelda 3, maxflowvalue=4 + /* node_iterator s=flow_test.add_node(); + node_iterator v1=flow_test.add_node(); + node_iterator v2=flow_test.add_node(); + node_iterator t=flow_test.add_node(); + node_iterator w=flow_test.add_node(); + + node_property_vector node_name(flow_test); + node_name.put(s, "s"); + node_name.put(v1, "v1"); + node_name.put(v2, "v2"); + node_name.put(t, "t"); + node_name.put(w, "w"); + + edge_iterator s_v1=flow_test.add_edge(s, v1); + edge_iterator v1_v2=flow_test.add_edge(v1, v2); + edge_iterator v2_t=flow_test.add_edge(v2, t); + edge_iterator v1_v1=flow_test.add_edge(v1, v1); + edge_iterator s_w=flow_test.add_edge(s, w); + + + edge_property_vector cap(flow_test); + + cap.put(s_v1, 16); + cap.put(v1_v2, 10); + cap.put(v2_t, 4); + cap.put(v1_v1, 3); + cap.put(s_w, 5); + + */ + + + + std::cout << "Testing reverse_bfs..." << std::endl; + + reverse_bfs bfs_test(flow_test, t); + + bfs_test.run(); + + for (each_node_iterator w=flow_test.first_node(); w.is_valid(); ++w) { + std::cout <<"The distance of " << w << " is " << bfs_test.dist(w) < preflow_push_test(flow_test, s, t, cap); + + preflow_push_test.run(); + + std::cout << "Maximum flow value is: " << preflow_push_test.maxflow() << "."< flow=preflow_push_test.allflow(); + for (each_edge_iterator e=flow_test.first_edge(); e.is_valid(); ++e) { + std::cout <<"Flow on edge " << flow_test.tail(e) <<"-" << flow_test.head(e)<< " is " < mincut=preflow_push_test.mincut(); + + for (each_node_iterator v=flow_test.first_node(); v.is_valid(); ++v) { + if (mincut.get(v)) std::cout < max_flow_test(flow_test, s, t, cap); + + max_flow_test.run(); + + std::cout << "Maximum flow value is: " << max_flow_test.maxflow() << "."<< std::endl; + + std::cout << "A minimum cut: " < mincut2=max_flow_test.mincut(); + + for (each_node_iterator v=flow_test.first_node(); v.is_valid(); ++v) { + if (mincut2.get(v)) std::cout <