UVALive 3305Tour（双调DP）

题目大意：一个平面，有n个点，每个点的x坐标都是不一样的，输入时按照点的x坐标从小到大输入，有一个人从1号点出发，从左往右走，到达最右边的点，然后又从最右边的点回到1号点，问你这段路程的最小总路程是多少？

思路：双调DP。设d[ i ][ j ] 表示从点 i 到 1 再从 1 到 j 的最小距离，且规定，i  > j ，由于前面那个规定，那么点 i 必须上面去走，求 d[ i ][ j ] 时分两种情况：（1）如果它是直接从 i - 1 过来的，也就是 i - 1 在上面那条路，那么就是 d[ i ][ j ] = d[ i - 1][ j ] ，i- 1 < j ；（2）i - 1 也可以在下面那条路，这就要找了，i 这个点它可以直接从 i - 1 之前点连过来，下面那条路的端点是 i - 1 ，d[ i ][ i- 1] = min( d[ j ][ i - 1 ] )，1 < j < i - 1，由于 d[ j ][ i - 1 ] 没有，因为 j < i - 1 ，但是自己想想就知道，d[ j ][ i - 1 ] = d[ i - 1 ][ j ]。综上，状态转移方程为：d[ i ][ j ] = d[ i - 1][ j ] + dis( i - 1 , i ) ，j < i - 1，d[ i ][ i - 1] = min(d[ i - 1 ][ j ])，j < i - 1。最后的答案当然就是 min( d[ n ][ i ] ) ,1 <= i < n。

自己反正是没想出来，想到前 i 个围成圈的最小值，由于是闭合的，这样就搞不出来，又想不出第二维。。。 原来把闭合的搞成不闭合的就好了，厉害哈！

代码如下：

#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
using namespace std;

const double INF = 1e9 ;

typedef pair<double,double> pdd;

#define MP make_pair

const int MAXN = 1111;

vector <pdd> point;

double dis(int a,int b)
{
    return sqrt((point[a].first - point[b].first)*(point[a].first - point[b].first) + (point[a].second - point[b].second)*(point[a].second - point[b].second));
}

double d[MAXN][MAXN];

int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        point.clear();
        point.push_back(MP(0,0));
        for(int i = 1;i<=n;i++)
        {
            double a,b;
            scanf("%lf%lf",&a,&b);
            point.push_back(MP(a,b));
        }
        if(n == 1)
        {
            puts("0");
            continue;
        }
        d[2][1] = dis(2,1);
        for(int i = 3;i<=n;i++)
        {
            for(int j = 1;j < i-1;j++)
            {
                d[i][j] = d[i - 1][j] + dis(i-1,i);
            }
            d[i][i - 1] = INF;
            for(int j = 1;j < i - 1;j++)
            {
                d[i][i - 1] = min(d[i][i - 1],d[i - 1][j] + dis(j,i));
            }
        }
        double ans = INF;
        for(int i = 1;i<n;i++)
            ans = min(ans,d[n][i] + dis(n,i));
        printf("%.2f\n",ans);
    }
    return 0;
}