BZOJ 1305: [CQOI2009]dance跳舞

1305: [CQOI2009]dance跳舞

啊啊啊noi同步赛弃疗来刷水题了不知道会不会被老师打死

。。。

曾经用几行贪心水过vijos和luogu的题。。终于跪在了bzoj面前

显然答案具有单调性，考虑二分（暴力枚举也可以）

考虑对答案x的可行性验证：

对于每一个人拆成两个点，p,p'

p表示连向喜欢的，p'则是不喜欢的

那么对于喜欢的两个人连p，否则连p'

然后考虑k的限制，我们在同一个人的p和p'间连k的边

然后，S向男生连流量x的边，女生向t连流量x的边

然后跑dinic

如果满流，说明答案可行

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <queue>
#include <string>
#include <map>
#include <cstring>
#include <ctime>
#include <vector>
#define inf 1e9
#define ll long long
#define maxm 500010
#define maxn 1001
#define For(i,j,k) for(int i=j;i<=k;i++)
#define Dow(i,j,k) for(int i=k;i>=j;i--)
using namespace std;
inline ll read(){   ll x=0,f=1;char ch=getchar();   while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}  while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}  return x*f; }
inline void write(ll x){    if (x<0) putchar('-'),x=-x; if (x>=10) write(x/10);   putchar(x%10+'0');  }
inline void writeln(ll x){write(x);puts("");}
int dis[2001],q[2001],S,T,n,like[201][201],poi[200001],nxt[200001],v[200001],head[20001],ANS;
int k,ans,cnt;
char s[201];
inline void add(int x,int y,int z){poi[++cnt]=y;nxt[cnt]=head[x];v[cnt]=z;head[x]=cnt;poi[++cnt]=x;nxt[cnt]=head[y];v[cnt]=0;head[y]=cnt;}
inline bool bfs()
{
	For(i,1,T)	dis[i]=0;
	int l=1,r=1;
	q[1]=S;dis[S]=1;
	while(l<=r)
	{
		for(int i=head[q[l]];i;i=nxt[i])
			if(!dis[poi[i]]&&v[i])	dis[poi[i]]=dis[q[l]]+1,q[++r]=poi[i];
		l++;
	}
	return dis[T]!=0;
}
inline int dfs(int x,int flow)
{
	if(x==T)	return flow;
	int used=0;
	for(int i=head[x];i;i=nxt[i])
	{
		if(v[i]&&dis[poi[i]]==dis[x]+1)
		{
			int mi=min(v[i],flow-used);
			int tmp=dfs(poi[i],mi);
			v[i]-=tmp;v[i^1]+=tmp;used+=tmp;
			if(used==flow)	return flow;
		}
	}
	return used;
}
inline void dinic(){while(bfs())	ans+=dfs(S,inf);}
inline bool check(int x)
{
	S=4*n+1;T=S+1;
	cnt=1;
	For(i,1,T)	head[i]=0;
	For(i,1,n)
		For(j,1,n)
		{
			if(like[i][j])	add(i,j+3*n,1);
				else add(i+n,j+2*n,1);
		}
	For(i,1,n)	add(i,i+n,k),add(i+2*n,i+3*n,k);
	For(i,1,n)	add(S,i,x),add(3*n+i,T,x);
	ans=0;dinic();
	if(ans==x*n)	return 1;else return 0;
}
int main()
{
	n=read();k=read();
	For(i,1,n)
	{
		scanf("\n%s",s+1);
		For(j,1,n)
			if(s[j]=='Y')	like[i][j]=1;
	}	
	int l=1,r=n;
	while(l<=r)
	{
		int mid=(l+r)>>1;
		if(check(mid))	ANS=max(ANS,mid),l=mid+1;else r=mid-1;
	}
	writeln(ANS);
}