/*

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