1

可能是一个新手问题，我需要你的帮助来比较两个图，顶点的数量和名称相同。

我的主题大纲是：

Graph origG, computedG;
...
...
int nmbrSameEdges, nmbrExteraEdges, nmbrMissingEdges, nmbrIncorrectEdges;

void compareGraph(origG, computedG, nmbrSameEdges, nmbrExteraEdges, nmbrMissingEdges, nmbrIncorrectEdges);

提前致谢。

4

1 回答 1

4

我想我会喜欢在这里做一个手指练习：

警告 据我所见，输出并不完全符合我的预期。我认为这是因为为边缘描述符定义的比较没有按我的预期工作。

如果您遇到任何问题，您可以使用同构算法

使用演示

int main()
{
    typedef adjacency_list<> Graph;

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

    diff<Graph> comparison;
    compare(a,b, comparison);

    std::cout << comparison;

/*
 *    std::cout << "edge_stats.same:      " << comparison.edge_stats.same.size()      << std::endl;
 *    std::cout << "edge_stats.extra:     " << comparison.edge_stats.extra.size()     << std::endl;
 *    std::cout << "edge_stats.missing:   " << comparison.edge_stats.missing.size()   << std::endl;
 *    std::cout << "edge_stats.reversed:  " << comparison.edge_stats.reversed.size()  << std::endl;
 *    std::cout << "vertex_stats.same:    " << comparison.vertex_stats.same.size()    << std::endl;
 *    std::cout << "vertex_stats.extra:   " << comparison.vertex_stats.extra.size()   << std::endl;
 *    std::cout << "vertex_stats.missing: " << comparison.vertex_stats.missing.size() << std::endl;
 */

}

输出（使用 Boost Spirit）：

Diff stats: 
--------------
Vertices
    Same:    3 (0, 1, 2)
    Extra:   1 (3)
    Missing: 0 ()
Edges
    Same:    0
    Extra:   1
    Missing: 1
    Reversed:1

请注意，示例代码被注释了，它显示了如何在不使用 Boost Spirit 的情况下打印差异结果。

执行

#include <boost/graph/adjacency_list.hpp>
#include <boost/graph/reverse_graph.hpp>
#include <boost/graph/iteration_macros.hpp>
#include <boost/spirit/include/karma.hpp>
using namespace boost;

namespace setops
{
    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;
    }
}

template <typename descriptor>
    struct stats 
    { 
        std::set<descriptor> same, extra, missing; 
    };
template <typename descriptor>
    struct directed_stats : stats<descriptor>
    { 
        std::set<descriptor> reversed; 
    };

template <typename G>
    struct diff
    {
        typedef typename graph_traits<G>::vertex_descriptor vdesc;
        typedef typename graph_traits<G>::edge_descriptor   edesc;

        stats<vdesc> vertex_stats;
        directed_stats<edesc> edge_stats;

        template <typename G2>
        friend std::ostream& operator<<(std::ostream& os, const diff<G2>& d)
        {
            using namespace boost::spirit::karma;

            os << "Diff stats: \n"
                  "--------------\n"
                  "Vertices\n"
               << format("\tSame:    " << int_ << " (" << -(auto_ % ", ") << ")" << eol, d.vertex_stats.same.size(), d.vertex_stats.same)
               << format("\tExtra:   " << int_ << " (" << -(auto_ % ", ") << ")" << eol, d.vertex_stats.extra.size(), d.vertex_stats.extra)
               << format("\tMissing: " << int_ << " (" << -(auto_ % ", ") << ")" << eol, d.vertex_stats.missing.size(), d.vertex_stats.missing);
            os << "Edges\n"
               << format("\tSame:    " << int_ << eol, d.edge_stats.same.size())
               << format("\tExtra:   " << int_ << eol, d.edge_stats.extra.size())
               << format("\tMissing: " << int_ << eol, d.edge_stats.missing.size())
               << format("\tReversed:" << int_ << eol, d.edge_stats.missing.size());
            return os;
        }
    };

template <typename G>
void compare(const G& a, const G& b, diff<G>& result)
{
    std::set<typename diff<G>::vdesc> av, bv;
    std::set<typename diff<G>::edesc> ae, be, re;

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

    // reverse adapt b:
    reverse_graph<G> r(b);
    BGL_FORALL_EDGES_T(e, r, G) re.insert(e);

    using namespace setops;
    result.edge_stats.same    = (ae / be);
    result.edge_stats.extra   = (be - ae);
    result.edge_stats.missing = (ae - be);
    result.edge_stats.reversed= (ae / re);

    BGL_FORALL_VERTICES_T(v, a, G) av.insert(v);
    BGL_FORALL_VERTICES_T(v, b, G) bv.insert(v);
    result.vertex_stats.same    = (av / bv);
    result.vertex_stats.extra   = (bv - av);
    result.vertex_stats.missing = (av - bv);
}
于 2011-11-18T00:51:02.877 回答