【bzoj1925】【SCOI2009】【最长距离】【spfa】

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

题解：直接求最短路不好求,我们转化一下。

因为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);   
}