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

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

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

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

  1 #include<iostream>

  2 #include<cstdio>

  3 #include<cstring>

  4 #include<algorithm>

  5 #include<queue>

  6 #include<vector>

  7 using namespace std;

  8 #define MAXN 55

  9 #define inf 1<<30

 10 #define FILL(a,b) memset(a,b,sizeof(a))

 11 

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

 13 

 14 int n,dir[4][2]={{-1,0},{1,0},{0,-1},{0,1}};

 15 int dist[MAXN][MAXN],dd[MAXN][MAXN];

 16 bool mark[MAXN][MAXN];

 17 char map[MAXN][MAXN];

 18 

 19 void spfa(int vs)

 20 {

 21     FILL(mark,false);

 22     for(int i=0;i<MAXN;i++)

 23         for(int j=0;j<MAXN;j++)dd[i][j]=inf;

 24     queue<pair<int,int> >que;

 25     que.push(make_pair(g[vs].first,g[vs].second));

 26     dd[g[vs].first][g[vs].second]=0;

 27     while(!que.empty()){

 28         pair<int,int>p=que.front();

 29         que.pop();

 30         mark[p.first][p.second]=false;

 31         for(int i=0;i<4;i++){

 32             int x=p.first+dir[i][0],y=p.second+dir[i][1];

 33             if(x>=0&&x<n&&y>=0&&y<n&&map[x][y]!='#'){

 34                 if(dd[p.first][p.second]+1<dd[x][y]){

 35                     dd[x][y]=dd[p.first][p.second]+1;

 36                     if(!mark[x][y]){

 37                         mark[x][y]=true;

 38                         que.push(make_pair(x,y));

 39                     }

 40                 }

 41             }

 42         }

 43     }

 44 }

 45 

 46 vector<int>gg[MAXN];

 47 void Build(int limit)

 48 {

 49     for(int i=0;i<MAXN;i++)gg[i].clear();

 50     for(int i=0;i<g.size();i++){

 51         for(int j=0;j<h.size();j++){

 52             if(dist[i][j]<=limit)gg[i].push_back(j);

 53         }

 54     }

 55 }

 56 

 57 int ly[MAXN];

 58 bool vis[MAXN];

 59 

 60 int dfs(int u)

 61 {

 62     for(int i=0;i<gg[u].size();i++){

 63         int v=gg[u][i];

 64         if(!vis[v]){

 65             vis[v]=true;

 66             if(ly[v]==-1||dfs(ly[v])){

 67                 ly[v]=u;

 68                 return 1;

 69             }

 70         }

 71     }

 72     return 0;

 73 }

 74 

 75 int MaxMatch()

 76 {

 77     int res=0;

 78     FILL(ly,-1);

 79     for(int i=0;i<g.size();i++){

 80         FILL(vis,false);

 81         res+=dfs(i);

 82     }

 83     return res;

 84 }

 85 

 86 int main()

 87 {

 88     int _case,t=1;

 89     scanf("%d",&_case);

 90     while(_case--){

 91         scanf("%d",&n);

 92         g.clear();

 93         h.clear();

 94         for(int i=0;i<n;i++)scanf("%s",map[i]);

 95         for(int i=0;i<n;i++){

 96             for(int j=0;j<n;j++){

 97                 if(map[i][j]=='G')g.push_back(make_pair(i,j));

 98                 else if(map[i][j]=='H')h.push_back(make_pair(i,j));

 99             }

100         }

101         for(int i=0;i<MAXN;i++)

102             for(int j=0;j<MAXN;j++)dist[i][j]=inf;

103         for(int i=0;i<g.size();i++){

104             spfa(i);

105             for(int j=0;j<h.size();j++){

106                 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);

107             }

108         }

109         int low=inf,high=0,mid,ans=inf;

110         for(int i=0;i<g.size();i++){

111             for(int j=0;j<h.size();j++){

112                 if(dist[i][j]!=inf){

113                     low=min(low,dist[i][j]);

114                     high=max(high,dist[i][j]);

115                 }

116             }

117         }

118         while(low<=high){

119             mid=(low+high)>>1;

120             Build(mid);

121             if(MaxMatch()==(int)h.size()){

122                 ans=mid;

123                 high=mid-1;

124             }else

125                 low=mid+1;

126         }

127         printf("Case %d: ",t++);

128         if(ans!=inf){

129             printf("%d\n",ans);

130         }else

131             puts("Vuter Dol Kupokat");

132     }

133     return 0;

134 }

View Code