poj 2677 双调欧几里得 dp

这个题目的思路可以参考小紫书的269页，动态规划初步那一章的思路。
这里虽然和书中描述的不一样，但是大致的思路是一样的：
dp[i][j]表示两个人分别在i位置和j位置的最小花费，当前这一步有两个决策i选择或者j选择，所以状态转移方程如下：
$dp[i+1][j] = min(dp[i+1][j],dp[i][j]+dis(p[j],p[i+1]))$ 表示当前节点是选择了j节点
$dp[i+1][j]=min(dp[i+1][j],dp[i][j]+dis(p[i],p[i+1]));$ 表示当前节点选择了i节点

AC代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int MX = 500 + 5;
const int INF = 0x3f3f3f3f;
double dp[MX][MX];
struct point {
    double x,y;
}p[MX];
double dis(point a,point b){
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double solve(int n)
{
    dp[1][0]=dis(p[0],p[1]);
    for(int i=0;i<n;i++){
        for(int j=0;j<i;j++){
            dp[i+1][i]=min(dp[i+1][i],dp[i][j]+dis(p[j],p[i+1]));
            dp[i+1][j]=min(dp[i+1][j],dp[i][j]+dis(p[i],p[i+1]));
        }
    }
    double ans=INF;
    for(int i=0;i<n-1;i++)
        ans=min(ans,dp[n-1][i]+dis(p[i],p[n-1]));
    return ans;
}
int main(){
    int n;
    while(scanf("%d",&n)!=EOF){
        for(int i=0;i<n;i++)
            scanf("%lf%lf",&p[i].x,&p[i].y);
        memset(dp,127,sizeof(dp));
        printf("%.2lf\n",solve(n));
    }
    return 0;
}