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迭代向量（超级pedants ²会使用vector<int>::size_type）索引您当前感兴趣的向量元素。在读取P的最后一个元素时终止。

参考文献：

¹ Iterator invalidation rules

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

C ++中的Bron Kerbosch算法

1 个答案: