高密度大图中的最大加权团

时间:2014-04-25 17:20:57

标签: graph clique clique-problem

是否有任何软件或算法描述可以在图表中找到具有~17000个加权顶点和~75%密度的已知​​顶点数的最大团(大约)?我试图使用Cliquer,但它太慢了(让我几天得到结果)。

关于我的问题,以防万一 - 这是一个sceduling问题,我有18个时间段,每个时间段可以填充不同数量的替代品。每个变量代表一个插槽的一种替代方案。因此,一个插槽的所有替代方案都是互斥的,并且不同插槽的一些替代方案是不兼容的。如果两个替代方案兼容,那么它们之间存在边缘。重量反映了替代品的价值。

1 个答案:

答案 0 :(得分:2)

Boost Graph Library实际上似乎有一个clique搜索算法的实现,虽然我还没能找到它的文档。但是,您可以查看the implementation source code of the algorithmthis example,了解它是如何使用的。为了您的目的,您可以执行以下操作(此代码使用标记-std=c++11编译并使用Boost 1.55.0和g ++ 4.8.2):

#include <boost/graph/undirected_graph.hpp>
#include <boost/graph/bron_kerbosch_all_cliques.hpp>
#include <set>
#include <iostream>

// this is the visitor that will process each clique found by the algorithm
template <typename VertexWeight>
struct max_weighted_clique {
    typedef typename boost::property_traits<VertexWeight>::key_type Vertex;
    typedef typename boost::property_traits<VertexWeight>::value_type Weight;

    max_weighted_clique(const VertexWeight& weight_map, std::set<Vertex>& max_clique, Weight& max_weight)
        : weight_map(weight_map), max_weight(max_weight), max_clique(max_clique) {
        // max_weight = -inf
        max_weight = std::numeric_limits<Weight>::lowest();
    }

    // this is called each time a clique is found
    template <typename Clique, typename Graph>
    void clique(const Clique& c, const Graph& g) {
        // check the current clique value
        std::set<Vertex> current_clique;
        Weight current_weight = Weight(0);
        for(auto it = c.begin(); it != c.end(); ++it) {
            current_clique.insert(*it);
            current_weight += weight_map[*it];
        }
        // if it is a bigger clique, replace
        if (current_weight > max_weight) {
            max_weight = current_weight;
            max_clique = current_clique;
        }
    }

    const VertexWeight& weight_map;
    std::set<Vertex> &max_clique;
    Weight& max_weight;
};


// may replace with long, double...
typedef int Weight;

// this struct defines the properties of each vertex
struct VertexProperty {
    Weight weight;
    std::string name;
};

int main(int argc, char *argv[]) {
    // graph type
    typedef boost::undirected_graph<VertexProperty> Graph;
    // vertex descriptor type
    typedef boost::graph_traits<Graph>::vertex_descriptor Vertex;

    // create the graph
    Graph g;

    // add vertices
    Vertex v1 = boost::add_vertex(g);
    g[v1].weight = 5; g[v1].name = "v1";
    Vertex v2 = boost::add_vertex(g);
    g[v2].weight = 2; g[v2].name = "v2";
    Vertex v3 = boost::add_vertex(g);
    g[v3].weight = 6; g[v3].name = "v3";
    // add edges
    boost::add_edge(v1, v2, g);
    boost::add_edge(v1, v3, g);
    boost::add_edge(v2, v3, g);

    // instantiate the visitor
    auto vertex_weight = boost::get(&VertexProperty::weight, g);
    std::set<Vertex> max_clique;
    Weight max_weight;
    max_weighted_clique<decltype(vertex_weight)> visitor(vertex_weight, max_clique, max_weight);

    // use the Bron-Kerbosch algorithm to find all cliques
    boost::bron_kerbosch_all_cliques(g, visitor);

    // now max_clique holds a set with the max clique vertices and max_weight holds its weight

    auto vertex_name = boost::get(&VertexProperty::name, g);
    std::cout << "Max. clique vertices:" << std::endl;
    for (auto it = max_clique.begin(); it != max_clique.end(); ++it) {
        std::cout << "\t" << vertex_name[*it] << std::endl;
    }
    std::cout << "Max. clique weight = " << max_weight << std::endl;

    return 0;
}

我不知道这个代码是否会比Cliquer更好或更差(我在大图中没有做太多的搜索),但是你可以尝试一下(并分享你的结论:))< / p>

还有一个Parallel Boost Graph Library,但不幸的是它似乎没有实现一个clique搜索算法(甚至不是“隐藏”的算法)。

此外,我刚刚进入一个算法的并行实现,以找到名为MaxCliquePara的最大派系。我没有使用它,所以我不能谈论易用性或性能,但它看起来不错。但请注意,此实现搜索具有最大顶点数量的派系,这不是您正在寻找的 - 尽管您可以根据需要调整代码。

编辑:

似乎在quickbook sources of the library中存在关于BGL bron_kerbosch_all_cliques()的一些文档。虽然它是一个文档生成器,但它具有相当的可读性。就是这样。