本文主要是介绍旅行商问题(travelling salesman problem, TSP) 解题报告，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！
旅行商问题是个熟知的问题。这次是因为coursera上面选的算法课而写代码实现。这里做个简单总结。
测试程序：
25
20833.3333 17100.0000
20900.0000 17066.6667
21300.0000 13016.6667
21600.0000 14150.0000
21600.0000 14966.6667
21600.0000 16500.0000
22183.3333 13133.3333
22583.3333 14300.0000
22683.3333 12716.6667
23616.6667 15866.6667
23700.0000 15933.3333
23883.3333 14533.3333
24166.6667 13250.0000
25149.1667 12365.8333
26133.3333 14500.0000
26150.0000 10550.0000
26283.3333 12766.6667
26433.3333 13433.3333
26550.0000 13850.0000
26733.3333 11683.3333
27026.1111 13051.9444
27096.1111 13415.8333
27153.6111 13203.3333
27166.6667 9833.3333
27233.3333 10450.0000
说明：这是25个城市的二维坐标，城市之间的距离为欧几里得距离。求从一个城市出发，每个城市只走一遍，再回到原城市的最短路程。
背景：http://en.wikipedia.org/wiki/Travelling_salesman_problem。 wiki百科对TSP的历史，精确解法，近似解法等都做了归纳。
根据题目要求，我要求的是一个精确解。维基百科说得不是很具体。可能大家都知道的一个精确解法是用动态规划，时间复杂度为O(n^22^n)。对于这样一个时间和空间都是指数复杂度的解法，遇到20以上节点就非常吃力了。我对其具体实现想得也不是很清楚，所以没有采用。后来使用的方法是branch and bound。其实就是搜索加剪枝。如果一个部分解的下界已经超过了当前的最好解，则该部分解的分支就不必进行了。可以看出，对一个具有n!解空间的问题，剪枝是否有效（当然，前提是正确）对算法的时间和空间消耗有很大的影响。
我个人觉得很有用的两个资料是：
1.http://www.academic.marist.edu/~jzbv/algorithms/Branch%20and%20Bound.htm
2.http://lcm.csa.iisc.ernet.in/dsa/node187.html
简单说明，部分解指的是只确定了部分城市的路径，还有一部分城市没有探索。部分解的下界估计值得的已确定的路径长度加上没确定的城市的路径最小的可能值。合起来就作为这个部分解的下界。即这个部分解所有的可能的完整解的最小值的下界。
两种链接1中的方法很简单，对于一个部分解，估计它的下界使用的方法是，当前已确定的路径长度加上没有确定路径的城市的最短邻接路径长度之和。这种方法实现起来很简单，能给人信心，但是好处也就到这里。对这个25个城市的TSP问题，我跑了24个小时都没跑出结果。。。链接中提供了Java代码实现，我的C++实现见后面。因为没有跑出结果，所以不能保证我的程序的正确性。
#include <iostream>
#include <cassert>
#include <climits>
#include <vector>
#include <queue>
#include <ctime>
using namespace std;const int N = 25;
int n;
double cities[N][2];
double dis[N][N];
double minedge[N];struct TSPNode{vector<int> path;double plen;double bnd;TSPNode(){plen = 0;bnd = -1;}double pathlen(){return plen;}int size(){return path.size();}int lastcity(){if(size() < 1){fprintf(stdout, "error: path is empty.\n");return -1;}return path[path.size() - 1];}void add(int city){path.push_back(city);}void copypath(const vector<int> &p){path = p;}bool contains(int city){for(int i = 0; i < path.size(); ++i){if(path[i] == city){return true;}}return false;}double bound(){if(bnd >= 0){return bnd;}bnd = plen;bool marked[N];memset(marked, false, sizeof(marked));for(int i = 0; i < path.size() - 1; ++i){marked[path[i]] = true;}for(int i = 0; i < n; ++i){if(!marked[i]){bnd += minedge[i];}}return bnd;}bool operator<(TSPNode& other){return bound() < other.bound();} 这篇关于旅行商问题(travelling salesman problem, TSP) 解题报告的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！
