#ifndef DIMACS_HH

 #define DIMACS_HH

 #include <iostream>

 #include <string>

 #include <vector>

 namespace hugo {

   template<typename Graph, typename CapacityMap>

   void readDimacsMaxFlow(std::istream& is, Graph &G, typename Graph::TrivNodeIt &s, typename Graph::TrivNodeIt &t, CapacityMap& capacity) {

     G.clear();

     int cap;

     char d;

     std::string problem;

     char c;

     int i, j;

     std::string str;

     int n, m; 

     std::vector<typename Graph::TrivNodeIt> nodes;

     while (is>>c) {

       switch (c) {

       case 'c': //comment

 	getline(is, str);

 	break;

       case 't': //type

 	getline(is, str);

 	break;

       case 'p': //problem definition

 	is >> problem >> n >> m;

 	getline(is, str);

 	nodes.resize(n+1);

 	for (int k=1; k<=n; ++k) nodes[k]=G.addNode();

 	break;

       case 'n': //node definition

 	if (problem=="sp") { //shortest path problem

 	  is >> i;

 	  getline(is, str);

 	  s=nodes[i];

 	}

 	if (problem=="max") { //max flow problem

 	  is >> i >> d;

 	  getline(is, str);

 	  if (d=='s') s=nodes[i];

 	  if (d=='t') t=nodes[i];

 	}

 	break;

       case 'a':

 	is >> i >> j >> cap;

 	getline(is, str);

 	typename Graph::TrivEdgeIt e=G.addEdge(nodes[i], nodes[j]);

 	capacity.update();

 	capacity.set(e, cap);

 	break;

       }

     }

   }

 } //namespace marci

 #endif //DIMACS_HH