【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[]={,,,-},dy[]={,-,,};
 queue<int>q;
 char map[][];
 int n,en,m,K,first[],next[],v[],w[],dis[],num[][];
 bool inq[];
 double ans;
 void init(const int &s)
 {memset(dis,0x7f,sizeof(dis)); inq[s]=true; dis[s]=; 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]=='');
     for(int i=;i<n;i++)
       for(int j=;j<m;j++)
         if(dis[num[i][j]]<=K)
           ans=max(ans,Euclid_Dis(x,y,i,j));
     K+=(map[x][y]=='');
 }
 void BuildGraph()
 {
     for(int i=;i<n;i++)
       for(int j=;j<m;j++)
         for(int k=;k<;k++)
           {
             int tx=i+dx[k],ty=j+dy[k];
             if(tx>=&&tx<n&&ty>=&&ty<n)
               if(map[tx][ty]=='') AddEdge(num[i][j],num[tx][ty],);
               else AddEdge(num[i][j],num[tx][ty],);
           }
 }
 int main()
 {
     scanf("%d%d%d",&n,&m,&K);
     for(int i=;i<n;i++)
       scanf("%s",map[i]);
     for(int i=;i<n;i++)
       for(int j=;j<m;j++)
         num[i][j]=++en;
     en=;
     BuildGraph();
     for(int i=;i<n;i++)
       for(int j=;j<m;j++)
         {
           spfa(num[i][j]);
           calc(i,j);
         }
     printf("%.6lf\n",ans);
     return ;
 }