alpar@896: #include <iostream>
alpar@896: #include "test_tools.h"
alpar@896: #include <hugo/list_graph.h>
alpar@896: #include <hugo/min_cost_flows.h>
alpar@896: //#include <path.h>
alpar@896: //#include <maps.h>
alpar@896: 
alpar@896: using namespace std;
alpar@896: using namespace hugo;
alpar@896: 
alpar@896: 
alpar@896: 
alpar@896: bool passed = true;
alpar@896: /*
alpar@896: void check(bool rc, char *msg="") {
alpar@896:   passed = passed && rc;
alpar@896:   if(!rc) {
alpar@896:     std::cerr << "Test failed! ("<< msg << ")" << std::endl; \
alpar@896:  
alpar@896: 
alpar@896:   }
alpar@896: }
alpar@896: */
alpar@896: 
alpar@896: 
alpar@896: int main()
alpar@896: {
alpar@896: 
alpar@896:   typedef ListGraph::Node Node;
alpar@896:   typedef ListGraph::Edge Edge;
alpar@896: 
alpar@896:   ListGraph graph;
alpar@896: 
alpar@896:   //Ahuja könyv példája
alpar@896: 
alpar@896:   Node s=graph.addNode();
alpar@896:   Node v1=graph.addNode();  
alpar@896:   Node v2=graph.addNode();
alpar@896:   Node v3=graph.addNode();
alpar@896:   Node v4=graph.addNode();
alpar@896:   Node v5=graph.addNode();
alpar@896:   Node t=graph.addNode();
alpar@896: 
alpar@896:   Edge s_v1=graph.addEdge(s, v1);
alpar@896:   Edge v1_v2=graph.addEdge(v1, v2);
alpar@896:   Edge s_v3=graph.addEdge(s, v3);
alpar@896:   Edge v2_v4=graph.addEdge(v2, v4);
alpar@896:   Edge v2_v5=graph.addEdge(v2, v5);
alpar@896:   Edge v3_v5=graph.addEdge(v3, v5);
alpar@896:   Edge v4_t=graph.addEdge(v4, t);
alpar@896:   Edge v5_t=graph.addEdge(v5, t);
alpar@896:   
alpar@896: 
alpar@896:   ListGraph::EdgeMap<int> length(graph);
alpar@896: 
alpar@896:   length.set(s_v1, 6);
alpar@896:   length.set(v1_v2, 4);
alpar@896:   length.set(s_v3, 10);
alpar@896:   length.set(v2_v4, 5);
alpar@896:   length.set(v2_v5, 1);
alpar@896:   length.set(v3_v5, 4);
alpar@896:   length.set(v4_t, 8);
alpar@896:   length.set(v5_t, 8);
alpar@896: 
alpar@896:   ListGraph::EdgeMap<int> capacity(graph);
alpar@896: 
alpar@896:   capacity.set(s_v1, 2);
alpar@896:   capacity.set(v1_v2, 2);
alpar@896:   capacity.set(s_v3, 1);
alpar@896:   capacity.set(v2_v4, 1);
alpar@896:   capacity.set(v2_v5, 1);
alpar@896:   capacity.set(v3_v5, 1);
alpar@896:   capacity.set(v4_t, 1);
alpar@896:   capacity.set(v5_t, 2);
alpar@896: 
alpar@896:   //  ConstMap<Edge, int> const1map(1);
alpar@896:   std::cout << "Mincostflows algorithm test..." << std::endl;
alpar@896: 
alpar@896:   MinCostFlows< ListGraph, ListGraph::EdgeMap<int>, ListGraph::EdgeMap<int> >
alpar@896:     surb_test(graph, length, capacity);
alpar@896: 
alpar@896:   int k=1;
alpar@896: 
alpar@896:   check(  surb_test.run(s,t,k) == 1 && surb_test.totalLength() == 19,"One path, total length should be 19");
alpar@896: 
alpar@896:   check(surb_test.checkComplementarySlackness(), "Is the primal-dual solution pair really optimal?");
alpar@896:   
alpar@896:   k=2;
alpar@896:   
alpar@896:   check(  surb_test.run(s,t,k) == 2 && surb_test.totalLength() == 41,"Two paths, total length should be 41");
alpar@896: 
alpar@896:   check(surb_test.checkComplementarySlackness(), "Is the primal-dual solution pair really optimal?");
alpar@896:   
alpar@896:   
alpar@896:   k=4;
alpar@896: 
alpar@896:   check(  surb_test.run(s,t,k) == 3 && surb_test.totalLength() == 64,"Three paths, total length should be 64");
alpar@896: 
alpar@896:   check(surb_test.checkComplementarySlackness(), "Is the primal-dual solution pair really optimal?");
alpar@896: 
alpar@896: 
alpar@896:   cout << (passed ? "All tests passed." : "Some of the tests failed!!!")
alpar@896:        << endl;
alpar@896: 
alpar@896:   return passed ? 0 : 1;
alpar@896: 
alpar@896: }