bzoj1295（spfa+省选所需的灵活思维）

最新推荐文章于 2018-11-01 16:24:00 发布

原创最新推荐文章于 2018-11-01 16:24:00 发布 · 794 阅读

0 ·

CC 4.0 BY-SA版权

最短路专栏收录该内容

10 篇文章

订阅专栏

本文介绍了一种利用SPFA算法解决特定障碍物环境下寻找最远可达路径的方法。通过构建图模型并进行最短路径计算，文章详细阐述了如何在给定地图上找到在移除一定数量障碍物后的最大可达距离。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

事实证明。学过的，只要敢想，做出来就只是时间问题，注意思维的灵活性，正常方向既然行不通（省选就不会是正常的），那么可以尝试一下不一样的思维方向

会有新的体会

关于这一题：

它说移掉T块障碍物以后的最大长度，我们可以这样想，如果两个块可以在隔T块障碍物就能到达的话，那么他们就是可以连起来的

由此，从每个点出发做一次spfa，算出每个点到这个点经过最少的障碍物的数量，然后统计时枚举所有点对，如果dis[i][j]<=T，那么他们就是合法的，

更新ans即可

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
const int qx[4]={0,0,1,-1};
const int qy[4]={1,-1,0,0}; 
int a[9009];
int n,m,t,tot,id[350][305];
int map[1000][4],dis[1200][1200];
double  ans=0;
bool b[10009];
bool pan(int x,int y)
{
	if (x<1||x>n||y<1||y>m) return false;
	return true;
}
int main()
{ 
	scanf("%d%d%d",&n,&m,&t);
	char c;
	tot=0;
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++)
		{
			c='\0';
			while (c!='0'&&c!='1') scanf("%c",&c);
			id[i][j]=++tot;
			if (c=='1') a[tot]=1;
		}
	//////////
	memset(map,0,sizeof(map)); 
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++)
			for (int l=0;l<4;l++)
			if (pan(i+qx[l],j+qy[l])) map[id[i][j]][l]=id[i+qx[l]][j+qy[l]];
	//////////////////////////////////////////////////////
	memset(dis,0x3f3f3f3f,sizeof(dis));
	for (int i=1;i<=tot;i++) dis[i][i]=a[i];
	queue<int> q;
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++)
		{
			while (!q.empty()) q.pop();
			int from=id[i][j];
			memset(b,false,sizeof(b));
			b[from]=true;
			q.push(from);
			while (!q.empty())
			{
				int u=q.front();
				q.pop();
				for(int i=0;i<4;i++) if (map[u][i])
				{
					int v=map[u][i];
					if (dis[from][u]+a[v]<dis[from][v]) 
					{
						dis[from][v]=dis[from][u]+a[v];
						if (!b[v])
						{
							b[v]=true;
							q.push(v);
						}
					}
				}
				b[u]=false;
			}
		}
	for (int i=1;i<tot;i++)
	for (int j=i+1;j<=tot;j++)
	if (dis[i][j]<=t)
	{
		int x=i/m+1,y=i%m;
		if (i%m==0) x--,y=m;
		int xx=j/m+1,yy=j%m;
		if (j%m==0) xx--,yy=m;
		ans=max(ans,sqrt((x-xx)*(x-xx)*1.0+(y-yy)*(y-yy)*1.0));
	}	
	printf("%.6f",ans);
	return 0;	
}