BZOJ 1295 Scoi2009 最长距离

1295: [SCOI2009]最长距离

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 1701  Solved: 921
[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 。

很水的一道最短路问题

由于n和m很小，我们可以暴力对于每一个点求出它到其他点所需要最小移除的障碍

然后如果到这个点所需要的移除障碍小于T的化，就用他们之间的距离来更新答案

#include <bits/stdc++.h>
#define ll long long
using namespace std;
inline int read(){
    int x=0;int f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
const int MAXN=1e6+10;
const int dx[4]={1,0,-1,0};
const int dy[4]={0,-1,0,1};
namespace zhangenming{
    struct node{
        int x,y;
    }q[MAXN];
    int head,tail,n,m,mapp[40][40],T,f[40][40],vis[40][40];
    double ans;
    inline void spfa(int xx,int yy){
        head=tail=0;memset(vis,0,sizeof(vis));
        memset(f,10,sizeof(f));
        vis[xx][yy]=1;
        f[xx][yy]=mapp[xx][yy];
        q[++tail].x=xx;q[tail].y=yy;
        while(head<tail){
            node tn=q[++head];
            for(int i=0;i<4;i++){
                int x=tn.x+dx[i];int y=tn.y+dy[i];
                if(x==0||y==0||x==n+1||y==m+1) continue;
                if(f[tn.x][tn.y]+mapp[x][y]<f[x][y]&&f[tn.x][tn.y]<=T){
                    f[x][y]=f[tn.x][tn.y]+mapp[x][y];
                    if(!vis[x][y]){
                        q[++tail].x=x;q[tail].y=y;
                        vis[x][y]=1;
                    }
                }
            }
            vis[tn.x][tn.y]=0;
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                if(f[i][j]<=T){
                    ans=max(ans,sqrt(1.0*((i-xx)*(i-xx)+(j-yy)*(j-yy))));
                }
            }
        }
    }
    void init(){
        n=read();m=read();T=read();
        for(int i=1;i<=n;i++){
            char ch[40];
            scanf("%s",ch);
            for(int j=0;j<m;j++){
                if(ch[j]=='0') mapp[i][j+1]=0;
                else mapp[i][j+1]=1;
            }
        }
    }
    void solve(){
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                spfa(i,j);
            }
        }
        cout<<setiosflags(ios::fixed)<<setprecision(6);
        cout<<ans<<endl;
    }
}
int main(){
    //freopen("All.in","r",stdin);
    //freopen("zhang.out","w",stdout);
    using namespace zhangenming;
    init();
    solve();
    return 0;
}

//对拍代码
#include <bits/stdc++.h>
#define ll long long
using namespace std;
inline int read(){
    int x=0;int f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int main(){
    freopen("All.in","w",stdout);
    srand(int(time(NULL)));
    int n=rand()%29+1;int m=rand()%29+1;int T=rand()%107;
    cout<<n<<' '<<m<<' '<<T<<endl;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            printf("%d",rand()%2);
        }
        printf("
");
    }
    return 0;
}