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