BZOJ 1295 SCOI2009 最长距离 SPFA+暴力

来源：互联网发布：房价数据库编辑：程序博客网时间：2024/05/04 06:07

题目大意：给定一个棋盘，一些格子上有障碍物，可以移除T个障碍物，求移除后所有能互相到达的点对中的最大欧几里得距离

m,n<=30，900个点，我们可以枚举起始点，跑一遍SPFA，求出所有经过不超过T个障碍物可达的点，统计ans即可

#include<cmath>#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>#define M 40using namespace std;typedef pair<int,int> abcd;const int dx[]={0,0,1,-1};const int dy[]={1,-1,0,0};int n,m,T;double ans;char s[M];int map[M][M];abcd q[65540];unsigned short r,h;bool v[M][M];int f[M][M];void SPFA(abcd x){int i,j;memset(f,0x3f,sizeof f);q[++r]=x;f[x.first][x.second]=map[x.first][x.second];while(r!=h){abcd temp=q[++h];v[temp.first][temp.second]=0;for(i=0;i<4;i++){int xx=temp.first+dx[i];int yy=temp.second+dy[i];if(xx<=0||yy<=0||xx>m||yy>n)continue;if( f[xx][yy]>f[temp.first][temp.second]+map[xx][yy] && f[temp.first][temp.second]+map[xx][yy]<=T ){f[xx][yy]=f[temp.first][temp.second]+map[xx][yy];if(!v[xx][yy])v[xx][yy]=1,q[++r]=make_pair(xx,yy);}}}for(i=1;i<=m;i++)for(j=1;j<=n;j++)if(f[i][j]<=T)ans=max(ans,sqrt( (i-x.first)*(i-x.first) + (j-x.second)*(j-x.second) ));}int main(){int i,j;cin>>m>>n>>T;for(i=1;i<=m;i++){scanf("%s",s+1);for(j=1;j<=n;j++)map[i][j]=s[j]-'0';}for(i=1;i<=m;i++)for(j=1;j<=n;j++)SPFA(make_pair(i,j));printf("%.6lf\n", ans );}

0 0