[bzoj1295][SCOI2009]SPFA

1295: [SCOI2009]最长距离

这是我BZOJ第一个AC的!!!!!!!!!

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

用SPFA暴力枚举，算出到每个点需要移动的障碍物的个数，挺简单的，可以当做SPFA裸题做。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
const int N = 100005;
int n,m,T;
double ans = 0;
int dx[4] = {1,0,-1,0};
int dy[4] = {0,-1,0,1};
bool map[35][35],vis[35][35];
int dis[35][35];
struct spf{
    int x;int y;
}q[N];
void put( int x, int y ) {
    for( int i = 1; i <= n; i++ )
        for( int j = 1; j <= m; j++ ) 
            if( dis[i][j] <= T && sqrt((x-i)*(x-i)+(y-j)*(y-j)) > ans )
                ans = sqrt((x-i)*(x-i)+(y-j)*(y-j));
}
void spfa( int a, int b ) {
    int tail = 1,head = 1;
    memset(vis,0,sizeof(vis));
    memset(dis,127,sizeof(dis));
    vis[a][b] = 1;
    dis[a][b] = map[a][b];
    q[1].x = a;q[1].y = b;
    while( tail <= head) {
        int x = q[tail].x, y = q[tail].y; tail++;
        vis[x][y] = 0;
        for( int i = 0; i < 4; i++ )
            if( x+dx[i] && x+dx[i]<=n && y+dy[i] && y+dy[i]<=m ) {
                int xx = x+dx[i],yy = y+dy[i];
                if( map[xx][yy] && dis[x][y]+1 < dis[xx][yy] ) {
                    dis[xx][yy] = dis[x][y]+1;
                    if( !vis[xx][yy] ) {
                        q[++head].x = xx;
                        q[head].y = yy;
                        vis[xx][yy] = 1;
                    }
                } else if( !map[xx][yy] && dis[x][y] < dis[xx][yy] ){
                    dis[xx][yy] = dis[x][y];
                    if( !vis[xx][yy] ) {
                        q[++head].x = xx;
                        q[head].y = yy;
                        vis[xx][yy] = 1;
                    }
                }
            }
    }
    put( a, b );
}
int main() {
    scanf("%d%d%d", &n, &m, &T);
    string s;
    for( int i = 1; i <= n; i++ ) {
        cin>>s;
        for( int j = 0; j < m; j++ )
            map[i][j+1] = s[j]-'0';
    }
    for( int i = 1; i <= n; i++ )
        for( int j = 1; j <= m; j++ )
            spfa( i, j );
    printf("%.6lf", ans);
    return 0;
}