UVa #1347 Tour (例题9-3)

来源：互联网发布：java定义一个方法show 编辑：程序博客网时间：2024/05/29 08:35

其实书里没有给出题目里全部的要求，导致一开始我感觉状态转移方程的设计不是很直观

原题中给出了几个条件/要求：

1. 点的横坐标distinct

2. Input按照点的横坐标由小到大的顺序给出

3. 从左往右走的过程中，不能“回头”，从右往左同理

其实这些条件给出之后，书里的状态转移方程就比较容易理解了：

1. 因为不能回头，所以从左往右时没走的点，回来的时候必须要走。

2. 从左往右和从右往左是对称的，可以等价为第二个人从左往右走了一遍。

3. 因为每个点都要走，所以从第二个点开始往后扫描，则每个点要么第一个人走，要么第二个人走。这样扫描到某一个程度后，若第一个人站在第i个点上，第二个人站在第j个点上，则max(i,j)之前的点肯定都已经走过了。据此可以设计出书中给出的状态和状态转移方程

4. d(i,j) = d(j,i)，因此可以设i > j，即任何时候都认为走在前面的人是第一个人。（其实第一个人还是第二个人不重要，他们走的路径随时可以对换，从而使第一个人走到横坐标比较大的点）。方便了递推法，另外减少了不必要的重复的运算。

刚开始学动态规划，还是觉得记忆化搜索做起来比较简单直接。不过递推法的效率要高一些，是学习的目标。

递推：Run Time: 0.022s

#define UVa  "LT9-3.1347.cpp"#include<cstdio>#include<vector>#include<cstring>#include<cmath>using namespace std;struct Pt {    int x, y;    Pt(){}    Pt(int a, int b):x(a),y(b){}};float dist(Pt&a, Pt&b) {    return sqrt((a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y));}//Global Variables.const int maxn = 1000 + 10;float dp[maxn][maxn];////int main() {    int n;    while(scanf("%d", &n) != EOF) {        vector<Pt> pts;        int a, b;        for(int i = 0; i < n; i ++) {            scanf("%d%d", &a, &b);            pts.push_back(Pt(a,b));        }        for(int i = 0; i < n - 2; i ++) {            int a = n - 1 - 1;            dp[a][i] = dist(pts[a], pts[n-1]) + dist(pts[i], pts[n-1]);        }        for(int i = n - 3; i >= 0; i --)            for(int j = 0; j < i; j ++)                dp[i][j] = min(dp[i + 1][j] + dist(pts[i],pts[i+1]), dp[i+1][i] + dist(pts[j], pts[i+1]));        printf("%.2f\n", dp[1][0] + dist(pts[1], pts[0]));    }    return 0;}

记忆化搜索：Run Time: 0.032s

#define UVa  "LT9-3.1347.cpp"#include<cstdio>#include<vector>#include<cstring>#include<cmath>using namespace std;struct Pt {    int x, y;    Pt(){}    Pt(int a, int b):x(a),y(b){}};float dist(Pt&a, Pt&b) {    return sqrt((a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y));}//Global Variables.const int maxn = 1000 + 10;float dp[maxn][maxn];vector<Pt> pts;////float solve(int i, int j) {    if(dp[i][j]) return dp[i][j];    dp[i][j] = min(solve(max(i,j) + 1, j) + dist(pts[i], pts[max(i,j)+1]), solve(i, max(i,j) + 1) + dist(pts[j], pts[max(i,j)+1]));    return dp[i][j];}int main() {    int n;    while(scanf("%d", &n) != EOF) {        pts.clear();        memset(dp, 0, sizeof(dp));        int a, b;        for(int i = 0; i < n; i ++) {            scanf("%d%d", &a, &b);            pts.push_back(Pt(a,b));        }        for(int i = 0; i < n - 2; i ++) {            int a = n - 1 - 1;            dp[a][i] = dist(pts[a], pts[n-1]) + dist(pts[i], pts[n-1]);            dp[i][a] = dist(pts[a], pts[n-1]) + dist(pts[i], pts[n-1]);        }        printf("%.2f\n", solve(0,0));    }    return 0;}

0 0