关键词:dp问题的转移
题意:求n个点的最短双调旅行商路径
解法:将环形问题分解为两个结尾点的问题,转换dp状态,易于dp转移!
1.dp[i][j]:1~i点与i相连的点是j(除了i-1外)的最短双调TSP路径
dp[i][j]=f[j+1][j]+sumd[i]-sumd[j+1]+DIS(i,j);
f[i][j]:以1开头,i、j为结尾,且从2,3…j-1作为中间节点的最短路径和
f[i][j]=min(f[i][j-1]+DIS(j,j-1),f[j][j-1]+DIS(i,j-1));
技巧:将dp问题转化为另一个dp问题,利于状态转移。
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<set>
#include<vector>
#include<map>
#define ll long long
using namespace std;
const int maxn = 3000+10;
int n;
double x[maxn],y[maxn];
double sumd[maxn],f[maxn][maxn],dp[maxn][maxn];
double ans;
double DIS(int i,int j){
return sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
}
int main(){
while(scanf("%d",&n)!=EOF){
memset(sumd,0,sizeof(sumd));
memset(f,0,sizeof(f));
memset(dp,0,sizeof(dp));
for(int i=1;i<=n;i++){
scanf("%lf%lf",&x[i],&y[i]);
if(i>=2){
sumd[i]=sumd[i-1]+DIS(i,i-1);
}
}
if(n==1) { printf("0.00\n"); continue; }
if(n==2) { printf("%.2lf\n",DIS(1,2)*2.0); continue; }
//f初始化
for(int i=1;i<=n;i++) f[i][1]=DIS(i,1);
for(int i=3;i<=n;i++){
for(int j=2;j<=i-1;j++){
f[i][j]=min(f[i][j-1]+DIS(j,j-1),f[j][j-1]+DIS(i,j-1));
}
}
for(int i=3;i<=n;i++){
for(int j=1;j<=i-2;j++){
dp[i][j]=f[j+1][j]+sumd[i]-sumd[j+1]+DIS(i,j);
}
}
for(int i=1;i<=n-2;i++){
if(i==1) ans=dp[n][1];
else ans=min(ans,dp[n][i]);
}
printf("%.2lf\n",ans);
}
return 0;
}