// -*- c++ -*-

 #include <iostream>

 #include <fstream>

 #include <sage_graph.h>

 #include <lemon/smart_graph.h>

 #include <lemon/dimacs.h>

 #include <lemon/time_measure.h>

 //#include <lemon/for_each_macros.h>

 #include <bfs_dfs.h>

 using namespace lemon;

 using std::cout;

 using std::endl;

 int main() {

   typedef SageGraph Graph; 

   typedef Graph::Node Node;

   typedef Graph::NodeIt NodeIt;

   typedef Graph::Edge Edge;

   typedef Graph::EdgeIt EdgeIt;

   typedef Graph::OutEdgeIt OutEdgeIt;

   Graph g;

   //Node s;

   //Graph::EdgeMap<int> cap(g);

   //readDimacsMaxFlow(std::cin, g, s, t, cap);

   readDimacs(std::cin, g);

   NodeIt s;

   g.first(s);

   cout << g.nodeNum() << endl;

   cout << g.edgeNum() << endl;

   Graph::NodeMap<Edge> pred(g);

   cout << "iteration time of bfs written by hand..." << endl;

   Timer ts;

   {

     ts.reset();

     Graph::NodeMap<bool> reached(g);

     reached.set(s, true);

     pred.set(s, INVALID);

     std::queue<Node> bfs_queue;

     bfs_queue.push(s);

     while (!bfs_queue.empty()) {

       Node v=bfs_queue.front();	

       bfs_queue.pop();

       OutEdgeIt e;

       for(g.first(e,v); g.valid(e); g.next(e)) {

 	Node w=g.head(e);

 	if (!reached[w]) {

 	  bfs_queue.push(w);

 	  reached.set(w, true);

 	  pred.set(w, e);

 	}

       }

     }

     std::cout << ts << std::endl;

   }

   cout << "iteration time with bfs iterator..." << endl;

   {

     ts.reset();      

     BfsIterator< Graph, Graph::NodeMap<bool> > bfs(g);

     bfs.pushAndSetReached(s);

     pred.set(s, INVALID);

     while (!bfs.finished()) { 

       ++bfs; 

       if (g.valid(bfs) && bfs.isBNodeNewlyReached()) 

 	pred.set(bfs.head(), Graph::Edge(bfs));

     }

     std::cout << ts << std::endl;

   }

   return 0;

 }