【luoguP3701】「伪模板」主席树

description

byx和诗乃酱都非常都非常喜欢种树。有一天，他们得到了两颗奇怪的树种，于是各自取了一颗回家种树，并约定几年后比一比谁种出来的树更加牛x。

很快，这棵树就开花结果了。byx和诗乃酱惊讶的发现，这是一棵主席树，树上长满了主席和主席的朋友们。这棵树上一共有五种人，主席(J)，记者(HK)，高人(W)，女王(E)和膜法师(YYY)。他们发现，他们的主席树上的人数相同，都为N。

研究发现，这五种人的输赢如上图所示（一样的人不能PK），箭头指向输的人。至于为什么，留给同学们自己思考。

比赛如期进行。

byx和诗乃酱要进行M场比赛，每一场比赛他们会选出树上的两个人来比较看谁更牛x。

第i个人寿命为Lifei秒，每次比完赛他们就会-1s。当他们生命为0s时他们就不能再比赛了。

同时，当J的寿命为0时，同一棵树上的YYY可以为他+1s。每个YYY只能给.每个J续一次。

那么问题来了

现在给定N,M（1≤N≤100,1≤M≤1000），诗乃酱和byx每一个人所属种类(J,HK,W,YYY或E)以及每一个人的生命，生命不超过50.请你算算A最多能够赢得多少场比赛呢。

数据保证每一场一定都有人用。两个人之间只能比一场。

---

analysis

• 我是本来想复习主席树结果被骗进去复习了一波dinic你们给我搞的这道题啊excited！

• 由题面可以构建网络流的模型，首先把两边膜法师的数量统计一下提前先给两边主席续命

• 每个人生命 （秒数） 为\(0\)就不可以再贡献，而且对战 （谈笑风生） 成功可以获得\(1\)的贡献

• 那么源点向\(A\)方连生命 （秒数） 的权值的边，\(B\)方向汇点连生命 （秒数） 的权值的边

• 中间能赢的关系就连\(1\)边，最后判断一下最大流和\(m\)哪个更小就吼了

---

code

#pragma GCC optimize("O3")
#pragma G++ optimize("O3")
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<queue>
#define MAXN 105
#define MAXM 20005
#define INF 1000000007
#define ll long long
#define reg register ll
#define fo(i,a,b) for (reg i=a;i<=b;++i)
#define fd(i,a,b) for (reg i=a;i>=b;--i)
#define rep(i,a) for (reg i=last[a];i;i=next[i])

using namespace std;

ll last[MAXM],next[MAXM],tov[MAXM],len[MAXM];
ll depth[MAXN*2],life[MAXN*2];
char a[MAXN*2],s[5];
ll n,m,S,T,tot=1,sum1,sum2;
queue<ll>q;

inline ll read()
{
	ll x=0,f=1;char ch=getchar();
	while (ch<'0' || '9'<ch){if (ch=='-')f=-1;ch=getchar();}
	while ('0'<=ch && ch<='9')x=x*10+ch-'0',ch=getchar();
	return x*f;
}
inline void link(ll x,ll y,ll z)
{
	next[++tot]=last[x],last[x]=tot,tov[tot]=y,len[tot]=z;
	next[++tot]=last[y],last[y]=tot,tov[tot]=x,len[tot]=0;
}
inline bool bfs()
{
	while (!q.empty())q.pop();
	memset(depth,0,sizeof(depth));
	q.push(S),depth[S]=1;
	while (!q.empty())
	{
		ll now=q.front();q.pop();
		rep(i,now)if (len[i] && !depth[tov[i]])depth[tov[i]]=depth[now]+1,q.push(tov[i]);
	}
	return depth[T];
}
inline ll dfs(ll x,ll flow)
{
	ll ans=0;
	if (x==T)return flow;
	rep(i,x)if (len[i] && depth[tov[i]]==depth[x]+1)
	{
		ll tmp=dfs(tov[i],min(flow-ans,len[i]));
		if (!tmp)depth[tov[i]]=0;
		ans+=tmp,len[i]-=tmp,len[i^1]+=tmp;
		if (ans==flow)break;
	}
	return ans;
}
inline ll dinic()
{
	ll ans=0;
	while (bfs())
	{
		while (ll x=dfs(S,INF))ans+=x;
	}
	return ans;
}
int main()
{
	freopen("P3701.in","r",stdin);
	n=read(),m=read(),S=2*n+1,T=2*n+2;
	fo(i,1,n)scanf("%s",&s),a[i]=s[0],sum1+=a[i]=='Y';
	scanf("\n");
	fo(i,n+1,n*2)scanf("%s",&s),a[i]=s[0],sum2+=a[i]=='Y';

	fo(i,1,n)life[i]=read()+(a[i]=='J'?sum1:0);
	fo(i,n+1,n*2)life[i]=read()+(a[i]=='J'?sum2:0);

	fo(i,1,n)link(S,i,life[i]);
	fo(i,n+1,n*2)link(i,T,life[i]);

	fo(i,1,n)fo(j,n+1,n*2)
	{
		if (a[i]=='J' && (a[j]=='W' || a[j]=='H'))link(i,j,1);
		if (a[i]=='W' && (a[j]=='Y' || a[j]=='E'))link(i,j,1);
		if (a[i]=='H' && (a[j]=='W' || a[j]=='E'))link(i,j,1);
		if (a[i]=='Y' && (a[j]=='H' || a[j]=='J'))link(i,j,1);
		if (a[i]=='E' && (a[j]=='Y' || a[j]=='J'))link(i,j,1);
	}
	printf("%lld\n",min(dinic(),m));
	return 0;
}