旅行商简化版 DP

传送门：http://www.cqoi.net:2012/problem.php?id=1508

题目描述：

旅行商简化版欧几里德旅行商(Euclidean Traveling Salesman)问题也就是货郎担问题一直是困扰全世界数学家、计算机学家的著名问题。现有的算法都没有办法在确定型机器上在多项式时间内求出最优解，但是有办法在多项式时间内求出一个较优解。 为了简化问题，而且保证能在多项式时间内求出最优解，J.L.Bentley提出了一种叫做bitonic tour的哈密尔顿环游。它的要求是任意两点(a,b)之间的相互到达的代价dist(a,b)=dist(b,a)且任意两点之间可以相互到达，并且环游的路线只能是从最西端单向到最东端，再单向返回最西端，并且是一个哈密尔顿回路。现在笛卡尔平面上有n(n<=1000)个点，每个点的坐标为(x,y)(-2^31<=x<=2^31,-2^31<=y<=2^31).

输入：

第一行一个整数n 。

接下来n行，每行两个整数x,y，表示某个点的坐标。

 输入中保证没有重复的两点，且保证没有两点的x坐标相同，并且保证最西端和最东端都只有一个点。

输出：

一行，即最短回路的长度，保留2位小数。

样例输入：

7
0 6
1 0
2 3
5 4
6 1
7 5
8 2

样例输出：

25.58

分析：

一道DP。

我们可以将一个人走一个经过所有点的回路的最小路径长度转化为两个人走不重复的路并经过所有的点的最小路径长度。

状态转移方程为f[i][j]=min(f[i-1][j]+getdis(i-1,i),f[i-1][k]+getdis(k,i))。

其中f为当前i个点都被经过一次，且靠后的人在j这个位置时所经过的路径长度。

看起来这种做法好像会漏掉很多情况，但其实并不会。某一靠前的人在点i1直接走到i2经过的路程比较短，在i1时可能并没有更新，但是它会在f[i2-1][k]这一步时更新。

code:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define MAXN 1000
#define LL long long
#define INF 1e18
using namespace std;
inline double min(double a,double b){return a<b?a:b;}
pair<LL,LL>node[MAXN+5];
int n;
double getdis(int i,int j){return sqrt((node[i].first-node[j].first)*(node[i].first-node[j].first)+(node[i].second-node[j].second)*(node[i].second-node[j].second));}
double f[MAXN+5][MAXN+5],ans=INF;
void dp()
{
	int i,j;
	for(i=3;i<=n;i++)
		for(j=1;j<i;j++)
		{
			if(j==i-1)
				for(int k=1;k<i-1;k++)
					f[i][j]=min(f[i][j],f[i-1][k]+getdis(i,k));
			else
				f[i][j]=min(f[i][j],f[i-1][j]+getdis(i-1,i));
		}
}
int main()
{
	int i;
	scanf("%d",&n);
	for(i=1;i<=n;i++)
		scanf("%I64d%I64d",&node[i].first,&node[i].second);
	sort(node+1,node+n+1);
	memset(f,0x42,sizeof f);
	f[1][1]=0;
	f[2][1]=getdis(1,2);
	dp();
	for(i=1;i<n;i++)
		ans=min(ans,f[n][i]+getdis(n,i));
	printf("%.2lf\n",ans);
}