就是一道最短路的裸题，直接跑spfa就行了。（spfa死了）
最后在答案处判断是否障碍物太多，然后就直接找最大值就行。
（数据特别水，我错误算法60）
题干：
Description windy有一块矩形土地，被分为 N*M 块 1*1 的小格子。 有的格子含有障碍物。 如果从格子A可以走到格子B，那么两个格子的距离就为两个格子中心的欧几里德距离。 如果从格子A不可以走到格子B，就没有距离。 如果格子X和格子Y有公共边，并且X和Y均不含有障碍物，就可以从X走到Y。 如果windy可以移走T块障碍物，求所有格子间的最大距离。 保证移走T块障碍物以后，至少有一个格子不含有障碍物。 Input 输入文件maxlength.in第一行包含三个整数，N M T。 接下来有N行，每行一个长度为M的字符串，'0'表示空格子，'1'表示该格子含有障碍物。 Output 输出文件maxlength.out包含一个浮点数，保留6位小数。 Sample Input 【输入样例一】 3 3 0 001 001 110 【输入样例二】 4 3 0 001 001 011 000 【输入样例三】 3 3 1 001 001 001 Sample Output 【输出样例一】 1.414214 【输出样例二】 3.605551 【输出样例三】 2.828427 HINT 20%的数据，满足 1 <= N,M <= 30 ； 0 <= T <= 0 。 40%的数据，满足 1 <= N,M <= 30 ； 0 <= T <= 2 。 100%的数据，满足 1 <= N,M <= 30 ； 0 <= T <= 30 。
代码：
#include<iostream> #include<cstdio> #include<cmath> #include<ctime> #include<queue> #include<algorithm> #include<cstring> using namespace std; #define duke(i,a,n) for(int i = a;i <= n;i++) #define lv(i,a,n) for(int i = a;i >= n;i--) #define clean(a) memset(a,0,sizeof(a)) const int INF = 1 << 30; typedef long long ll; typedef double db; template <class T> void read(T &x) { char c; bool op = 0; while(c = getchar(), c < '0' || c > '9') if(c == '-') op = 1; x = c - '0'; while(c = getchar(), c >= '0' && c <= '9') x = x * 10 + c - '0'; if(op) x = -x; } template <class T> void write(T x) { if(x < 0) putchar('-'), x = -x; if(x >= 10) write(x / 10); putchar('0' + x % 10); } int n,m,t; int dx[5] = {0,0,-1,1}; int dy[5] = {1,-1,0,0}; int dis[40][40],a[40][40],mp[40][40]; bool inq[40][40]; bool p[40][40]; char s[40]; db ans = 0; struct node { int x,y; }q[100010]; void getans(int x,int y) { duke(i,x,n) { duke(j,1,m) { if(dis[i][j] <= t && (y - j) * (y - j) + (x - i) * (x - i) > ans) ans = (y - j) * (y - j) + (x - i) * (x - i); } } } void spfa(int x,int y) { int nowx,nowy,t = 1,w = 2,nx,ny; q[1].x = x;q[1].y = y; clean(inq); memset(dis,127,sizeof(dis)); inq[x][y] = 1; dis[x][y] = mp[x][y]; while(t <= w) { nowx = q[t].x;nowy = q[t].y; t++; for(int i = 0;i < 4;i++) { nx = nowx + dx[i]; ny = nowy + dy[i]; if(nx > n || nx < x || ny > m || ny < 1) continue; if(dis[nowx][nowy] + mp[nx][ny] < dis[nx][ny]) { dis[nx][ny] = dis[nowx][nowy] + mp[nx][ny]; if(!inq[nx][ny]) { q[++w].x = nx; q[w].y = ny; inq[nx][ny] = 1; } } } inq[nowx][nowy] = 0; } getans(x,y); } int main() { read(n);read(m);read(t); duke(i,1,n) { scanf("%s",s); duke(j,0,m - 1) mp[i][j + 1] = s[j] - '0'; } duke(i,1,n) { duke(j,1,m) { spfa(i,j); } } printf("%.6lf",sqrt(ans)); return 0; } /* 3 3 0 001 001 110 */