[bzoj1305]跳舞

考虑如果没有k个人，那么就是裸的二分答案+最大流
对于这k个人，再在原来的每一个点裂点，中间的流量为k，然后裂出来的点向所有不能匹配的点连边，再二分答案+最大流即可

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 #define N 205
 4 struct ji{
 5     int nex,to,len;
 6 }edge[N*N];
 7 queue<int>q;
 8 int E,n,k,head[N],work[N],d[N];
 9 char s[N][N];
10 void add(int x,int y,int z){
11     edge[E].nex=head[x];
12     edge[E].to=y;
13     edge[E].len=z;
14     head[x]=E++;
15     if (E&1)add(y,x,0);
16 }
17 bool bfs(){
18     memset(d,-1,sizeof(d));
19     q.push(0);
20     d[0]=0;
21     while (!q.empty()){
22         int k=q.front();
23         q.pop();
24         for(int i=head[k];i!=-1;i=edge[i].nex)
25             if ((edge[i].len)&&(d[edge[i].to]<0)){
26                 d[edge[i].to]=d[k]+1;
27                 q.push(edge[i].to);
28             }
29     }
30     return d[4*n+1]>=0;
31 }
32 int dfs(int k,int s){
33     if (k>4*n)return s;
34     for(int &i=work[k];i!=-1;i=edge[i].nex)
35         if ((edge[i].len)&&(d[edge[i].to]==d[k]+1)){
36             int p=dfs(edge[i].to,min(s,edge[i].len));
37             if (p){
38                 edge[i].len-=p;
39                 edge[i^1].len+=p;
40                 return p;
41             }
42         }
43     return 0;
44 }
45 int dinic(){
46     int k,ans=0;
47     while (bfs()){
48         memcpy(work,head,sizeof(head));
49         while (k=dfs(0,0x3f3f3f3f))ans+=k;
50     }
51     return ans;
52 }
53 bool pd(int mid){
54     memset(head,-1,sizeof(head));
55     E=0;
56     for(int i=1;i<=n;i++){
57         add(0,i,mid);
58         add(i,i+n,k);
59         add(i+2*n,4*n+1,mid);
60         add(i+3*n,i+2*n,k);
61     }
62     for(int i=1;i<=n;i++)
63         for(int j=0;j<n;j++)
64             if (s[i][j]=='Y')add(i,j+2*n+1,1);
65             else add(i+n,j+3*n+1,1);
66     return dinic()==mid*n;
67 }
68 int main(){
69     scanf("%d%d",&n,&k);
70     for(int i=1;i<=n;i++)scanf("%s",s[i]);
71     int l=0,r=n;
72     while (l<r){
73         int mid=(l+r+1>>1);
74         if (pd(mid))l=mid;
75         else r=mid-1;
76     }
77     printf("%d",l);
78 }