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 。

解法：

水水更健康系列，直接枚举起点，spfa跑一下

//BZOJ 1295 spfa

#include <bits/stdc++.h>
using namespace std;

const int inf = 0x3f3f3f3f;
int mp[40][40];
char s[40];
int n, m, t;
double ans;
int inq[40][40], dis[40][40];
int dir[4][2]={{0,1},{0,-1},{1,0},{-1,0}};
struct node{int x, y;};
void spfa(int x, int y){
    memset(inq, 0, sizeof(inq));
    node now;
    now.x = x, now.y = y;
    for(int i=1; i<=n; i++){
        for(int j=1; j<=m; j++){
            dis[i][j]=inf;
        }
    }
    dis[x][y]=mp[x][y];
    queue<node>que;
    que.push(now);
    inq[x][y]=1;
    while(!que.empty()){
        now = que.front();
        que.pop();
        inq[now.x][now.y]=0;
        for(int i=0; i<4; i++){
            int tx = now.x+dir[i][0];
            int ty = now.y+dir[i][1];
            if(tx<1||tx>n) continue;
            if(ty<1||ty>m) continue;
            if(dis[tx][ty]>dis[now.x][now.y]+mp[tx][ty]){
                if(dis[now.x][now.y]+mp[tx][ty]>t) continue;
                dis[tx][ty]=dis[now.x][now.y]+mp[tx][ty];
                if(!inq[tx][ty]){
                    que.push(node{tx,ty});
                    inq[tx][ty]=1;
                }
            }
        }
    }
    double tmp=0;
    for(int i=1; i<=n; i++){
        for(int j=1; j<=m; j++){
            if(dis[i][j]<=t){
                tmp = max(tmp, sqrt((x-i)*(x-i)+(y-j)*(y-j)));
            }
        }
    }
    ans = max(ans, tmp);
}

int main(){
    scanf("%d%d%d", &n,&m,&t);
    for(int i=1; i<=n; i++){
        scanf("%s", s+1);
        for(int j=1; j<=m; j++){
            mp[i][j]=s[j]-'0';
        }
    }
    for(int i=1; i<=n; i++){
        for(int j=1; j<=m; j++){
            spfa(i, j);
        }
    }
    printf("%.6f\n", ans);
    return 0;
}