loj 1150(spfa预处理+二分+最大匹配)

题目链接：http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=26864

思路：首先是spfa预处理出每个'G'到'H'的最短距离，然后就是二分最大距离，最大匹配验证即可。

PS:这道题一开始没什么思路，然后想先最简单的spfa预处理写一下吧，然后写着写着就突然豁然开朗，有思路了！然后就AC了!

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<queue>
 #include<vector>
 using namespace std;
 #define MAXN 55
 #define inf 1<<30
 #define FILL(a,b) memset(a,b,sizeof(a))

 vector<pair<int,int> >g,h;

 int n,dir[][]={{-,},{,},{,-},{,}};
 int dist[MAXN][MAXN],dd[MAXN][MAXN];
 bool mark[MAXN][MAXN];
 char map[MAXN][MAXN];

 void spfa(int vs)
 {
     FILL(mark,false);
     for(int i=;i<MAXN;i++)
         for(int j=;j<MAXN;j++)dd[i][j]=inf;
     queue<pair<int,int> >que;
     que.push(make_pair(g[vs].first,g[vs].second));
     dd[g[vs].first][g[vs].second]=;
     while(!que.empty()){
         pair<int,int>p=que.front();
         que.pop();
         mark[p.first][p.second]=false;
         for(int i=;i<;i++){
             int x=p.first+dir[i][],y=p.second+dir[i][];
             if(x>=&&x<n&&y>=&&y<n&&map[x][y]!='#'){
                 if(dd[p.first][p.second]+<dd[x][y]){
                     dd[x][y]=dd[p.first][p.second]+;
                     if(!mark[x][y]){
                         mark[x][y]=true;
                         que.push(make_pair(x,y));
                     }
                 }
             }
         }
     }
 }

 vector<int>gg[MAXN];
 void Build(int limit)
 {
     for(int i=;i<MAXN;i++)gg[i].clear();
     for(int i=;i<g.size();i++){
         for(int j=;j<h.size();j++){
             if(dist[i][j]<=limit)gg[i].push_back(j);
         }
     }
 }

 int ly[MAXN];
 bool vis[MAXN];

 int dfs(int u)
 {
     for(int i=;i<gg[u].size();i++){
         int v=gg[u][i];
         if(!vis[v]){
             vis[v]=true;
             if(ly[v]==-||dfs(ly[v])){
                 ly[v]=u;
                 return ;
             }
         }
     }
     return ;
 }

 int MaxMatch()
 {
     int res=;
     FILL(ly,-);
     for(int i=;i<g.size();i++){
         FILL(vis,false);
         res+=dfs(i);
     }
     return res;
 }

 int main()
 {
     int _case,t=;
     scanf("%d",&_case);
     while(_case--){
         scanf("%d",&n);
         g.clear();
         h.clear();
         for(int i=;i<n;i++)scanf("%s",map[i]);
         for(int i=;i<n;i++){
             for(int j=;j<n;j++){
                 if(map[i][j]=='G')g.push_back(make_pair(i,j));
                 else if(map[i][j]=='H')h.push_back(make_pair(i,j));
             }
         }
         for(int i=;i<MAXN;i++)
             for(int j=;j<MAXN;j++)dist[i][j]=inf;
         for(int i=;i<g.size();i++){
             spfa(i);
             for(int j=;j<h.size();j++){
                 if(dd[h[j].first][h[j].second]!=inf)dist[i][j]=min(dist[i][j],*dd[h[j].first][h[j].second]+);
             }
         }
         int low=inf,high=,mid,ans=inf;
         for(int i=;i<g.size();i++){
             for(int j=;j<h.size();j++){
                 if(dist[i][j]!=inf){
                     low=min(low,dist[i][j]);
                     high=max(high,dist[i][j]);
                 }
             }
         }
         while(low<=high){
             mid=(low+high)>>;
             Build(mid);
             if(MaxMatch()==(int)h.size()){
                 ans=mid;
                 high=mid-;
             }else
                 low=mid+;
         }
         printf("Case %d: ",t++);
         if(ans!=inf){
             printf("%d\n",ans);
         }else
             puts("Vuter Dol Kupokat");
     }
     return ;
 }