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;
}