Question

我无法识别为什么此算法没有返回TSP的最短路径。

vector<int> tsp(int n, vector< vector<float> >& cost)
{
  long nsub = 1 << n;
  vector< vector<float> > opt(nsub, vector<float>(n));

  for (long s = 1; s < nsub; s += 2)
    for (int i = 1; i < n; ++i) {
      vector<int> subset;
      for (int u = 0; u < n; ++u)
        if (s & (1 << u))
          subset.push_back(u);

      if (subset.size() == 2)
        opt[s][i] = cost[0][i];

      else if (subset.size() > 2) {
        float min_subpath = FLT_MAX;
        long t = s & ~(1 << i);
        for (vector<int>::iterator j = subset.begin(); j != subset.end(); ++j)
          if (*j != i && opt[t][*j] + cost[*j][i] < min_subpath)
            min_subpath = opt[t][*j] + cost[*j][i];
        opt[s][i] = min_subpath;
      }
    }

  vector<int> tour;
  tour.push_back(0);

  bool selected[n];
  fill(selected, selected + n, false);
  selected[0] = true;

  long s = nsub - 1;

  for (int i = 0; i < n - 1; ++i) {
    int j = tour.back();
    float min_subpath = FLT_MAX;
    int best_k;
    for (int k = 0; k < n; ++k)
      if (!selected[k] && opt[s][k] + cost[k][j] < min_subpath) {
        min_subpath = opt[s][k] + cost[k][j];
        best_k = k;
      }
    tour.push_back(best_k);
    selected[best_k] = true;
    s -= 1 << best_k;
  }
  tour.push_back(0);

  return tour;
}

例如，在仅5个点（图中的5个不同节点）的距离cost矩阵上，算法返回的路径不是最优的。任何帮助识别明显或小错误将不胜感激。或者任何关于出错的有用提示。

Answer 1

有一点看起来很奇怪的是，即使我不是子集s的一部分，主循环也能完成。

换句话说，opt [17] [8]将设置为cost [0] [8]。 opt [17] [8]表示在节点8处的状态，并且访问了节点0和4（因为5 = 2 ^ 0 + 2 ^ 4）。

这应该被标记为不可能，因为如果我们在节点8，我们当然必须访问节点8！

我建议通过改变来防止这些情况发生：

for (int i = 1; i < n; ++i) {
  vector<int> subset;

到

 for (int i = 1; i < n; ++i) {
  vector<int> subset;
  if ((s&(1<<i))==0) {
     opt[s][i]=FLT_MAX;
     continue;
  }

Answer 2

嵌套循环for(j=遍历subset中的所有节点，包括起始节点。这导致使用未初始化的值opt[t][0]，因此使用不正确的最佳路径长度计算。

最简单的解决方法是从subset排除起始节点：

for (int u = 1; u < n; ++u)
...
if (subset.size() == 1)
...
else if (subset.size() > 1)

TSP的动态方法

2 个答案: