Boost Graph边缘比较不起作用

时间:2016-01-21 10:33:13

标签: c++ boost graph boost-graph

我正在尝试用图表进行一些统计(比较2个图表)。但是当我比较边缘时会遇到问题。

因此,我创建了两个带有边的图,然后是顶点和边操作的一些模板。对于顶点,它的工作效果非常好,但对于边缘,它无法正常工作。

#include <boost/graph/adjacency_list.hpp>
#include <boost/graph/reverse_graph.hpp>
#include <boost/graph/iteration_macros.hpp>
#include <boost/graph/graph_utility.hpp>
typedef boost::adjacency_list < boost::vecS, boost::vecS,
                            boost::undirectedS > Graph;

template <typename T> std::set<T> operator-(const std::set<T>& a,
                                            const std::set<T>& b)
{
    std::set<T> r;
    std::set_difference(
                        a.begin(), a.end(),
                        b.begin(), b.end(),
                        std::inserter(r, r.end()));

    return r;
}

template <typename T> std::set<T> operator/(const std::set<T>& a,
                                            const std::set<T>& b)
{
    std::set<T> r;
    std::set_intersection(
                          a.begin(), a.end(),
                          b.begin(), b.end(),
                          std::inserter(r, r.end()));

    return r;
}



void compare(const Graph& a, const Graph& b)
{
    std::set<Graph::vertex_descriptor > av, bv, samev, extrav, missingv;
    std::set<Graph::edge_descriptor> ae, be, re, samee, extrae, missinge;

    BGL_FORALL_VERTICES_T(v, a, Graph) av.insert(v);
    BGL_FORALL_VERTICES_T(v, b, Graph) bv.insert(v);
    samev    = (av / bv);
    extrav   = (bv - av);
    missingv = (av - bv);


    BGL_FORALL_EDGES_T(e, a, Graph) ae.insert(e);
    BGL_FORALL_EDGES_T(e, b, Graph) be.insert(e);

    samee    = (ae / be);
    extrae   = (be - ae);
    missinge = (ae - be);

    // TODO(silgon): reverse_graph
    // boost::reverse_graph<Graph> r(b);
    // BGL_FORALL_EDGES_T(e, r, Graph) re.insert(e);
    std::cout << "---- Vertices ----\n"
              << "same: " << samev.size() << std::endl
              << "extra: " << extrav.size() << std::endl
              << "missing: " << missingv.size() << std::endl;

    std::cout << "---- Edges ----\n"
              << "same: " << samee.size() << std::endl
              << "extra: " << extrae.size() << std::endl
              << "missing: " << missinge.size() << std::endl;
}

int main()
{
    Graph a;
    {
        boost::graph_traits<Graph>::vertex_descriptor v, u, y;
        u = boost::vertex(1, a);
        v = boost::vertex(2, a);
        y = boost::vertex(3, a);
        boost::add_edge(u, v, a);
    }
    Graph b;
    {
        boost::graph_traits<Graph>::vertex_descriptor v, u, y;
        u = vertex(1, b);
        v = vertex(2, b);
        y = vertex(3, b);
        boost::add_edge(u, v, b);
        boost::add_edge(y, v, b);
    }

    const char* name = "123456";
    std::cout << "---Graph1---" << "\n";
    boost::print_graph(a);
    std::cout << "Edges: ";
    boost::print_edges(a,name);
    std::cout << "---Graph2---" << "\n";
    boost::print_graph(b);
    std::cout << "Edges: ";
    boost::print_edges(b,name);

    compare(a,b);
}

关于程序的结果,如下:

---Graph1---
0 <--> 
1 <--> 2 
2 <--> 1 
Edges: (2,3) 
---Graph2---
0 <--> 
1 <--> 2 
2 <--> 1 3 
3 <--> 2 
Edges: (2,3) (4,3) 
---- Vertices ----
same: 3
extra: 1
missing: 0
---- Edges ----
same: 0
extra: 2
missing: 1

你可以看到其中一条边是相同的(2,3),但是在进行操作时它没有考虑到它,所以在结果中它说的相同:0,并且额外或缺少的结果不起作用

2 个答案:

答案 0 :(得分:3)

深入研究Boost的edge_descriptor实现(在boost/graph/detail/edge.hpp中,我们发现每个边缘描述符都存储源和目标vertex_descriptor,但也是一个指向edge属性的指针。这个指针对于vertex_descriptor s相同的两条边是不同的。

这意味着您需要定义自己的edge_descriptor比较器,用于std::setset_intersection / set_difference操作,例如像这样:

struct edge_cmp
{
  bool operator () (const Graph::edge_descriptor& a, const Graph::edge_descriptor& b) const
  {
    if (a.m_source < b.m_source) { return true; }
    if (a.m_source > b.m_source) { return false; }
    return (a.m_target < b.m_target);
  }
};

此类比较器对象必须传递给您构造的所有集合以及交集/差异调用。

I've forked your code on ideone并相应修改。输出:

---- Vertices ----
same: 3
extra: 1
missing: 0
---- Edges ----
same: 1
extra: 1
missing: 0

答案 1 :(得分:1)

作为我问题的临时解决方案,而不是:

std::set<Graph::edge_descriptor> ae
BGL_FORALL_EDGES_T(e, a, Graph) ae.insert(e)

我为边缘做了以下代码。

std::set<std::pair<unsigned long, unsigned long> > ae;
BGL_FORALL_EDGES_T(e, ga, Graph) gae.insert(std::make_pair(e.m_source, 
                                                           e.m_target));

问题是boost图的边缘有“m_eproperty”,我不知道为什么会出现这种情况,它包含0x125d0c0之类的值。然后我根据边缘的源和目标创建了一对,然后我以与上面相同的方式(使用模板)进行评估。