POJ_S1E02_Recursion&DP(3)

OUTLINE:

TOUR//@@@双调欧几里得商旅问题@@@

the triangle

help Jimmy

最佳加法表达式

POJ1037 A Decorative Fen ce（难，见后……假期再见ヾ(￣▽￣)Bye~Bye~）

OpenJudge 2755 神奇的口袋

POJ1014 Dividing

POJ2373 Dividing the Path

Basic idea:

tour:

1.【转】找状态转移方程——

首先将给出的点排序，关键字x，重新编号，从左至右1，2，3，…，n。

定义p[i][j]，表示结点i到结点j之间的距离。

定义d[i][j]，表示从i连到1，再从1连到j，（注意，i>j，且并没有相连。）

对于任意一个点i来说，有两种连接方法，一种是如图（a）所示，i与i-1相连，另一种呢是如图（b），i与i-1不相连。

根据双调旅程，我们知道结点n一定与n相连，那么，如果我们求的d[n][n-1]，只需将其加上p[n-1][n]就是最短双调闭合路线。

根据上图，很容易写出方程式：

dp[i][j]=dp[i-1][j]+dist[i][i-1];//dp[i][j],i to j across the 0th point

dp[i][i-1]=min(dp[i][i-1],dp[i-1][j]+dist[j][i]);

2. @@类似于迷宫BFS模拟标注@@——从起点双调延伸

传说中的双调DP问题。

要走完全部的点，且要来回的话，我们可以简化为有两个小人，都是从最右边出发的。

我们假设小人中有快的小人和慢的小人，如果慢的小人在一步中超过快的小人，则只是交换了角色。

dp[i][j]指的就是快的人走到i，慢的人走到j的最短距离。

假如现在快的小人走到i-1，慢的小人走到j，则此时就只有两种情况

1. 快的小人走到i，则有dp[i][j]=min(dp[i][j],dp[i-1][j]+dis(i-1,i));

2. 慢的小人走到i，则有dp[i][i-1]=min(dp[i][i-1],dp[i-1][j]+dis(i,j));

最后，必然是快的小人走到终点n-1，我们还要搜寻慢的小人在那个位置时总的路径最短。

#include <stdio.h>
#include <math.h>
#define MIN(a,b)	((a<b)?a:b)

struct point
{
	int x, y;
};
inline double dist(const point p1, const point p2){
	return sqrt(pow(p1.x - p2.x*1.0, 2) + pow(p1.y - p2.y*1.0, 2));
}
int main()
{
	const int max_n = 50;
	const double INF = 20000 * 2 * 50;

	int n = 0;
	point points[max_n + 10];
	while (scanf("%d", &n) != EOF){
		int i = 0, j = 0;
		for (i = 0; i < n; ++i)
			scanf("%d%d", &points[i].x, &points[i].y);

		double dp[max_n + 10][max_n + 10] = { 0 };
		for (i = 0; i < n; ++i)
		for (j = 0; j < n; ++j)	
				dp[i][j] = INF;//注意初始化时都为INF
		

		dp[0][0] = 0;
		for (i = 1; i < n; ++i)
		for (j = 0; j < i; ++j){
			dp[i][i - 1] = MIN(dp[i][i - 1], dp[i - 1][j] + dist(points[i], points[j]));
			dp[i][j] = MIN(dp[i][j], dp[i - 1][j] + dist(points[i], points[i - 1]));
		}
		dp[n - 1][n - 1] = dp[n - 1][n - 2] + dist(points[n - 1], points[n - 2]);

		printf("%.2f\n", dp[n - 1][n - 1]);

	}

	return 0;
}