bzoj1305[CQOI2009]dance跳舞

网络流是真的难想到正解。。

%po姐姐的题解：
将一个人拆成两个点，点1向点2连一条流量为k的边，两个人若互相喜欢则点1之间连边，不喜欢则点2之间连边
对于每一个要验证的x值 将每个人的点1向源或汇连一条流量为x的边
然后二分答案跑最大流即可

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define S 0
#define T 1001
using namespace std;
const int N=3e5+5;
const int M=2e3+5;
const int inf=1e9;
int head[N],go[N],next[N],val[N],tot=1;
int mp[M][M];
int n,m;
int q[N],dis[M],ans;
inline bool bfs()
{
    int t=0,w=1;
    //int w,used=0;
    memset(dis,-1,sizeof(dis));
    dis[0]=0;
    q[1]=0;
    while (t<w)
    {
        int x=q[++t];
        for(int i=head[x];i;i=next[i])
        {
            int v=go[i];
            if(val[i]&&dis[v]==-1)
            {
                dis[v]=dis[x]+1;
                q[++w]=v;
            }
        }
    }
    return dis[T]!=-1;
}
inline int dfs(int x,int f)
{
    if (x==T)return f;
    int w,used=0;
    for(int i=head[x];i;i=next[i])
    {
        int v=go[i];
        if (val[i]&&dis[v]==dis[x]+1)
        {
            w=f-used;
            w=dfs(v,min(w,val[i]));
            val[i]-=w;
            val[i^1]+=w;
            used+=w;
            if (used==f)return f;
        }
    }
    if (!used)dis[x]=-1;
    return used;
}
inline void dinic()
{
    while(bfs())ans+=dfs(0,inf);
}
inline void add(int x,int y,int z)
{
    go[++tot]=y;
    val[tot]=z;
    next[tot]=head[x];
    head[x]=tot;
}
inline void ins(int x,int y,int z)
{
    add(x,y,z);
    add(y,x,0);
}
inline void build(int x)
{
    tot=1;
    memset(head,0,sizeof(head));
    fo(i,1,n)ins(0,i,x);
    fo(i,1,n)ins(i,i+500,m);
    fo(i,1,n)ins(i+n+500,n+i,m);
    fo(i,1,n)ins(n+i,T,x);
    fo(i,1,n)
    fo(j,1,n)
    if (mp[i][j])ins(i,n+j,1);
    else ins(i+500,n+j+500,1);
}
int main()
{
    scanf("%d%d",&n,&m);
    char ch[60];
    fo(i,1,n)
    {
        scanf("%s",ch+1);
        int k=m;
        fo(j,1,n)
        {
            if (ch[j]=='Y')mp[i][j]=1;
        }
    }

    int l=0,r=50;
    int mx=0;
    while (l<=r)
    {
        int mid=(l+r)>>1;
        build(mid);
        ans=0;
        dinic();
        if (ans>=n*mid)
        {
            mx=mid,l=mid+1;
        }
        else r=mid-1;
    }
    printf("%d\n",mx);
}