[SCOI2009]最长距离

Description
windy有一块矩形土地，被分为 NM 块 11 的小格子。 有的格子含有障碍物。 如果从格子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 1
3 3 0
001
001
110

Sample Output 1
1.414214

Sample Output 2
3.605551

Sample Input 2
4 3 0
001
001
011
000

Sample Output 3
2.828427

Sample Input 3
3 3 1
001
001
001

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 。

---

N,M很小，我们可以直接跑N×M次SPFA，求出任意两点之间经过最少障碍数，然后\(N^2\)枚举，判断障碍数，直接统计答案

/*program from Wolfycz*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
#define sqr(x) ((x)*(x))
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline char gc(){
	static char buf[1000000],*p1=buf,*p2=buf;
	return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;
}
inline int frd(){
	int x=0,f=1;char ch=gc();
	for (;ch<'0'||ch>'9';ch=gc())	if (ch=='-')    f=-1;
	for (;ch>='0'&&ch<='9';ch=gc())	x=(x<<1)+(x<<3)+ch-'0';
	return x*f;
}
inline int read(){
	int x=0,f=1;char ch=getchar();
	for (;ch<'0'||ch>'9';ch=getchar())	if (ch=='-')	f=-1;
	for (;ch>='0'&&ch<='9';ch=getchar())	x=(x<<1)+(x<<3)+ch-'0';
	return x*f;
}
inline void print(int x){
	if (x<0)    putchar('-'),x=-x;
	if (x>9)	print(x/10);
	putchar(x%10+'0');
}
const int N=30;
const int dx[4]={-1,1,0,0};
const int dy[4]={0,0,-1,1};
int dist[N+10][N+10][N+10][N+10];
int n,m,T;
struct S1{
	int x,y;
	void insert(int _x,int _y){x=_x,y=_y;}
}h[N*N+10];
char map[N+10][N+10];
bool vis[N+10][N+10];
bool in_map(int x,int y){return x>0&&x<=n&&y>0&&y<=m;}
void SPFA(int x,int y,int dis[N+10][N+10]){
	memset(vis,0,sizeof(vis));
	h[1].insert(x,y),dis[x][y]=(map[x][y]=='1'),vis[x][y]=1;
	int head=0,tail=1;
	while (head!=tail){
		if (++head>N*N)	head=1;
		int Nx=h[head].x,Ny=h[head].y;
		for (int k=0;k<4;k++){
			int tx=Nx+dx[k],ty=Ny+dy[k];
			if (in_map(tx,ty)&&dis[tx][ty]>dis[Nx][Ny]+(map[tx][ty]=='1')){
				dis[tx][ty]=dis[Nx][Ny]+(map[tx][ty]=='1');
				if (!vis[tx][ty]){
					if (++tail>N*N)	tail=1;
					h[tail].insert(tx,ty);
					vis[tx][ty]=1;
				}
			}
		}
		vis[Nx][Ny]=0;
	}
}
int main(){
	n=read(),m=read(),T=read();
	memset(dist,63,sizeof(dist));
	for (int i=1;i<=n;i++)	scanf("%s",map[i]+1);
	for (int i=1;i<=n;i++)
		for (int j=1;j<=n;j++)
			SPFA(i,j,dist[i][j]);
	double Ans=0;
	for (int ox=1;ox<=n;ox++)
		for (int oy=1;oy<=m;oy++)
			for (int ex=1;ex<=n;ex++)
				for (int ey=1;ey<=m;ey++)
					if (dist[ox][oy][ex][ey]<=T)
						Ans=max(Ans,sqrt(sqr(ex-ox)+sqr(ey-oy)));
	printf("%.6lf\n",Ans);
	return 0;
}