[BZOJ1295][SCOI2009]最长距离 最短路+枚举

1295: [SCOI2009]最长距离

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 1683  Solved: 912
[Submit][Status][Discuss]

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<cstring>
 #include<cstdio>
 #include<cstdlib>
 #include<cmath>
 #include<algorithm>
 #include<queue>
 using namespace std;
 const int tx[]={,,,-};
 const int ty[]={,-,,};
 struct data {
     int x,y,val;
     bool operator <(const data &t)const {
         return val>t.val;
     }
 };
 int n,m,T;
 char s[][];
 bool vis[][];
 int dis[][];
 void dij(int stx,int sty) {
     priority_queue<data> q;
     memset(vis,,sizeof(vis));
     memset(dis,,sizeof(dis));
     dis[stx][sty]=s[stx][sty]==''?:;
     q.push((data){stx,sty,dis[stx][sty]});
     while(!q.empty()) {
         data now=q.top();q.pop();
         if(vis[now.x][now.y]) continue;
 
         vis[now.x][now.y]=;
         for(int i=;i<;i++) {
             int tox=now.x+tx[i],toy=now.y+ty[i];
             if(tox>n||toy>m||tox<||toy<) continue;
             if(dis[tox][toy]>now.val+(s[tox][toy]==''?:)) {
                 dis[tox][toy]=now.val+(s[tox][toy]==''?:);
                 q.push((data){tox,toy,dis[tox][toy]});
             }
         }
     }
 }
 int main() {
     int ans=;
     scanf("%d%d%d",&n,&m,&T);
     for(int i=;i<=n;i++) scanf("%s",s[i]+);
     for(int i=;i<=n;i++) {
         for(int j=;j<=m;j++) {
             dij(i,j);
             int mx=;
             for(int k=;k<=n;k++)
                 for(int w=;w<=m;w++)
                     if(dis[k][w]<=T) mx=max(mx,(i-k)*(i-k)+(j-w)*(j-w));
             ans=max(ans,mx);
         }
     }
     double tmp=sqrt(ans);
     printf("%.6lf",tmp);
 }