HDU-2224-双调TSP

  双调tsp的模板题，暑假时看没看懂，现在很好理解方程。

  f[i][j]表示dis[1...i]+dis[1...j]的最短路径长度(max(i,j)之前的点全部经过且仅经过一次)，f[i][j]=f[j][i],只需要计算一半即可，边界就是f[1][2]=dis(P[1],P[2]);

if(i<j-1) f[i][j]=f[i][j-1]+d[j-1][j];

if(i==j-1) f[i][j]=MIN{f[k][j-1]+d[k][j] | k<j-1 }

if(i==j) 两点相遇，只有在i==n时才能计算,f[n][n]=f[n-1][n]+d[n-1][n];

#include<bits/stdc++.h>
using namespace std;
struct Point
{
    double x,y;
}P[220];
double f[210][210];
double d[210][210];
double dis(Point A,Point B){return sqrt((A.x-B.x)*(A.x-B.x)+(A.y-B.y)*(A.y-B.y));}
int main()
{
    int n,i,j,k;
    while(cin>>n){
        for(i=1;i<=n;++i)
            cin>>P[i].x>>P[i].y;
        for(i=1;i<=n;++i)
            for(j=i+1;j<=n;++j)
        d[i][j]=d[j][i]=dis(P[i],P[j]);
        f[1][2]=d[1][2];
        for(j=3;j<=n;++j)
        {
            for(i=1;i<j-1;++i)
                f[i][j]=f[i][j-1]+d[j-1][j];
            f[j-1][j]=9999999999;
            for(k=1;k<j-1;++k)
                f[j-1][j]=min(f[j-1][j],f[k][j-1]+d[k][j]);
        }
        printf("%.2f
",f[n-1][n]+d[n-1][n]);
    }
    return 0;
}