C++ 中的 Bron Kerbosch 算法

Question

我一直在练习我的 C++ 算法知识，但被标准 BK 实现所困。该算法输出的派系太多，我似乎不明白为什么。我将图表示为邻接表：

vector< list<int> > adjacency_list;

我的 BK 函数如下所示：

void graph::BronKerbosch(vector<int> R, vector<int> P, vector<int> X){

  if (P.empty() && X.empty()){
    result_cliques.insert(R);
  }

  for (int node : P){

    vector<int> intersection = {}, intersectionX = {};

    //N(P)
    for (int nodeP : adjacency_list[node]){
      for (int node2 : P){
        if (nodeP == node2){
          intersection.push_back(nodeP);
        }   
      }

      //N(X)
      for (int node3 : X){
        if (nodeP == node3){
          intersectionX.push_back(nodeP);
        }
      }
    }

    R.push_back(node);
    BronKerbosch(R,intersection,intersectionX);
    P.erase(remove(P.begin(),P.end(),node),P.end());
    X.push_back(node);    

  }
}

我称之为：

void graph::run_BronKerbosch(){

  vector<int> R,P,X;

  for (int i=1; i < adjacency_list.size(); i++) {
    P.push_back(i);
  }

  BronKerbosch(R,P,X);

  cout << "................\nClassic: " << result_cliques.size() << endl;
  for (auto clique : result_cliques){
    cout << "(";
    for (int node : clique){
      cout << node <<" ";
    }    
    cout << ")\n";    
  } 

}

我正在尝试实现该算法的基本版本，但我似乎在这里遗漏了一个细节。问题出在：

for (int node : P){

我应该以某种方式在第一个循环中使用 P 的副本吗？ （我在相关问题中看到了这一点）

感谢您的帮助。

Answer 1

是的，您应该复制一份，除非您提前保留空间，以便保证不会重新分配¹。 （请注意，C++ foreach 实现归结为底层的一堆迭代器。）

如果我是你，我会改写为老式的 for 循环，使用 std::size_t 迭代向量（超级学究² 将使用 std::vector<int>::size_type）为您当前感兴趣的向量元素编制索引。在读取P 的最后一个元素时终止。

---

参考资料：

¹Iterator invalidation rules

²C++ for-loop - size_type vs. size_t