[CQOI2009]DANCE跳舞（ISAP写法）

https://daniu.luogu.org/problemnew/show/3153

#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>

using namespace std;

#define N 2501
#define M 3001
#define inf 2e9

int n,k;
char s[];
bool mp[][];

int tot;
int front[N],to[M<<],nxt[M<<],val[M<<],from[M<<];

int src,decc;

int cur[N];

int path[N],num[N],lev[N];

void read(int &x)
{
    x=; char c=getchar();
    while(!isdigit(c))  c=getchar();
    while(isdigit(c)) { x=x*+c-''; c=getchar();  }
}

void add(int u,int v,int w)
{
    to[++tot]=v; nxt[tot]=front[u]; front[u]=tot; val[tot]=w; from[tot]=u;
    to[++tot]=u; nxt[tot]=front[v]; front[v]=tot; val[tot]=; from[tot]=v;
    //cout<<u<<' '<<v<<'\n';
}

void rebuild(int mid)
{
    tot=;
    memset(front,,sizeof(front));
    for(int i=;i<=n;++i)
        for(int j=;j<=n;++j)
            if(mp[i][j]) add(i,n*+j,);
            else add(n+i,n*+j,);
    for(int i=;i<=n;++i) add(src,i,mid);
    for(int i=;i<=n;++i) add(n*+i,decc,mid);
    for(int i=;i<=n;++i) add(i,n+i,k);
    for(int i=;i<=n;++i) add(n*+i,n*+i,k);
}

bool bfs()
{
    queue<int>q;
    for(int i=src;i<=decc;++i) lev[i]=decc;
    q.push(decc);
    lev[decc]=;
    int now,t;
    while(!q.empty())
    {
        now=q.front();
        q.pop();
        for(int i=front[now];i;i=nxt[i])
        {
            t=to[i];
            if(lev[t]==decc && val[i^])
            {
                lev[t]=lev[now]+;
                q.push(t);
            }
        }
    }
    return lev[src]!=decc;
}

int augment()
{
    int now=decc,flow=inf;
    int i;
    while(now!=src)
    {
        i=path[now];
        flow=min(flow,val[i]);
        now=from[i];
    }
    now=decc;
    while(now!=src)
    {
        val[path[now]]-=flow;
        val[path[now]^]+=flow;
        now=from[path[now]];
    }
    return flow;
}

int max_flow(int mid)
{
    int flow=;
    if(!bfs()) return ;
    memset(num,,sizeof(num));
    for(int i=src;i<=decc;++i) num[lev[i]]++,cur[i]=front[i];
    int now=src; int t;
    while(lev[src]<=decc)
    {
        if(now==decc)
        {
            flow+=augment();
            now=src;
        }
        int advanced=false;
        for(int i=cur[now];i;i=nxt[i])
        {
            t=to[i];
            if(val[i]> && lev[t]==lev[now]-)
            {
                advanced=true;
                path[t]=i;
                cur[now]=i;
                now=t;
                break;
            }
        }
        if(!advanced)
        {
            int m=decc-;
            for(int i=front[now];i;i=nxt[i])
                if(val[i]>) m=min(m,lev[to[i]]);
            if(!--num[lev[now]]) break;
            num[lev[now]=m+]++;
            cur[now]=front[now];
            if(now!=src) now=from[path[now]];
        }
    }
    return flow;
}

bool check(int mid)
{
    rebuild(mid);
    return max_flow(mid)==n*mid;
}

int main()
{
    read(n); read(k);
    decc=n*+;
    for(int i=;i<=n;++i)
    {
        scanf("%s",s+);
        for(int j=;j<=n;++j) mp[i][j]=s[j]=='Y' ? true : false;
    }
    int l=,r=n,mid,ans;
    while(l<=r)
    {
        mid=l+r>>;
        if(check(mid)) ans=mid,l=mid+;
        else r=mid-;
    }
    cout<<ans;
}