B1295 [SCOI2009]最长距离 最短路

就是一道最短路的裸题，直接跑spfa就行了。（spfa死了）

最后在答案处判断是否障碍物太多，然后就直接找最大值就行。

（数据特别水，我错误算法60）

题干：

Description

windy有一块矩形土地，被分为 N*M 块 * 的小格子。 有的格子含有障碍物。 如果从格子A可以走到格子B，那么两个格子的距离就为两个格子中心的欧几里德距离。 如果从格子A不可以走到格子B，就没有距离。 如果格子X和格子Y有公共边，并且X和Y均不含有障碍物，就可以从X走到Y。 如果windy可以移走T块障碍物，求所有格子间的最大距离。 保证移走T块障碍物以后，至少有一个格子不含有障碍物。
Input
输入文件maxlength.in第一行包含三个整数，N M T。 接下来有N行，每行一个长度为M的字符串，''表示空格子，''表示该格子含有障碍物。
Output
输出文件maxlength.out包含一个浮点数，保留6位小数。
Sample Input
【输入样例一】

【输入样例二】

【输入样例三】

Sample Output
【输出样例一】
1.414214
【输出样例二】
3.605551
【输出样例三】
2.828427
HINT
%的数据，满足  <= N,M <=  ；  <= T <=  。 %的数据，满足  <= N,M <=  ；  <= T <=  。 %的数据，满足  <= N,M <=  ；  <= T <=  。

代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<ctime>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;
#define duke(i,a,n) for(int i = a;i <= n;i++)
#define lv(i,a,n) for(int i = a;i >= n;i--)
#define clean(a) memset(a,0,sizeof(a))
const int INF =  << ;
typedef long long ll;
typedef double db;
template <class T>
void read(T &x)
{
    char c;
    bool op = ;
    while(c = getchar(), c < '' || c > '')
        if(c == '-') op = ;
    x = c - '';
    while(c = getchar(), c >= '' && c <= '')
        x = x *  + c - '';
    if(op) x = -x;
}
template <class T>
void write(T x)
{
    if(x < ) putchar('-'), x = -x;
    if(x >= ) write(x / );
    putchar('' + x % );
}
int n,m,t;
int dx[] = {,,-,};
int dy[] = {,-,,};
int dis[][],a[][],mp[][];
bool inq[][];
bool p[][];
char s[];
db ans = ;
struct node
{
    int x,y;
}q[];
void getans(int x,int y)
{
    duke(i,x,n)
    {
        duke(j,,m)
        {
            if(dis[i][j] <= t && (y - j) * (y - j) + (x - i) * (x - i) > ans)
            ans = (y - j) * (y - j) + (x - i) * (x - i);
        }
    }
}
void spfa(int x,int y)
{
    int nowx,nowy,t = ,w = ,nx,ny;
    q[].x = x;q[].y = y;
    clean(inq);
    memset(dis,,sizeof(dis));
    inq[x][y] = ;
    dis[x][y] = mp[x][y];
    while(t <= w)
    {
        nowx = q[t].x;nowy = q[t].y;
        t++;
        for(int i = ;i < ;i++)
        {
            nx = nowx + dx[i];
            ny = nowy + dy[i];
            if(nx > n || nx < x || ny > m || ny < )
            continue;
            if(dis[nowx][nowy] + mp[nx][ny] < dis[nx][ny])
            {
                dis[nx][ny] = dis[nowx][nowy] + mp[nx][ny];
                if(!inq[nx][ny])
                {
                    q[++w].x = nx;
                    q[w].y = ny;
                    inq[nx][ny] = ;
                }
            }
        }
        inq[nowx][nowy] = ;
    }
    getans(x,y);
}
int main()
{
    read(n);read(m);read(t);
    duke(i,,n)
    {
        scanf("%s",s);
        duke(j,,m - )
        mp[i][j + ] = s[j] - '';
    }
    duke(i,,n)
    {
        duke(j,,m)
        {
            spfa(i,j);
        }
    }
    printf("%.6lf",sqrt(ans));
    return ;
}
/*
3 3 0
001
001
110
*/