是否有任何软件或算法描述可以在图表中找到具有~17000个加权顶点和~75%密度的已知顶点数的最大团(大约)?我试图使用Cliquer,但它太慢了(让我几天得到结果)。
关于我的问题,以防万一 - 这是一个sceduling问题,我有18个时间段,每个时间段可以填充不同数量的替代品。每个变量代表一个插槽的一种替代方案。因此,一个插槽的所有替代方案都是互斥的,并且不同插槽的一些替代方案是不兼容的。如果两个替代方案兼容,那么它们之间存在边缘。重量反映了替代品的价值。
答案 0 :(得分:2)
Boost Graph Library实际上似乎有一个clique搜索算法的实现,虽然我还没能找到它的文档。但是,您可以查看the implementation source code of the algorithm和this 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()
的一些文档。虽然它是一个文档生成器,但它具有相当的可读性。就是这样。