网络流入门 bzoj1305

当初以为网络流是金牌题，从来都没接触过，后来才知道原来网络流只是一个模板，真正难的是构图问题，而且如果一眼看出是模板题的话只要能构出来图问题就不大，大家的Dinic模板都差不多。来一道简单的构图问题吧。
Description
一次舞会有n个男孩和n个女孩。每首曲子开始时，所有男孩和女孩恰好配成n对跳交谊舞。每个男孩都不会和同一个女孩跳两首（或更多）舞曲。有一些男孩女孩相互喜欢，而其他相互不喜欢（不会“单向喜欢”）。每个男孩最多只愿意和k个不喜欢的女孩跳舞，而每个女孩也最多只愿意和k个不喜欢的男孩跳舞。给出每对男孩女孩是否相互喜欢的信息，舞会最多能有几首舞曲？

Input
第一行包含两个整数n和k。以下n行每行包含n个字符，其中第i行第j个字符为’Y’当且仅当男孩i和女孩j相互喜欢。

Output
仅一个数，即舞曲数目的最大值。

Sample Input
3 0

YYY

YYY

YYY

Sample Output
3

建一个S,T。S连接每个男孩，T连接每个女孩。每个男孩有两个点，喜欢和不喜欢，女孩也有两个点，喜欢和不喜欢，喜欢的点连接喜欢的点，不喜欢的点连接不喜欢的点，从源点出来的点的权值为答案，每个人和自己喜欢的属性点之间的权值为答案，和自己不喜欢的属性点之间的权值为k，然后和自己喜欢或不喜欢的人之间的权值为1。二分答案跑Dinic就行

#include<bits/stdc++.h>
using namespace std;
const int INF=0x3f3f3f3f;
const int maxn=1e5+5;
int head[maxn],ver[maxn],edge[maxn],Next[maxn],d[maxn];
int n,k,s,t,tot,maxflow;
queue<int>Q;

void add(int x,int y,int z) {
    ver[++tot]=y,edge[tot]=z,Next[tot]=head[x],head[x]=tot;
    ver[++tot]=x,edge[tot]=0,Next[tot]=head[y],head[y]=tot;
}

bool bfs() {
    memset(d,0,sizeof(d));
    while(Q.size()) Q.pop();
    Q.push(s);d[s]=1;
    while(!Q.empty()) {
        int x=Q.front();Q.pop();
        for(int i=head[x];i;i=Next[i]) {
            if(edge[i]&&!d[ver[i]]) {
                Q.push(ver[i]);
                d[ver[i]]=d[x]+1;
                if(ver[i]==t) return 1;
            }
        }
    }
    return 0;
}

int dinic(int x,int flow) {
    if(x==t) return flow;
    int rest=flow,k;
    for(int i=head[x];i&&rest;i=Next[i]) {
        if(edge[i]&&d[ver[i]]==d[x]+1) {
            k=dinic(ver[i],min(rest,edge[i]));
            if(!k) d[ver[i]]=0;
            edge[i]-=k;
            edge[i^1]+=k;
            rest-=k;
        }
    }
    return flow-rest;
}
char node[100][100];

void build(int mid) {
    tot=1;memset(head,0,sizeof(head));
    for(int i=1;i<=n;i++)add(s,i,mid);//S点向男孩建边
    for(int i=1;i<=n;i++)add(i,i+500,k);//男孩向男孩分出来的不喜欢的点建边
    for(int i=1;i<=n;i++)add(n+i+500,n+i,k);//女孩向女孩分出来的不喜欢的点建边
    for(int i=1;i<=n;i++)add(n+i,t,mid);//女孩向T点建边
    for(int i=1;i<=n;i++)
       for(int j=1;j<=n;j++)
          if(node[i][j]=='Y')add(i,n+j,1);//i和j相互喜欢，两人建边
          else add(i+500,n+j+500,1);//i和j相互不喜欢，两人建边
}

int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    cin>>n>>k;
    s=0;t=1001;
    for(int i=1;i<=n;i++) {
        for(int j=1;j<=n;j++){
            cin>>node[i][j];
        }
    }
    int l=0,r=50,ans=0;
    while(l<r) {
        int mid=(l+r)/2;
        build(mid);
        int flow=0;maxflow=0;
        while(bfs()) {
            while(flow=dinic(s,INF)) maxflow+=flow;
        }
        if(maxflow>=n*mid) ans=mid,l=mid+1;
        else r=mid;
    }
    cout<<ans<<endl;
    return 0;
}