首先,将问题换个问法。
即Ans支舞曲是否可以跳完。
如果Ans支舞曲可以跳完,那么匹配数==Ans*n
我们只需要不停地枚举Ans,直到跳不完为止,那么Ans-1即为答案。
建模分析:
首先 我们把一个人分为喜欢和不喜欢两个点。
从源点向男生喜欢连一条容量为Ans的边,模拟跳Ans支舞曲。
对于女生也是如此。
不喜欢上限为K的限制,我们只需要从不喜欢向喜欢边连一条容量为K的边即可。
男女生连边时,喜欢的连喜欢的,不喜欢的连不喜欢的,容量均为1。
跑出的最大流即为最大匹配数。
BZOJ
Luogu
#include <cstdio>
#include <iostream>
#include <queue>
#include <cstring>
#define il inline
using namespace std;
const int inf=0x7fffffff;
const int maxm=110000;
int head[maxm],to[maxm*2],cap[maxm*2],net[maxm*2],deep[maxm],cnt=1;
il void add(int x,int y,int c){cnt++,to[cnt]=y,cap[cnt]=c,net[cnt]=head[x],head[x]=cnt;}
queue <int> dl;
char like[90][90];
int s=0,t=1000;
il bool BFS(int s,int t)
{
while(!dl.empty()) dl.pop();
memset(deep,-1,sizeof(deep));
dl.push(s),deep[s]=0;
while(!dl.empty())
{
int x=dl.front();dl.pop();
for(int i=head[x];i;i=net[i])
if(cap[i]>0&&deep[to[i]]==-1)
dl.push(to[i]),deep[to[i]]=deep[x]+1;
}
return deep[t]==-1?0:1;
}
int dfs(int now,int flow,int t)
{
if(now==t) return flow;
int w,used=0;
for(int i=head[now];i;i=net[i])
{
int v=to[i];
if(deep[v]==deep[now]+1&&cap[i])
{
w=dfs(v,min(flow-used,cap[i]),t);
cap[i]-=w;
cap[i^1]+=w;
used+=w;
if(used==flow) return flow;
}
}
if(!used) deep[now]=-1;
return used;
}
il int dinic(int s,int t)
{
int maxflow=0;
while(BFS(s,t)) maxflow+=dfs(s,inf,t);
return maxflow;
}
inline void adx(int x,int y,int cax)
{
add(x,y,cax),add(y,x,0);
}
il int read()
{
int x=0,w=1;
char ch=0;
while(ch<'0'||ch>'9')
{
if(ch=='-') w=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
return x*w;
}
il void build(int n,int k,int lit)
{
memset(head,0,sizeof(head));
cnt=1;
for(int i=1;i<=n;i++) adx(s,i,lit),adx(i+2*n,t,lit),adx(i,i+n,k),adx(i+3*n,i+2*n,k);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
{
if(like[i][j]=='Y') adx(i,j+2*n,1);
else adx(i+n,j+3*n,1);
}
}
int main()
{
int n=read(),k=read();
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
cin>>like[i][j];
int ans=0;
for(;;)
{
ans++;
build(n,k,ans);
int mf=dinic(s,t);
if(ans*n!=mf) break;
}
printf("%d",ans-1);
return 0;
}