hdu5093(二分图最大匹配)

链接：点击打开链接

题意：给你一张 n*m 的矩阵,'*'表示海,'o'表示浮冰,'#'表示冰山,同一行或同一列只能放一个战舰。若两个战舰想放在同一行或同一列，那么它们之间必须要有冰山相隔。然后问在这张图上最多可以放多少个战舰

代码：

#include <queue>
#include <vector>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <iostream>
using namespace std;
const int INF=0x3f3f3f3f;
struct node{
int u,v,cap;
node(){}
node(int u,int v,int cap):u(u),v(v),cap(cap){}
}es[500005];
int R,S,T;
int dis[5005],iter[5005];
vector<int> tab[5005];
void addedge(int u, int v, int cap){
tab[u].push_back(R);
es[R++]=node(u,v,cap);
tab[v].push_back(R);
es[R++]=node(v,u,0);
}
int bfs(){
int i,h;
queue<int> q;
q.push(S);
memset(dis,INF,sizeof(dis));
dis[S]=0;
while(q.size()){
h=q.front();
q.pop();
for(i=0;i<tab[h].size();i++){
node &e=es[tab[h][i]];
if(e.cap>0&&dis[e.v]==INF){
dis[e.v]=dis[h]+1;
q.push(e.v);
}
}
}
return dis[T]<INF;
}
int dfs(int x,int maxflow){
int flow;
if(x==T)
return maxflow;
for(int &i=iter[x];i<tab[x].size();i++){
node &e=es[tab[x][i]];
if(dis[e.v]==dis[x]+1&&e.cap>0){
flow=dfs(e.v,min(maxflow,e.cap));
if(flow){
e.cap-=flow;
es[tab[x][i]^1].cap+=flow;
return flow;
}
}
}
return 0;
}
int dinic(){
int ans,flow;
ans=0;
while(bfs()){
memset(iter,0,sizeof(iter));
while(flow=dfs(S,INF))
ans+=flow;
}
return ans;
}                                               //网络流跑二分图最大匹配
char s[55][55];
int x[55][55],y[55][55]                                     ;
int main(){
int t,i,j,k,n,m,num1,num2;
scanf("%d",&t);
while(t--){
scanf("%d%d",&n,&m);
for(i=0;i<n;i++)
for(j=0;j<m;j++)
cin>>s[i][j];
memset(x,0,sizeof(x));
memset(y,0,sizeof(y));
num1=num2=1;
for(i=0;i<n;i++)                        //拆点，相当于将满足行的点进行编号和满足
for(j=0;j<m;j++){                       //列的点进行编号
if(s[i][j]=='*'&&!x[i][j]){
for(k=j;k<m;k++){
if(s[i][k]=='*')
x[i][k]=num1;
else if(s[i][k]=='#')           //只有有#出现时才可以增加满足的点
break;
}
num1++;
}
}
for(i=0;i<m;i++)
for(j=0;j<n;j++){
if(s[j][i]=='*'&&!y[j][i]){
for(k=j;k<n;k++){
if(s[k][i]=='*')
y[k][i]=num2;
else if(s[k][i]=='#')
break;
}
num2++;
}                                   //*ooo      1ooo   1ooo
}                                       //o###拆成->o###和 o###
num1--,num2--;                          //**#*      22#3   23#4
S=R=0,T=num1+num2+1;                    //ooo*      ooo4   ooo4
for(i=0;i<=T;i++)
tab[i].clear();
for(i=1;i<=num1;i++)
addedge(S,i,1);
for(i=1;i<=num2;i++)
addedge(i+num1,T,1);
for(i=0;i<n;i++)
for(j=0;j<m;j++){
if(x[i][j]&&y[i][j])
addedge(x[i][j],y[i][j]+num1,1);
}                                       //建二分图跑网络流
printf("%d\n",dinic());
}
return 0;
}