hdu 5093 最大流或二分图

点击打开链接

题意：给个m行n列的矩阵，其中#代表冰山，*代表大海，o代表浮冰，现在要停一些船，船只能停在大海上，但是每一行每一列最多可以停一艘船，但是如果中间有冰山隔离的话，一行或一列可以停靠多艘，问最多可以停多少艘船

思路：之前做过类似的，但是没有冰山这个限制，那个二分图直接用源点连行，列连汇点就行，这道题用二分图应该才是正解，我看了题意之后，感觉不到二分图的存在，根本没往那边想，但是还是过了，不知道有没有神犇之处错误，不过我写的最大流非常愚蠢而已，建边也比较麻烦......这里就不想多说了，有兴趣的看代码把

#include <queue>
#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int maxn=10010;
struct edge{
    int to,cap,rev;
    edge(){}
    edge(int a,int b,int c){to=a;cap=b;rev=c;}
};
vector<edge> G[maxn];
int level[maxn],iter[maxn];
void add_edge(int from,int to,int cap){
    G[from].push_back(edge(to,cap,G[to].size()));
    G[to].push_back(edge(from,0,G[from].size()-1));
}
void bfs(int s){
    memset(level,-1,sizeof(level));
    queue<int>que;
    level[s]=0;que.push(s);
    while(!que.empty()){
        int v=que.front();que.pop();
        for(unsigned int i=0;i<G[v].size();i++){
            edge &e=G[v][i];
            if(e.cap>0&&level[e.to]<0){
                level[e.to]=level[v]+1;
                que.push(e.to);
            }
        }
    }
}
int dfs(int v,int t,int f){
    if(v==t) return f;
    for(int &i=iter[v];i<G[v].size();i++){
        edge &e=G[v][i];
        if(e.cap>0&&level[v]<level[e.to]){
            int d=dfs(e.to,t,min(f,e.cap));
            if(d>0){
                e.cap-=d;
                G[e.to][e.rev].cap+=d;
                return d;
            }
        }
    }
    return 0;
}
int max_flow(int s,int t){
    int flow=0;
    while(1){
        bfs(s);
        if(level[t]<0) return flow;
        memset(iter,0,sizeof(iter));
        int f;
        while((f=dfs(s,t,inf))>0) flow+=f;
    }
}
char str[60][60];
int num[60][60],cnt[60];
int main(){
    int T,n,m;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        for(int i=1;i<=m;i++) scanf("%s",str[i]+1);
        for(int i=0;i<maxn;i++) G[i].clear();
        for(int i=1;i<=m;i++){
            int k=0;
            num[i][k++]=0;
            for(int j=1;j<=n;j++){
                if(str[i][j]=='#') num[i][k++]=j;
            }
            num[i][k++]=n+1;
            cnt[i]=k;
        }
        int k=1;
        for(int i=1;i<=m;i++){
            for(int j=1;j<cnt[i];j++){
                add_edge(0,k,1);
                for(int l=num[i][j-1]+1;l<num[i][j];l++){
                    if(str[i][l]=='*')
                    add_edge(k,3000+(i-1)*n+l,1);
                }
                k++;
            }
        }
        for(int i=1;i<=n;i++){
            int k=0;
            num[i][k++]=0;
            for(int j=1;j<=m;j++){
                if(str[j][i]=='#') num[i][k++]=j;
            }
            num[i][k++]=m+1;
            cnt[i]=k;
        }
        int k1=6000;
        for(int i=1;i<=n;i++){
            for(int j=1;j<cnt[i];j++){
                add_edge(k1,8050,1);
                for(int l=num[i][j-1]+1;l<num[i][j];l++){
                    if(str[l][i]=='*')
                    add_edge(3000+(l-1)*m+i,k1,1);
                }
                k1++;
            }
        }
        printf("%d\n",max_flow(0,8050));
    }
    return 0;
}