/*
reverse_bfs
by jacint
Performs a bfs on the out edges. It does not count predecessors, 
only the distances, but one can easily modify it to know the pred as well.

Constructor: 

reverse_bfs(graph_type& G, node_iterator t)



Member functions:

void run(): runs a reverse bfs from t

  The following function should be used after run() was already run.

int dist(node_iterator v) : returns the distance from v to t. It is the number of nodes if t is not reachable from v.

*/
#ifndef REVERSE_BFS_HH
#define REVERSE_BFS_HH

#include <queue>

//#include <marci_list_graph.hh>
//#include <marci_property_vector.hh>



namespace  marci {

  template <typename graph_type>
  class reverse_bfs {

    typedef typename graph_type::NodeIt NodeIt;
    typedef typename graph_type::EdgeIt EdgeIt;
    typedef typename graph_type::EachNodeIt EachNodeIt;
    typedef typename graph_type::EachEdgeIt EachEdgeIt;
    typedef typename graph_type::OutEdgeIt OutEdgeIt;
    typedef typename graph_type::InEdgeIt InEdgeIt;
    typedef typename graph_type::SymEdgeIt SymEdgeIt;



    graph_type& G;
    NodeIt t;
//    node_property_vector<graph_type, edge_iterator> pred;
    //node_property_vector<graph_type, int>
    typename graph_type::NodeMap<int> distance;
    

  public :

    /*
      The distance of the nodes is n, except t for which it is 0.
    */
    reverse_bfs(graph_type& _G, NodeIt _t) : 
      G(_G), t(_t), 
      distance(G, G.nodeNum()) {
      distance.set(t,0);
    }
    
    void run() {

      //node_property_vector<graph_type, bool> 
      typename graph_type::NodeMap<bool> reached(G, false); 
      reached.set(t, true);

      std::queue<NodeIt> bfs_queue;
      bfs_queue.push(t);

      while (!bfs_queue.empty()) {

        NodeIt v=bfs_queue.front();	
	bfs_queue.pop();

	for(InEdgeIt e=G.template first<InEdgeIt>(v); e.valid(); ++e) {
	  NodeIt w=G.tail(e);
	  if (!reached.get(w)) {
	    bfs_queue.push(w);
	    distance.set(w, distance.get(v)+1);
	    reached.set(w, true);
	  }
	}
      }
    }



    int dist(NodeIt v) {
      return distance.get(v);
    }


  };

} // namespace hugo

#endif //REVERSE_BFS_HH