[SCOI2009]最长距离

题目

BZOJ 1295: [SCOI2009]最长距离

分析

本来以为是一道DP，想了很久都想不到咋做。于是考虑最短路，发现一些操作。
将有障碍的格子的权值标为1，没有障碍为0；
对于一个起点，跑一遍最长路，寻找$dis[i][j]≤t$的最大长度的(i,j)；
格子(i,j)即为距离起点的最远点，求出欧几里得距离；
枚举所有的格子作为起点，不断更新答案最大值。
对于每一个起点，只需一次spfa，然后枚举各点求得当前最远距离，复杂度为$O(n^2\cdot m^2)$。

代码

#include<cstdio>
#include<queue>
#include<cmath>
#include<cstring>
#include<iostream>
using namespace std;
const int maxn=33,inf=2147483647;
struct pos{
    int x,y;
};
int a[maxn][maxn],dis[maxn][maxn];
bool vis[maxn][maxn];
//方向数组，个人习惯，用于向四周扩展：
int dx[]={0,1,0,-1};
int dy[]={1,0,-1,0};
queue<pos>q;
int n,m,t;
double ans=0.0;

double calc(int i,int j,int x,int y)
{
    return (double)sqrt((x-i)*(x-i)+(y-j)*(y-j));
}

void spfa(int x,int y)
{
//每次spfa都要初始化：
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            dis[i][j]=inf,vis[i][j]=0;
    dis[x][y]=a[x][y],vis[x][y]=1;
    q.push(pos{x,y});
//最长路：
    while(!q.empty())
    {
        pos now=q.front();
        q.pop();
        int i=now.x,j=now.y;
        vis[i][j]=0;
        for(int k=0;k<4;k++)
        {
            int u=i+dx[k],v=j+dy[k];
            if(u<1||u>n||v<1||v>m)  continue;
            if(dis[u][v]>dis[i][j]+a[u][v])
            {
                dis[u][v]=dis[i][j]+a[u][v];
                if(!vis[u][v])
                {
                    q.push(pos{u,v});
                    vis[u][v]=1;
                }
            }
        }
    }
//求最远距离 
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            if(dis[i][j]<=t)
                ans=max(ans,calc(x,y,i,j)); 
}

int main()
{
    scanf("%d%d%d",&n,&m,&t);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            char c=getchar();
            while(c!='1'&&c!='0')   c=getchar();
            if(c=='1')  a[i][j]=1;
        }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            spfa(i,j);
    printf("%.6lf",ans);
    return 0;
}