kruskal-telephone.cpp 2.4 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566
  1. //=======================================================================
  2. // Copyright 2001 Jeremy G. Siek, Andrew Lumsdaine, Lie-Quan Lee,
  3. //
  4. // Distributed under the Boost Software License, Version 1.0. (See
  5. // accompanying file LICENSE_1_0.txt or copy at
  6. // http://www.boost.org/LICENSE_1_0.txt)
  7. //=======================================================================
  8. /*
  9. IMPORTANT!!!
  10. ~~~~~~~~~~~~
  11. This example uses interfaces that have been deprecated and removed from Boost.Grpah.
  12. Someone needs to update it, as it does NOT compile.
  13. */
  14. #include <boost/config.hpp>
  15. #include <iostream>
  16. #include <fstream>
  17. #include <boost/lexical_cast.hpp>
  18. #include <boost/graph/graphviz.hpp>
  19. #include <boost/graph/kruskal_min_spanning_tree.hpp>
  20. int
  21. main()
  22. {
  23. using namespace boost;
  24. GraphvizGraph g_dot;
  25. read_graphviz("figs/telephone-network.dot", g_dot);
  26. typedef adjacency_list < vecS, vecS, undirectedS, no_property,
  27. property < edge_weight_t, int > > Graph;
  28. Graph g(num_vertices(g_dot));
  29. property_map < GraphvizGraph, edge_attribute_t >::type
  30. edge_attr_map = get(edge_attribute, g_dot);
  31. graph_traits < GraphvizGraph >::edge_iterator ei, ei_end;
  32. for (boost::tie(ei, ei_end) = edges(g_dot); ei != ei_end; ++ei) {
  33. int weight = lexical_cast < int >(edge_attr_map[*ei]["label"]);
  34. property < edge_weight_t, int >edge_property(weight);
  35. add_edge(source(*ei, g_dot), target(*ei, g_dot), edge_property, g);
  36. }
  37. std::vector < graph_traits < Graph >::edge_descriptor > mst;
  38. typedef std::vector < graph_traits < Graph >::edge_descriptor >::size_type size_type;
  39. kruskal_minimum_spanning_tree(g, std::back_inserter(mst));
  40. property_map < Graph, edge_weight_t >::type weight = get(edge_weight, g);
  41. int total_weight = 0;
  42. for (size_type e = 0; e < mst.size(); ++e)
  43. total_weight += get(weight, mst[e]);
  44. std::cout << "total weight: " << total_weight << std::endl;
  45. typedef graph_traits < Graph >::vertex_descriptor Vertex;
  46. for (size_type i = 0; i < mst.size(); ++i) {
  47. Vertex u = source(mst[i], g), v = target(mst[i], g);
  48. edge_attr_map[edge(u, v, g_dot).first]["color"] = "black";
  49. }
  50. std::ofstream out("figs/telephone-mst-kruskal.dot");
  51. graph_property < GraphvizGraph, graph_edge_attribute_t >::type &
  52. graph_edge_attr_map = get_property(g_dot, graph_edge_attribute);
  53. graph_edge_attr_map["color"] = "gray";
  54. graph_edge_attr_map["style"] = "bold";
  55. write_graphviz(out, g_dot);
  56. return EXIT_SUCCESS;
  57. }