题目大意：给定一个棋盘，一些格子上有障碍物，可以移除T个障碍物，求移除后所有能互相到达的点对中的最大欧几里得距离
m,n<=30，900个点，我们可以枚举起始点，跑一遍SPFA，求出所有经过不超过T个障碍物可达的点，统计ans即可
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 40
using 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 );
}
BZOJ 1295 SCOI2009 最长距离 SPFA+暴力
原文地址：http://blog.csdn.net/popoqqq/article/details/40542871
