【spfa】bzoj1295 [SCOI2009]最长距离

题意：给你一个n*m的点阵、有些点是障碍，求一个欧几里得距离最大的点对(A,B)，使得在移走的障碍≤T的情况下，可以从A走到B。

建图，跑n*m次spfa，求出从 每个点 出发到 其他所有点 的 经过的障碍数。若这个值<=T，则可以用来更新答案。

#include<cstdio>
#include<queue>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int dx[]={0,0,1,-1},dy[]={1,-1,0,0};
queue<int>q;
char map[31][31];
int n,en,m,K,first[50001],next[50001],v[50001],w[50001],dis[50001],num[31][31];
bool inq[50001];
double ans;
void init(const int &s)
{memset(dis,0x7f,sizeof(dis)); inq[s]=true; dis[s]=0; q.push(s);}
void AddEdge(const int &U,const int &V,const int &W)
{v[++en]=V; w[en]=W; next[en]=first[U]; first[U]=en;}
int sqr(const int &x){return x*x;}
double Euclid_Dis(const int &x1,const int &y1,const int &x2,const int &y2)
{return sqrt(sqr(x1-x2)+sqr(y1-y2));}
void spfa(const int &s)
{
    init(s);
    while(!q.empty())
      {
          int cur=q.front();
          for(int i=first[cur];i;i=next[i])
            if(dis[v[i]]>dis[cur]+w[i])
              {
                dis[v[i]]=dis[cur]+w[i];
                if(!inq[v[i]])
                  {
                    q.push(v[i]);
                    inq[v[i]]=true;
                  }
              }
          q.pop(); inq[cur]=false;
      }
}
void calc(const int &x,const int &y)
{
    K-=(map[x][y]=='1');
    for(int i=0;i<n;i++)
      for(int j=0;j<m;j++)
        if(dis[num[i][j]]<=K)
          ans=max(ans,Euclid_Dis(x,y,i,j));
    K+=(map[x][y]=='1');
}
void BuildGraph()
{
    for(int i=0;i<n;i++)
      for(int j=0;j<m;j++)
        for(int k=0;k<4;k++)
          {
            int tx=i+dx[k],ty=j+dy[k];
            if(tx>=0&&tx<n&&ty>=0&&ty<n)
              if(map[tx][ty]=='1') AddEdge(num[i][j],num[tx][ty],1);
              else AddEdge(num[i][j],num[tx][ty],0);
          }
}
int main()
{
    scanf("%d%d%d",&n,&m,&K);
    for(int i=0;i<n;i++)
      scanf("%s",map[i]);
    for(int i=0;i<n;i++)
      for(int j=0;j<m;j++)
        num[i][j]=++en;
    en=0;
    BuildGraph();
    for(int i=0;i<n;i++)
      for(int j=0;j<m;j++)
        {
          spfa(num[i][j]);
          calc(i,j);
        }
    printf("%.6lf\n",ans);
    return 0;
}