src/work/jacint/reverse_bfs.h
author jacint
Sun, 22 Feb 2004 12:17:16 +0000
changeset 113 cf7b01232d86
parent 78 ecc1171307be
permissions -rw-r--r--
*** empty log message ***
jacint@78
     1
/*
jacint@78
     2
reverse_bfs
jacint@78
     3
by jacint
jacint@78
     4
Performs a bfs on the out Edges. It does not count predecessors, 
jacint@78
     5
only the distances, but one can easily modify it to know the pred as well.
jacint@78
     6
jacint@78
     7
Constructor: 
jacint@78
     8
jacint@78
     9
reverse_bfs(Graph& G, NodeIt t)
jacint@78
    10
jacint@78
    11
jacint@78
    12
jacint@78
    13
Member functions:
jacint@78
    14
jacint@78
    15
void run(): runs a reverse bfs from t
jacint@78
    16
jacint@78
    17
  The following function should be used after run() was already run.
jacint@78
    18
jacint@78
    19
int dist(NodeIt v) : returns the distance from v to t. It is the number of nodes if t is not reachable from v.
jacint@78
    20
jacint@78
    21
*/
jacint@78
    22
#ifndef REVERSE_BFS_H
jacint@78
    23
#define REVERSE_BFS_H
jacint@78
    24
jacint@78
    25
#include <queue>
jacint@78
    26
#include <list_graph.hh>
jacint@78
    27
jacint@78
    28
jacint@78
    29
namespace  marci {
jacint@78
    30
jacint@78
    31
  template <typename Graph>
jacint@78
    32
  class reverse_bfs {
jacint@78
    33
    typedef typename Graph::NodeIt NodeIt;
jacint@78
    34
    typedef typename Graph::EachNodeIt EachNodeIt;
jacint@78
    35
    typedef typename Graph::InEdgeIt InEdgeIt;
jacint@78
    36
jacint@78
    37
jacint@78
    38
    Graph& G;
jacint@78
    39
    NodeIt t;
jacint@78
    40
    typename Graph::NodeMap<int> distance;
jacint@78
    41
    
jacint@78
    42
jacint@78
    43
  public :
jacint@78
    44
jacint@78
    45
    /*
jacint@78
    46
      The distance of the Nodes is n, except t for which it is 0.
jacint@78
    47
    */
jacint@78
    48
    reverse_bfs(Graph& _G, NodeIt _t) : G(_G), t(_t), distance(G, G.nodeNum()) {
jacint@78
    49
      distance.set(t,0);
jacint@78
    50
    }
jacint@78
    51
    
jacint@78
    52
    void run() {
jacint@78
    53
jacint@78
    54
      typename Graph::NodeMap<bool> reached(G, false); 
jacint@78
    55
      reached.set(t, true);
jacint@78
    56
jacint@78
    57
      std::queue<NodeIt> bfs_queue;
jacint@78
    58
      bfs_queue.push(t);
jacint@78
    59
jacint@78
    60
      while (!bfs_queue.empty()) {
jacint@78
    61
jacint@78
    62
        NodeIt v=bfs_queue.front();	
jacint@78
    63
	bfs_queue.pop();
jacint@78
    64
jacint@78
    65
	for(InEdgeIt e=G.template first<InEdgeIt>(v); e.valid(); ++e) {
jacint@78
    66
	  NodeIt w=G.tail(e);
jacint@78
    67
	  if (!reached.get(w)) {
jacint@78
    68
	    bfs_queue.push(w);
jacint@78
    69
	    distance.set(w, distance.get(v)+1);
jacint@78
    70
	    reached.set(w, true);
jacint@78
    71
	  }
jacint@78
    72
	}
jacint@78
    73
      }
jacint@78
    74
    }
jacint@78
    75
jacint@78
    76
jacint@78
    77
jacint@78
    78
    int dist(NodeIt v) {
jacint@78
    79
      return distance.get(v);
jacint@78
    80
    }
jacint@78
    81
jacint@78
    82
jacint@78
    83
  };
jacint@78
    84
alpar@105
    85
} // namespace hugo
jacint@78
    86
jacint@78
    87
#endif //REVERSE_BFS_HH
jacint@78
    88
jacint@78
    89