n<=50,不用说些什么了吧,网络流呀。
怎么跑呢?注意到,两个人之间的关系有喜欢与不喜欢,所以可以想到拆点,然后用关系建边,然后每人的两个点之间建一条容量为k的边,源点连向男生,女生连向汇点,这样就可以了。
日常语文不好,画个图吧。
3 1
YYN
NYY
YNY
图丑233
能看吧。。。
之后如何计算答案呢?
显然满足二分性质。。
源汇点的边连容量为二分的mid,然后满足这些边满流即可。
#include<bits/stdc++.h>
#define inf ~0u>>1
#define N 205
#define M 6000
using namespace std;
int n,k,mx,ans,S,T;
int like[N+1][N+1];
char opt;
int first[N+1],nex[M+1],to[M+1],flow[M+1],siz;
int dis[N+1];
inline char nc()
{
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline char getc()
{
char c=nc();
for(;!(c<='Z'&&c>='A');)c=nc();
return c;
}
inline int read()
{
int x=0,b=1;
char c=nc();
for(;!(c<='9'&&c>='0');c=nc())if(c=='-')b=-1;
for(;c<='9'&&c>='0';c=nc())x=x*10+c-'0';
return x*b;
}
inline void add(int x,int y,int z)
{
nex[siz]=first[x];
first[x]=siz;
to[siz]=y;
flow[siz]=z;
siz++;
}
inline void insert(int x,int y,int z)
{
add(x,y,z),add(y,x,0);
}
inline void build(int val)
{
siz=0,memset(first,-1,sizeof(first));
for(int i=1;i<=n;i++)insert(S,i,val);
for(int i=n+1;i<=2*n;i++)insert(i,T,val);
for(int i=1;i<=n;i++)insert(i,i+2*n,k);
for(int i=n+1;i<=2*n;i++)insert(i+2*n,i,k);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if(like[i][j])insert(i,j+n,1);
else insert(i+2*n,j+3*n,1);
}
inline bool BFS(int s,int t)
{
queue <int> Q;
for(int i=S;i<=T;i++)dis[i]=inf;
dis[s]=0,Q.push(s);
while(!Q.empty())
{
int x=Q.front();Q.pop();
for(int i=first[x];i!=-1;i=nex[i])
{
int u=to[i];
if(flow[i]&&dis[x]+1<dis[u])
{
dis[u]=dis[x]+1;
if(u==t)return true;
Q.push(u);
}
}
}
if(dis[t]==inf)return false;
return true;
}
inline int DFS(int s,int f,int t)
{
if(s==t)return f;
int rest=f;
for(int i=first[s];i!=-1;i=nex[i])
{
int u=to[i];
if(flow[i]&&dis[u]==dis[s]+1&&rest)
{
int tmp=DFS(u,min(rest,flow[i]),t);
if(tmp==0)dis[u]=0;
flow[i]-=tmp;
flow[i^1]+=tmp;
rest-=tmp;
}
}
return f-rest;
}
inline int Dinic(int s,int t)
{
int ans=0;
while(BFS(s,t))
ans+=DFS(s,~0u>>1,t);
return ans;
}
int main()
{
//freopen("in.txt","r",stdin);
n=read(),k=read();
S=0,T=4*n+1;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
{
opt=getc();
if(opt=='Y')like[i][j]=1;
else like[i][j]=0;
}
int l=0,r=50;
while(l<=r)
{
int mid=(l+r)/2;
build(mid);
ans=Dinic(S,T);
if(ans>=n*mid)mx=mid,l=mid+1;
else r=mid-1;
}
cout<<mx;
return 0;
}
慢啊,其实建边可以不用每次重建,改一改源汇点的边就好了。