BZOJ 1295 [SCOI2009]最长距离

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<cstdio>
#include<queue>
#include<cmath>
#include<cstring>
#define maxn 1000
#define abs(x) (((x)>0)?(x):(-(x)))
using namespace std;
const int dx[10]={0,0,0,1,-1};
const int dy[10]={0,1,-1,0,0};
char ch[31][31];
int head[maxn],next[maxn],point[maxn],value[maxn];
int dist[maxn],now;
void add(int x,int y,int v)
{
    next[++now] = head[x];
    head[x] = now;
    point[now] = y;
    value[now] = v;
}
int dis(int x,int y,int xx,int yy)
{
    int u = abs(x-xx),v = abs(y-yy);
    return u*u+v*v;
}
void spfa(int s)
{
    memset(dist,-1,sizeof(dist));
    queue<int>q;
    q.push(s);
    dist[s]=0;
    int visit[maxn];
    visit[s]=1;
    while(!q.empty())
    {
        int k = q.front();
        q.pop();
        visit[k]=0;
        for(int i=head[k];i;i=next[i])
        {
            int u = point[i];
            if(dist[u]==-1 || dist[u]>dist[k]+value[i])
            {
                dist[u]=dist[k]+value[i];
                if(!visit[u])
                {
                    q.push(u);
                    visit[u] = 1;
                }
            }
        }
    }
}
int main()
{
    int n , m , t;
    scanf("%d%d%d",&n,&m,&t);
    for(int i=1;i<=n;i++)scanf("%s",ch[i]+1);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            for(int k=1;k<=4;k++)
            {
                int xx = i+dx[k],yy = j + dy[k];
                if(xx<1 || yy<1 || xx >n || yy>m)continue;
                add((i-1)*m+j,(xx-1)*m+yy,ch[xx][yy]=='0'?0:1);
            }
        }
    }
    int ans = 0;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            int u = (i-1)*m+j;
            spfa(u);
            for(int x=1;x<=n;x++)
            {
                for(int y = 1;y<=m;y++)
                {
                    int v = (x-1)*m+y;
                    int uu = dis(i,j,x,y);
                    if(dist[(x-1)*m+y]+ch[i][j]-'0'<=t && ans < uu)
                    {
                        ans = uu;
                    }
                }
            }
        }
    }
    printf("%.6f\n",sqrt(ans));
    return 0;
}