Blame view

3rdparty/boost_1_81_0/libs/graph/example/bfs-example.cpp 3.06 KB
977ed18d   Hu Chunming   提交三方库
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
  //=======================================================================
  // Copyright 2001 Jeremy G. Siek, Andrew Lumsdaine, Lie-Quan Lee,
  //
  // Distributed under the Boost Software License, Version 1.0. (See
  // accompanying file LICENSE_1_0.txt or copy at
  // http://www.boost.org/LICENSE_1_0.txt)
  //=======================================================================
  #include <boost/graph/adjacency_list.hpp>
  #include <boost/graph/breadth_first_search.hpp>
  #include <boost/pending/indirect_cmp.hpp>
  #include <boost/range/irange.hpp>
  
  #include <iostream>
  
  using namespace boost;
  template < typename TimeMap >
  class bfs_time_visitor : public default_bfs_visitor
  {
      typedef typename property_traits< TimeMap >::value_type T;
  
  public:
      bfs_time_visitor(TimeMap tmap, T& t) : m_timemap(tmap), m_time(t) {}
      template < typename Vertex, typename Graph >
      void discover_vertex(Vertex u, const Graph& g) const
      {
          put(m_timemap, u, m_time++);
      }
      TimeMap m_timemap;
      T& m_time;
  };
  
  int main()
  {
      using namespace boost;
      // Select the graph type we wish to use
      typedef adjacency_list< vecS, vecS, undirectedS > graph_t;
      // Set up the vertex IDs and names
      enum
      {
          r,
          s,
          t,
          u,
          v,
          w,
          x,
          y,
          N
      };
      const char* name = "rstuvwxy";
      // Specify the edges in the graph
      typedef std::pair< int, int > E;
      E edge_array[] = { E(r, s), E(r, v), E(s, w), E(w, r), E(w, t), E(w, x),
          E(x, t), E(t, u), E(x, y), E(u, y) };
      // Create the graph object
      const int n_edges = sizeof(edge_array) / sizeof(E);
  #if defined(BOOST_MSVC) && BOOST_MSVC <= 1300
      // VC++ has trouble with the edge iterator constructor
      graph_t g(N);
      for (std::size_t j = 0; j < n_edges; ++j)
          add_edge(edge_array[j].first, edge_array[j].second, g);
  #else
      typedef graph_traits< graph_t >::vertices_size_type v_size_t;
      graph_t g(edge_array, edge_array + n_edges, v_size_t(N));
  #endif
  
      // Typedefs
      typedef graph_traits< graph_t >::vertices_size_type Size;
  
      // a vector to hold the discover time property for each vertex
      std::vector< Size > dtime(num_vertices(g));
      typedef iterator_property_map< std::vector< Size >::iterator,
          property_map< graph_t, vertex_index_t >::const_type >
          dtime_pm_type;
      dtime_pm_type dtime_pm(dtime.begin(), get(vertex_index, g));
  
      Size time = 0;
      bfs_time_visitor< dtime_pm_type > vis(dtime_pm, time);
      breadth_first_search(g, vertex(s, g), visitor(vis));
  
      // Use std::sort to order the vertices by their discover time
      std::vector< graph_traits< graph_t >::vertices_size_type > discover_order(
          N);
      integer_range< int > range(0, N);
      std::copy(range.begin(), range.end(), discover_order.begin());
      std::sort(discover_order.begin(), discover_order.end(),
          indirect_cmp< dtime_pm_type, std::less< Size > >(dtime_pm));
  
      std::cout << "order of discovery: ";
      for (int i = 0; i < N; ++i)
          std::cout << name[discover_order[i]] << " ";
      std::cout << std::endl;
  
      return EXIT_SUCCESS;
  }