Chris@16
|
1 //=======================================================================
|
Chris@16
|
2 // Copyright 2002 Indiana University.
|
Chris@16
|
3 // Authors: Andrew Lumsdaine, Lie-Quan Lee, Jeremy G. Siek
|
Chris@16
|
4 //
|
Chris@16
|
5 // Distributed under the Boost Software License, Version 1.0. (See
|
Chris@16
|
6 // accompanying file LICENSE_1_0.txt or copy at
|
Chris@16
|
7 // http://www.boost.org/LICENSE_1_0.txt)
|
Chris@16
|
8 //=======================================================================
|
Chris@16
|
9
|
Chris@16
|
10 #ifndef BOOST_GRAPH_DAG_SHORTEST_PATHS_HPP
|
Chris@16
|
11 #define BOOST_GRAPH_DAG_SHORTEST_PATHS_HPP
|
Chris@16
|
12
|
Chris@16
|
13 #include <boost/graph/topological_sort.hpp>
|
Chris@16
|
14 #include <boost/graph/dijkstra_shortest_paths.hpp>
|
Chris@16
|
15
|
Chris@16
|
16 // single-source shortest paths for a Directed Acyclic Graph (DAG)
|
Chris@16
|
17
|
Chris@16
|
18 namespace boost {
|
Chris@16
|
19
|
Chris@16
|
20 // Initalize distances and call depth first search
|
Chris@16
|
21 template <class VertexListGraph, class DijkstraVisitor,
|
Chris@16
|
22 class DistanceMap, class WeightMap, class ColorMap,
|
Chris@16
|
23 class PredecessorMap,
|
Chris@16
|
24 class Compare, class Combine,
|
Chris@16
|
25 class DistInf, class DistZero>
|
Chris@16
|
26 inline void
|
Chris@16
|
27 dag_shortest_paths
|
Chris@16
|
28 (const VertexListGraph& g,
|
Chris@16
|
29 typename graph_traits<VertexListGraph>::vertex_descriptor s,
|
Chris@16
|
30 DistanceMap distance, WeightMap weight, ColorMap color,
|
Chris@16
|
31 PredecessorMap pred,
|
Chris@16
|
32 DijkstraVisitor vis, Compare compare, Combine combine,
|
Chris@16
|
33 DistInf inf, DistZero zero)
|
Chris@16
|
34 {
|
Chris@16
|
35 typedef typename graph_traits<VertexListGraph>::vertex_descriptor Vertex;
|
Chris@16
|
36 std::vector<Vertex> rev_topo_order;
|
Chris@16
|
37 rev_topo_order.reserve(num_vertices(g));
|
Chris@16
|
38
|
Chris@16
|
39 // Call 'depth_first_visit', not 'topological_sort', because we don't
|
Chris@16
|
40 // want to traverse the entire graph, only vertices reachable from 's',
|
Chris@16
|
41 // and 'topological_sort' will traverse everything. The logic below
|
Chris@16
|
42 // is the same as for 'topological_sort', only we call 'depth_first_visit'
|
Chris@16
|
43 // and 'topological_sort' calls 'depth_first_search'.
|
Chris@16
|
44 topo_sort_visitor<std::back_insert_iterator<std::vector<Vertex> > >
|
Chris@16
|
45 topo_visitor(std::back_inserter(rev_topo_order));
|
Chris@16
|
46 depth_first_visit(g, s, topo_visitor, color);
|
Chris@16
|
47
|
Chris@16
|
48 typename graph_traits<VertexListGraph>::vertex_iterator ui, ui_end;
|
Chris@16
|
49 for (boost::tie(ui, ui_end) = vertices(g); ui != ui_end; ++ui) {
|
Chris@16
|
50 put(distance, *ui, inf);
|
Chris@16
|
51 put(pred, *ui, *ui);
|
Chris@16
|
52 }
|
Chris@16
|
53
|
Chris@16
|
54 put(distance, s, zero);
|
Chris@16
|
55 vis.discover_vertex(s, g);
|
Chris@16
|
56 typename std::vector<Vertex>::reverse_iterator i;
|
Chris@16
|
57 for (i = rev_topo_order.rbegin(); i != rev_topo_order.rend(); ++i) {
|
Chris@16
|
58 Vertex u = *i;
|
Chris@16
|
59 vis.examine_vertex(u, g);
|
Chris@16
|
60 typename graph_traits<VertexListGraph>::out_edge_iterator e, e_end;
|
Chris@16
|
61 for (boost::tie(e, e_end) = out_edges(u, g); e != e_end; ++e) {
|
Chris@16
|
62 vis.discover_vertex(target(*e, g), g);
|
Chris@16
|
63 bool decreased = relax(*e, g, weight, pred, distance,
|
Chris@16
|
64 combine, compare);
|
Chris@16
|
65 if (decreased)
|
Chris@16
|
66 vis.edge_relaxed(*e, g);
|
Chris@16
|
67 else
|
Chris@16
|
68 vis.edge_not_relaxed(*e, g);
|
Chris@16
|
69 }
|
Chris@16
|
70 vis.finish_vertex(u, g);
|
Chris@16
|
71 }
|
Chris@16
|
72 }
|
Chris@16
|
73
|
Chris@16
|
74 namespace detail {
|
Chris@16
|
75
|
Chris@16
|
76 // Defaults are the same as Dijkstra's algorithm
|
Chris@16
|
77
|
Chris@16
|
78 // Handle Distance Compare, Combine, Inf and Zero defaults
|
Chris@16
|
79 template <class VertexListGraph, class DijkstraVisitor,
|
Chris@16
|
80 class DistanceMap, class WeightMap, class ColorMap,
|
Chris@16
|
81 class IndexMap, class Params>
|
Chris@16
|
82 inline void
|
Chris@16
|
83 dag_sp_dispatch2
|
Chris@16
|
84 (const VertexListGraph& g,
|
Chris@16
|
85 typename graph_traits<VertexListGraph>::vertex_descriptor s,
|
Chris@16
|
86 DistanceMap distance, WeightMap weight, ColorMap color, IndexMap /*id*/,
|
Chris@16
|
87 DijkstraVisitor vis, const Params& params)
|
Chris@16
|
88 {
|
Chris@16
|
89 typedef typename property_traits<DistanceMap>::value_type D;
|
Chris@16
|
90 dummy_property_map p_map;
|
Chris@16
|
91 D inf =
|
Chris@16
|
92 choose_param(get_param(params, distance_inf_t()),
|
Chris@16
|
93 (std::numeric_limits<D>::max)());
|
Chris@16
|
94 dag_shortest_paths
|
Chris@16
|
95 (g, s, distance, weight, color,
|
Chris@16
|
96 choose_param(get_param(params, vertex_predecessor), p_map),
|
Chris@16
|
97 vis,
|
Chris@16
|
98 choose_param(get_param(params, distance_compare_t()), std::less<D>()),
|
Chris@16
|
99 choose_param(get_param(params, distance_combine_t()), closed_plus<D>(inf)),
|
Chris@16
|
100 inf,
|
Chris@16
|
101 choose_param(get_param(params, distance_zero_t()),
|
Chris@16
|
102 D()));
|
Chris@16
|
103 }
|
Chris@16
|
104
|
Chris@16
|
105 // Handle DistanceMap and ColorMap defaults
|
Chris@16
|
106 template <class VertexListGraph, class DijkstraVisitor,
|
Chris@16
|
107 class DistanceMap, class WeightMap, class ColorMap,
|
Chris@16
|
108 class IndexMap, class Params>
|
Chris@16
|
109 inline void
|
Chris@16
|
110 dag_sp_dispatch1
|
Chris@16
|
111 (const VertexListGraph& g,
|
Chris@16
|
112 typename graph_traits<VertexListGraph>::vertex_descriptor s,
|
Chris@16
|
113 DistanceMap distance, WeightMap weight, ColorMap color, IndexMap id,
|
Chris@16
|
114 DijkstraVisitor vis, const Params& params)
|
Chris@16
|
115 {
|
Chris@16
|
116 typedef typename property_traits<WeightMap>::value_type T;
|
Chris@16
|
117 typename std::vector<T>::size_type n;
|
Chris@16
|
118 n = is_default_param(distance) ? num_vertices(g) : 1;
|
Chris@16
|
119 std::vector<T> distance_map(n);
|
Chris@16
|
120 n = is_default_param(color) ? num_vertices(g) : 1;
|
Chris@16
|
121 std::vector<default_color_type> color_map(n);
|
Chris@16
|
122
|
Chris@16
|
123 dag_sp_dispatch2
|
Chris@16
|
124 (g, s,
|
Chris@16
|
125 choose_param(distance,
|
Chris@16
|
126 make_iterator_property_map(distance_map.begin(), id,
|
Chris@16
|
127 distance_map[0])),
|
Chris@16
|
128 weight,
|
Chris@16
|
129 choose_param(color,
|
Chris@16
|
130 make_iterator_property_map(color_map.begin(), id,
|
Chris@16
|
131 color_map[0])),
|
Chris@16
|
132 id, vis, params);
|
Chris@16
|
133 }
|
Chris@16
|
134
|
Chris@16
|
135 } // namespace detail
|
Chris@16
|
136
|
Chris@16
|
137 template <class VertexListGraph, class Param, class Tag, class Rest>
|
Chris@16
|
138 inline void
|
Chris@16
|
139 dag_shortest_paths
|
Chris@16
|
140 (const VertexListGraph& g,
|
Chris@16
|
141 typename graph_traits<VertexListGraph>::vertex_descriptor s,
|
Chris@16
|
142 const bgl_named_params<Param,Tag,Rest>& params)
|
Chris@16
|
143 {
|
Chris@16
|
144 // assert that the graph is directed...
|
Chris@16
|
145 null_visitor null_vis;
|
Chris@16
|
146 detail::dag_sp_dispatch1
|
Chris@16
|
147 (g, s,
|
Chris@16
|
148 get_param(params, vertex_distance),
|
Chris@16
|
149 choose_const_pmap(get_param(params, edge_weight), g, edge_weight),
|
Chris@16
|
150 get_param(params, vertex_color),
|
Chris@16
|
151 choose_const_pmap(get_param(params, vertex_index), g, vertex_index),
|
Chris@16
|
152 choose_param(get_param(params, graph_visitor),
|
Chris@16
|
153 make_dijkstra_visitor(null_vis)),
|
Chris@16
|
154 params);
|
Chris@16
|
155 }
|
Chris@16
|
156
|
Chris@16
|
157 } // namespace boost
|
Chris@16
|
158
|
Chris@16
|
159 #endif // BOOST_GRAPH_DAG_SHORTEST_PATHS_HPP
|