我正在尝试用图表进行一些统计(比较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,并且额外或缺少的结果不起作用
答案 0 :(得分:3)
深入研究Boost的edge_descriptor
实现(在boost/graph/detail/edge.hpp
中,我们发现每个边缘描述符都存储源和目标vertex_descriptor
,但也是一个指向edge属性的指针。这个指针对于vertex_descriptor
s相同的两条边是不同的。
这意味着您需要定义自己的edge_descriptor
比较器,用于std::set
和set_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
之类的值。然后我根据边缘的源和目标创建了一对,然后我以与上面相同的方式(使用模板)进行评估。