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[4][2]={{-1,0},{1,0},{0,-1},{0,1}};
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=0;i<MAXN;i++)
        for(int j=0;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]=0;
    while(!que.empty()){
        pair<int,int>p=que.front();
        que.pop();
        mark[p.first][p.second]=false;
        for(int i=0;i<4;i++){
            int x=p.first+dir[i][0],y=p.second+dir[i][1];
            if(x>=0&&x<n&&y>=0&&y<n&&map[x][y]!='#'){
                if(dd[p.first][p.second]+1<dd[x][y]){
                    dd[x][y]=dd[p.first][p.second]+1;
                    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=0;i<MAXN;i++)gg[i].clear();
    for(int i=0;i<g.size();i++){
        for(int j=0;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=0;i<gg[u].size();i++){
        int v=gg[u][i];
        if(!vis[v]){
            vis[v]=true;
            if(ly[v]==-1||dfs(ly[v])){
                ly[v]=u;
                return 1;
            }
        }
    }
    return 0;
}

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

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