[BZOJ1305] dance 最大流+二分

据说这是一种套路qwq（雾
首先拆点，把每一个人拆成‘喜欢点’和‘不喜欢点’，从每个男生的‘喜欢点’向自己的‘不喜欢点’连接容量为k的边，从每个女生的‘不喜欢点’向每个女生连接容量为k的边。对于男生i喜欢男生j，由男生i的‘喜欢点’向女生j的‘喜欢点’连接容量为1的点（因为两人只能跳1首歌）。对于男生i不喜欢女生j，由男生i的’不喜欢点’向女生j的‘不喜欢点’连接一条容量为1的边。
接下来，从源点向每个男生的‘喜欢点’连接一条容量为a的边（表示跳a首曲子），从每个女生的‘喜欢点’向汇点连接一条容量为a的点。a为最小的不满流时的，所以答案时a-1。（不满流时即跳不出a首歌）二分a即可

#include<bits/stdc++.h>
#define maxn 1010
#define inf 0x3f3f3f3f
using namespace std;
struct node{
    int to,f,m;
    node*next,*rev;
}*con[maxn];
char sh[55][55];
int n,k,s,t;
void addedge(int x,int y,int m)
{
    node*p=new node;
    p->to=y;
    p->f=0;p->m=m;
    p->next=con[x];
    con[x]=p;
    p=new node;
    p->to=x;
    p->f=p->m=0;
    p->next=con[y];
    con[y]=p;
    con[x]->rev=con[y];
    con[y]->rev=con[x];
}
int dl[maxn],de[maxn],head,tail;
bool bfs()
{
    bool tmp=false;
    memset(de,-1,sizeof(de));
    head=tail=1;
    de[s]=1;dl[tail]=s;
    while(head<=tail)
    {
        int v=dl[head++];
        if(v==t) tmp=true;
        for(node*p=con[v];p;p=p->next)
            if(p->f<p->m&&de[p->to]==-1)
        de[p->to]=de[v]+1,dl[++tail]=p->to;
    }
    return tmp;
}
int dinic(int v,int e)
{
    int temp=0,o=0;
    if(de[v]==-1) return 0;
    if(v==t) return e;
    for(node*p=con[v];p;p=p->next)
    {
    if(de[p->to]==de[v]+1&&p->f<p->m)
    {
        o=dinic(p->to,min(e-temp,p->m-p->f));
        temp+=o;
        p->f+=o;
        p->rev->f-=o;
    }
    if(temp==e) break;
    }
    if(temp==0) de[v]=-1;
    return temp;
}
void build(int x)
{
    for(int i=0;i<=1000;++i) con[i]=NULL;
    for(int i=1;i<=n;++i)
        addedge(s,i,x),addedge(i+n,t,x);
    for(int i=1;i<=n;++i)
        addedge(i,i+300,k),addedge(i+300+n,i+n,k);
    for(int i=1;i<=n;++i)
        for(int j=1;j<=n;++j)
        if(sh[i][j]=='Y') addedge(i,j+n,1);
                    else addedge(i+300,j+n+300,1);
}

int check()
{
    int ans=0;
    while(bfs())
        ans+=dinic(s,inf);
    return ans;
}
void work()
{
    for(int i=1;i<=n;++i)
        scanf("%s",sh[i]+1);
    int l=0,r=50,ans;
    while(l<=r)
    {
        int mid=(l+r)>>1;
        build(mid);
        if(check()>=n*mid) ans=mid,l=mid+1;
        else r=mid-1;
    }
    cout<<ans;
}
int main()
{
    scanf("%d%d",&n,&k);
    s=0;t=1000;
    work();
    return 0;
}