Tour UVA - 1347

Tour UVA - 1347
紫书P269
给定平面上n个点的坐（按x递增的顺序给出，x各不相同，且均为正整数），设计一条路线，从最左边的点出发，走到最右边的点后再返回，路径中经过的点横坐标严格递增或递减（紫书第一个图不对），要求除了最左点和最右点外每个点恰好经过一次，且路径总长度（欧几里得距离）最短。

一个人向右走在向左走，抽象为两个人都从最左点出发，前往最右点，计算两个人路径和妙啊 。

dp[i][j]为一个人在第i个点，另一个在第j个点，且前max（i,j）个点都走过时，还要走的路径（本来是想用已经走过的路径来表示，但是不好转移状态）

易得dp[i][j]=dp[j][i]，所以我们令i≥j，这样可以只更新一半的空间.

转移的状态有：
1、走在前面的人走到下一个点，即为dp[i+1][j]
2、走在后面的人走到下一个点，即为dp[i][i+1]，也就是dp[i+1][i]

#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
const int maxn = 1010;
const double inf = 0x3f3f3f3f;
struct point
{
	int x, y;
}a[maxn];
int n;
//从左到右再到左
//抽象成两个人同时从左边出发，到达最右边
double dp[maxn][maxn];
//一个人在i点，另一个人在j点时,还要走的最短路程
//dp[i][j]=dp[j][i]，所以令i>=j
double dis[maxn][maxn];
inline double cal(int &p, int &q)
{
	int tempx = a[p].x - a[q].x;
	int tempy = a[p].y - a[q].y;
	return sqrt(tempx * tempx + tempy * tempy);
}
double dfs(int one, int two)
{
	if (one == n)
	{
		return dis[two][n];
	}
	if (dp[one][two] < inf)
		return dp[one][two];
	//两个转移状态
	//1、从dp[i][j]转移到dp[i+1][j]
	//2、从dp[i][j]转移到dp[i][i+1],也就是dp[i+1][i]
	double temp1 = dfs(one + 1, two) + dis[one][one + 1];
	double temp2 = dfs(one + 1, one) + dis[two][one + 1];
	return dp[one][two] = min(temp1, temp2);
}
int main()
{
	while (scanf("%d", &n) != EOF)
	{
		for (int i = 1; i <= n; i++)
		{
			scanf("%d%d", &a[i].x, &a[i].y);
		}
		mem(dis, 0);
		for (int i = 1; i <= n; i++)
		{
			for (int j = i + 1; j <= n; j++)
			{
				dis[i][j] = dis[j][i] = cal(i, j);
			}
		}
		for (int i = 1; i <= n; i++)
		{
			for (int j = 1; j <= i; j++)
			{
				dp[i][j] = 2 * inf;
			}
		}
		dp[n][n] = 0;
		printf("%.2f\n", dfs(1, 1));
	}
}