src/work/athos/min_cost_flow.cc
author athos
Tue, 25 May 2004 12:31:18 +0000
changeset 659 c5984e925384
child 661 d306e777117e
permissions -rw-r--r--
Almost ready.
     1 #include <iostream>
     2 #include "test_tools.h"
     3 #include <hugo/list_graph.h>
     4 #include <mincostflow.h>
     5 //#include <path.h>
     6 //#include <maps.h>
     7 
     8 using namespace std;
     9 using namespace hugo;
    10 
    11 
    12 
    13 bool passed = true;
    14 /*
    15 void check(bool rc, char *msg="") {
    16   passed = passed && rc;
    17   if(!rc) {
    18     std::cerr << "Test failed! ("<< msg << ")" << std::endl; \
    19  
    20 
    21   }
    22 }
    23 */
    24 
    25 
    26 int main()
    27 {
    28 
    29   typedef ListGraph::Node Node;
    30   typedef ListGraph::Edge Edge;
    31 
    32   ListGraph graph;
    33 
    34   //Ahuja könyv példája
    35 
    36   Node s=graph.addNode();
    37   Node v1=graph.addNode();  
    38   Node v2=graph.addNode();
    39   Node v3=graph.addNode();
    40   Node v4=graph.addNode();
    41   Node v5=graph.addNode();
    42   Node t=graph.addNode();
    43 
    44   Edge s_v1=graph.addEdge(s, v1);
    45   Edge v1_v2=graph.addEdge(v1, v2);
    46   Edge s_v3=graph.addEdge(s, v3);
    47   Edge v2_v4=graph.addEdge(v2, v4);
    48   Edge v2_v5=graph.addEdge(v2, v5);
    49   Edge v3_v5=graph.addEdge(v3, v5);
    50   Edge v4_t=graph.addEdge(v4, t);
    51   Edge v5_t=graph.addEdge(v5, t);
    52   
    53 
    54   ListGraph::EdgeMap<int> length(graph);
    55 
    56   length.set(s_v1, 6);
    57   length.set(v1_v2, 4);
    58   length.set(s_v3, 10);
    59   length.set(v2_v4, 5);
    60   length.set(v2_v5, 1);
    61   length.set(v3_v5, 4);
    62   length.set(v4_t, 8);
    63   length.set(v5_t, 8);
    64 
    65   /*
    66   ListGraph::EdgeMap<int> capacity(graph);
    67 
    68   capacity.set(s_v1, 2);
    69   capacity.set(v1_v2, 2);
    70   capacity.set(s_v3, 1);
    71   capacity.set(v2_v4, 1);
    72   capacity.set(v2_v5, 1);
    73   capacity.set(v3_v5, 1);
    74   capacity.set(v4_t, 1);
    75   capacity.set(v5_t, 2);
    76   */
    77 
    78   //  ConstMap<Edge, int> const1map(1);
    79   std::cout << "Enhanced capacity scaling algorithm test (for the mincostflow problem)..." << std::endl;
    80 
    81   MinCostFlow< ListGraph, ListGraph::EdgeMap<int>, ListGraph::NodeMap<int> >
    82     min_cost_flow_test(graph, length, supply_demand);
    83 
    84   int k=1;
    85 
    86   /*
    87   check(  min_cost_flow_test.run(s,t,k) == 1 && min_cost_flow_test.totalLength() == 19,"One path, total length should be 19");
    88 
    89   check(min_cost_flow_test.checkComplementarySlackness(), "Is the primal-dual solution pair really optimal?");
    90   
    91   k=2;
    92   
    93   check(  min_cost_flow_test.run(s,t,k) == 2 && min_cost_flow_test.totalLength() == 41,"Two paths, total length should be 41");
    94 
    95   check(min_cost_flow_test.checkComplementarySlackness(), "Is the primal-dual solution pair really optimal?");
    96   
    97   
    98   k=4;
    99 
   100   check(  min_cost_flow_test.run(s,t,k) == 3 && min_cost_flow_test.totalLength() == 64,"Three paths, total length should be 64");
   101 
   102   check(min_cost_flow_test.checkComplementarySlackness(), "Is the primal-dual solution pair really optimal?");
   103 
   104 
   105   cout << (passed ? "All tests passed." : "Some of the tests failed!!!")
   106        << endl;
   107 
   108   return passed ? 0 : 1;
   109   */
   110 }