RhodeCode

  if(report) std::cerr << "Run Dijkstra: " << t << '\n';

}

template<class Value>

void solve_max(ArgParser &ap, std::istream &is, std::ostream &,

               Value infty, DimacsDescriptor &desc)

{

  bool report = !ap.given("q");

  Digraph g;

  Node s,t;

  Digraph::ArcMap<Value> cap(g);

  Timer ti;

  ti.restart();

  readDimacsMax(is, g, cap, s, t, infty, desc);

  if(report) std::cerr << "Read the file: " << ti << '\n';

  ti.restart();

  Preflow<Digraph, Digraph::ArcMap<Value> > pre(g,cap,s,t);

  if(report) std::cerr << "Setup Preflow class: " << ti << '\n';

  ti.restart();

  pre.run();

  if(report) std::cerr << "Run Preflow: " << ti << '\n';

  if(report) std::cerr << "\nMax flow value: " << pre.flowValue() << '\n';  

}

void solve_min(ArgParser &ap, std::istream &is, std::ostream &,

               Value infty, DimacsDescriptor &desc)

{

  bool report = !ap.given("q");

  Digraph g;

  Digraph::ArcMap<Value> lower(g), cap(g), cost(g);

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

  Timer ti;

  ti.restart();

  readDimacsMin(is, g, lower, cap, cost, sup, infty, desc);

  ti.stop();

  Value sum_sup = 0;

  for (Digraph::NodeIt n(g); n != INVALID; ++n) {

    sum_sup += sup[n];

  }

  if (report) {

    std::cerr << "Sum of supply values: " << sum_sup << "\n";

    if (sum_sup <= 0)

      std::cerr << "GEQ supply contraints are used for NetworkSimplex\n\n";

    else

      std::cerr << "LEQ supply contraints are used for NetworkSimplex\n\n";

  }

  if (report) std::cerr << "Read the file: " << ti << '\n';

  ti.restart();

  NetworkSimplex<Digraph, Value> ns(g);

  ns.lowerMap(lower).upperMap(cap).costMap(cost).supplyMap(sup);

  if (sum_sup > 0) ns.supplyType(ns.LEQ);

  if (report) std::cerr << "Setup NetworkSimplex class: " << ti << '\n';

  ti.restart();

  bool res = ns.run();

  if (report) {

    std::cerr << "Run NetworkSimplex: " << ti << "\n\n";

    std::cerr << "Feasible flow: " << (res ? "found" : "not found") << '\n';

  }

}

void solve_mat(ArgParser &ap, std::istream &is, std::ostream &,

              DimacsDescriptor &desc)

{

  bool report = !ap.given("q");

  Graph g;

  Timer ti;

  ti.restart();

  readDimacsMat(is, g, desc);

  if(report) std::cerr << "Read the file: " << ti << '\n';

  ti.restart();

  MaxMatching<Graph> mat(g);

  if(report) std::cerr << "Setup MaxMatching class: " << ti << '\n';

  ti.restart();

  mat.run();

  if(report) std::cerr << "Run MaxMatching: " << ti << '\n';

  if(report) std::cerr << "\nCardinality of max matching: "

                       << mat.matchingSize() << '\n';  

}

void solve(ArgParser &ap, std::istream &is, std::ostream &os,

           DimacsDescriptor &desc)

{

  std::stringstream iss(static_cast<std::string>(ap["infcap"]));

  Value infty;

  iss >> infty;

  if(iss.fail())

    {

      std::cerr << "Cannot interpret '"

                << static_cast<std::string>(ap["infcap"]) << "' as infinite"

                << std::endl;

      exit(1);

    }

  switch(desc.type)

    {

    case DimacsDescriptor::MIN:

      break;

    case DimacsDescriptor::MAX:

      solve_max<Value>(ap,is,os,infty,desc);

      break;

    case DimacsDescriptor::SP:

      solve_sp<Value>(ap,is,os,desc);

      break;

    case DimacsDescriptor::MAT:

      solve_mat(ap,is,os,desc);

      break;

    default:

      break;

    }

}

int main(int argc, const char *argv[]) {

  typedef SmartDigraph Digraph;

  typedef Digraph::Arc Arc;

  std::string inputName;

  std::string outputName;

  ArgParser ap(argc, argv);

      std::cout << "Problem type: ";

      switch(desc.type)

        {

        case DimacsDescriptor::MIN:

          std::cout << "min";

          break;

        case DimacsDescriptor::MAX:

          std::cout << "max";

          break;

        case DimacsDescriptor::SP:

          std::cout << "sp";

        case DimacsDescriptor::MAT:

          std::cout << "mat";

          break;

        default:

          exit(1);

          break;

        }

      std::cout << "\nNum of nodes: " << desc.nodeNum;

      std::cout << "\nNum of arcs:  " << desc.edgeNum;

      std::cout << "\n\n";

    }

  if(ap.given("double"))

  else if(ap.given("ldouble"))

#ifdef LEMON_HAVE_LONG_LONG

  else if(ap.given("long"))

#endif

  return 0;

}