[gmoj] [NOIP2013模拟11.5B组] [动态规划] 最短路径

[gmoj] [NOIP2013模拟11.5B组] [动态规划] 最短路径

题目

---

解题思路

考试时
用dfs暴力求出去时经过的点
取最优答案
妥妥的超时，20分

改正时
设f[i][j]为第一个人走到第i个点，第二个人走到第j个点，所走的最短路径
k为max(i,j)+1
可以理解为下一个要走的点，0-k的点已经走完
转移 f[k][j]=min(f[k][j],f[i][j]+dis(i,k));
f[i][k]=min(f[i][k],f[i][j]+dis(j,k));
f[0][0]=0
特判i==j的情况
注意b1和b2两个特殊点

---

代码

20分

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
struct lzf{
	int x,y;
}a[1020];
double da=0x7fffffffffff,f[1020][1020];
bool p[1020];
int n,b1,b2; 
void jl()  //求任意两点距离
{
	 for (int i=1;i<n;i++)
	     for (int j=i+1;j<=n;j++)
	         f[i][j]=f[j][i]=sqrt((a[i].x-a[j].x)*(a[i].x-a[j].x)+(a[i].y-a[j].y)*(a[i].y-a[j].y));
}
double dfsh()  //归程已定，直接累计 
{ 
	   double ans=0;
	   int dq=n;
	   for (int i=n-1;i>0;i--)
	       if (!p[i])
	          ans+=f[dq][i],dq=i;
	   return ans;
}
void dfsq(int fa,int dep,double sum)  //暴搜
{ 
	   if (dep==n)  //已经走完来程 
	   {
	      sum+=f[fa][dep];
	   	  da=min(da,sum+dfsh());  
	   	  return;
	   }
	   if (dep==b1)  //必须走
	   {
	   	  p[dep]=1;
	   	  dfsq(dep,dep+1,sum+f[fa][dep]);
	   	  p[dep]=0;
	   	  return;
	   }
	   if (dep==b2)  //不能走
	   {
	   	  dfsq(fa,dep+1,sum);
	   	  return;
	   }
	   p[dep]=1;
	   dfsq(dep,dep+1,sum+f[fa][dep]);
	   p[dep]=0;
	   dfsq(fa,dep+1,sum);
}
int main()
{
	freopen("path.in","r",stdin);
	freopen("path.out","w",stdout);
	scanf("%d%d%d",&n,&b1,&b2);
	b1++,b2++;
	for (int i=1;i<=n;i++)
		scanf("%d%d",&a[i].x,&a[i].y);
	jl();
	dfsq(1,2,0);
	printf("%0.2f\n",da);
	fclose(stdin);
	fclose(stdout);
	return 0;
}

100分

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
const int maxx=200001020;
struct lzf{
	int x,y;
}a[1020];
double da=0x7ffffff,f[1020][1020],d[1020][1020];
bool p[1020];
int n,b1,b2; 
void jl()
{
	 for (int i=0;i<n;i++)
	     for (int j=i+1;j<n;j++)
	         f[i][j]=f[j][i]=sqrt((a[i].x-a[j].x)*(a[i].x-a[j].x)+(a[i].y-a[j].y)*(a[i].y-a[j].y));
}
int main()
{
	freopen("path.in","r",stdin);
	freopen("path.out","w",stdout);
	scanf("%d%d%d",&n,&b1,&b2);
	memset(d,0x7f,sizeof(d));
	for (int i=0;i<n;i++)
		scanf("%d%d",&a[i].x,&a[i].y);
	jl();
	d[0][0]=0;
    for (int i=0;i<n;i++)
        for (int j=0;j<n;j++)
        { 
            if ((i!=0)&&(i==j)) continue;
            int k=max(i,j)+1;
            if (k==n) 
            {
            	if (j==n-1) d[n-1][n-1]=min(d[n-1][n-1],d[i][j]+f[i][n-1]);
            	if (i==n-1) d[n-1][n-1]=min(d[n-1][n-1],d[i][j]+f[j][n-1]); 
            	continue;
			}
			if (k!=b2) d[k][j]=min(d[k][j],d[i][j]+f[i][k]);
		    if (k!=b1) d[i][k]=min(d[i][k],d[i][j]+f[j][k]);
		}
	printf("%.2f\n",d[n-1][n-1]);
	fclose(stdin);
	fclose(stdout);
	return 0;
}