| 
athos@671
 | 
     1  | 
// -*- c++ -*-
  | 
| 
athos@671
 | 
     2  | 
#include <iostream>
  | 
| 
athos@671
 | 
     3  | 
#include <fstream>
  | 
| 
athos@671
 | 
     4  | 
  | 
| 
athos@671
 | 
     5  | 
#include <sage_graph.h>
  | 
| 
athos@671
 | 
     6  | 
//#include <smart_graph.h>
  | 
| 
athos@671
 | 
     7  | 
#include <hugo/dimacs.h>
  | 
| 
athos@671
 | 
     8  | 
#include <hugo/time_measure.h>
  | 
| 
athos@671
 | 
     9  | 
#include <hugo/for_each_macros.h>
  | 
| 
athos@671
 | 
    10  | 
#include <bfs_dfs.h>
  | 
| 
athos@671
 | 
    11  | 
  | 
| 
athos@671
 | 
    12  | 
using namespace hugo;
  | 
| 
athos@671
 | 
    13  | 
  | 
| 
athos@671
 | 
    14  | 
int main() {
 | 
| 
athos@671
 | 
    15  | 
  typedef SageGraph Graph; 
  | 
| 
athos@671
 | 
    16  | 
  typedef Graph::Node Node;
  | 
| 
athos@671
 | 
    17  | 
  typedef Graph::NodeIt NodeIt;
  | 
| 
athos@671
 | 
    18  | 
  typedef Graph::Edge Edge;
  | 
| 
athos@671
 | 
    19  | 
  typedef Graph::EdgeIt EdgeIt;
  | 
| 
athos@671
 | 
    20  | 
  typedef Graph::OutEdgeIt OutEdgeIt;
  | 
| 
athos@671
 | 
    21  | 
  | 
| 
athos@671
 | 
    22  | 
  Graph g;
  | 
| 
athos@671
 | 
    23  | 
  Node s, t;
  | 
| 
athos@671
 | 
    24  | 
  Graph::EdgeMap<int> cap(g);
  | 
| 
athos@671
 | 
    25  | 
  //readDimacsMaxFlow(std::cin, g, s, t, cap);
  | 
| 
athos@671
 | 
    26  | 
  readDimacs(std::cin, g);
  | 
| 
athos@671
 | 
    27  | 
  | 
| 
athos@671
 | 
    28  | 
  Graph::NodeMap<OutEdgeIt> pred(g);
  | 
| 
athos@671
 | 
    29  | 
  | 
| 
athos@671
 | 
    30  | 
  Timer ts;
  | 
| 
athos@671
 | 
    31  | 
  /*
  | 
| 
athos@671
 | 
    32  | 
  {
 | 
| 
athos@671
 | 
    33  | 
    ts.reset();
  | 
| 
athos@671
 | 
    34  | 
    Graph::NodeMap<bool> reached(g);
  | 
| 
athos@671
 | 
    35  | 
    reached.set(s, true);
  | 
| 
athos@671
 | 
    36  | 
    pred.set(s, INVALID);
  | 
| 
athos@671
 | 
    37  | 
    std::queue<Node> bfs_queue;
  | 
| 
athos@671
 | 
    38  | 
    bfs_queue.push(t);
  | 
| 
athos@671
 | 
    39  | 
    while (!bfs_queue.empty()) {
 | 
| 
athos@671
 | 
    40  | 
      Node v=bfs_queue.front();	
  | 
| 
athos@671
 | 
    41  | 
      bfs_queue.pop();
  | 
| 
athos@671
 | 
    42  | 
      OutEdgeIt e;
  | 
| 
athos@671
 | 
    43  | 
      for(g.first(e,v); g.valid(e); g.next(e)) {
 | 
| 
athos@671
 | 
    44  | 
	Node w=g.head(e);
  | 
| 
athos@671
 | 
    45  | 
	if (!reached[w]) {
 | 
| 
athos@671
 | 
    46  | 
	  bfs_queue.push(w);
  | 
| 
athos@671
 | 
    47  | 
	  reached.set(w, true);
  | 
| 
athos@671
 | 
    48  | 
	  pred.set(w, e);
  | 
| 
athos@671
 | 
    49  | 
	}
  | 
| 
athos@671
 | 
    50  | 
      }
  | 
| 
athos@671
 | 
    51  | 
    }
  | 
| 
athos@671
 | 
    52  | 
  | 
| 
athos@671
 | 
    53  | 
    std::cout << ts << std::endl;
  | 
| 
athos@671
 | 
    54  | 
  }
  | 
| 
athos@671
 | 
    55  | 
  */
  | 
| 
athos@671
 | 
    56  | 
  | 
| 
athos@671
 | 
    57  | 
  {
 | 
| 
athos@671
 | 
    58  | 
    ts.reset();
  | 
| 
athos@671
 | 
    59  | 
    Graph::NodeMap<bool> bfs_reached(g);
  | 
| 
athos@671
 | 
    60  | 
    Graph::NodeMap<Edge> bfs_pred(g); 
  | 
| 
athos@671
 | 
    61  | 
    Graph::NodeMap<int> bfs_dist(g);
  | 
| 
athos@671
 | 
    62  | 
      
  | 
| 
athos@671
 | 
    63  | 
    Bfs< Graph, Graph::NodeMap<bool>, 
  | 
| 
athos@671
 | 
    64  | 
      Graph::NodeMap<Edge>, Graph::NodeMap<int> > 
  | 
| 
athos@671
 | 
    65  | 
      bfs(g,bfs_reached, bfs_pred, bfs_dist );
  | 
| 
athos@671
 | 
    66  | 
    bfs.run(s);
  | 
| 
athos@671
 | 
    67  | 
    /*
  | 
| 
athos@671
 | 
    68  | 
    pred.set(s, INVALID);
  | 
| 
athos@671
 | 
    69  | 
    while (!bfs.finished()) { 
 | 
| 
athos@671
 | 
    70  | 
      ++bfs; 
  | 
| 
athos@671
 | 
    71  | 
      if (g.valid(bfs) && bfs.isBNodeNewlyReached()) 
  | 
| 
athos@671
 | 
    72  | 
	pred.set(bfs.bNode(), bfs);
  | 
| 
athos@671
 | 
    73  | 
    }
  | 
| 
athos@671
 | 
    74  | 
    */
  | 
| 
athos@671
 | 
    75  | 
    std::cout << ts << std::endl;
  | 
| 
athos@671
 | 
    76  | 
  }
  | 
| 
athos@671
 | 
    77  | 
  | 
| 
athos@671
 | 
    78  | 
  return 0;
  | 
| 
athos@671
 | 
    79  | 
}
  |