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
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
1.414214
【输出样例二】
3.605551
【输出样例三】
2.828427
题解:直接求最短路不好求,我们转化一下。
因为n和m都比较小。我们可以暴力枚举两个点
看一下这两个点能否通过移动小于等于t个障碍物联通
所以我们从每个点跑一遍spfa
预处理一下从这个点到其他点联通最少需要移动几个障碍物即可。
代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
int a[50][50],l[100001][4],n,m,dis[50][50][50][50],t,dx[5],dy[5];
double ans;
bool f[50][50];
char ch[1001];
double calc(double x,double y,double xx,double yy){
return sqrt((x-xx)*(x-xx)+(y-yy)*(y-yy));
}
void spfa(int x,int y){
int h(0),t(1),xx,yy;
memset(f,false,sizeof(f));
dis[x][y][x][y]=a[x][y];f[x][y]=true;l[t][1]=x;l[t][2]=y;
while(h<=t){
int p=l[++h][1],q=l[h][2];
f[p][q]=false;
for (int i=1;i<=4;i++){
xx=p+dx[i];yy=q+dy[i];
if (xx>=1&&xx<=n&&yy>=1&&yy<=m&&dis[x][y][xx][yy]>dis[x][y][p][q]+a[xx][yy]){
dis[x][y][xx][yy]=dis[x][y][p][q]+a[xx][yy];
if (!f[xx][yy]){
l[++t][1]=xx;l[t][2]=yy;
f[xx][yy]=true;
}
}
}
}
}
int main(){
scanf("%d%d%d",&n,&m,&t);
for (int i=1;i<=n;i++){
scanf("%*c%s",&ch);
for (int j=0;j<strlen(ch);j++) a[i][j+1]=ch[j]-'0';
}
dx[1]=-1;dx[2]=0;dx[3]=1;dx[4]=0;
dy[1]=0;dy[2]=-1;dy[3]=0;dy[4]=1;
memset(dis,127/3,sizeof(dis));
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
spfa(i,j);
for(int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
for (int k=1;k<=n;k++)
for (int p=1;p<=m;p++){
if (dis[i][j][k][p]<=t){
ans=max(ans,calc((double)i,(double)j,double(k),double(p)));
}
}
printf("%.6lf",ans);
}