#include <lemon/circulation.h>

 #include <lemon/smart_graph.h>

 #include <lemon/lgf_reader.h>

 #include <lemon/dimacs.h>

 #include <lemon/time_measure.h>

 #include <lemon/error.h>

 const char test_name[] = "circulation";

 using namespace lemon;

 int main(int argc, char **argv)

 {

   if(argc!=2) exit(1);

   typedef int Value;

   SmartDigraph g;

   SmartDigraph::ArcMap<Value> lo_cap(g);

   SmartDigraph::ArcMap<Value> up_cap(g);

   SmartDigraph::ArcMap<Value> cost(g);

   SmartDigraph::NodeMap<Value> sup(g);

   Timer ti;

   try {

     std::ifstream input;

     input.open(argv[1]);

     readDimacsMin(input, g, lo_cap, up_cap, cost, sup);

   } catch (FormatError& error) {

     std::cerr << error.what() << std::endl;

     return 1;

   }

   std::cerr << "Read the file: " << ti << '\n';

   ti.restart();

   Circulation<SmartDigraph,

     SmartDigraph::ArcMap<Value>,SmartDigraph::ArcMap<Value>,

     SmartDigraph::NodeMap<Value> > circ(g,lo_cap,up_cap,sup);

   std::cerr << "Setup Circulation class: " << ti << '\n';

   ti.restart();

   bool res = circ.run();

   if(res)

     {

       std::cerr << "A feasible circulation is found: " << ti << "\n";

       ti.restart();

       bool res2 = circ.checkFlow();

       std::cerr << "Checked in time " << ti << "\n";

       if(res2)

         std::cerr << "Success!\nn";

       else

         std::cerr << "Oops!!!!\n\n";

     }

   else

     {

       std::cerr << "A dual solution is found: " << ti << "\n";

       ti.restart();

       bool res2 = circ.checkBarrier();

       std::cerr << "Checked in time " << ti << "\n";

       if(res2)

         std::cerr << "Success!\nn";

       else

         std::cerr << "Dual-Oops!!!!\n\n";

     }

 }