P3701 「伪模板」主席树-（网络流最大流）

题目链接：点击进入

思路

对于每个byx的人，从源点向人连边，容量为此人的寿命。
对于每个诗乃酱的人，从人向汇点连边，容量为此人的寿命。
注：对于主席，他的寿命要多加上sum，sum为本方膜法师的人数,代表续命。
对于每个byx的人与诗乃酱的人，如果byx能够用此人赢得诗乃酱的人，从byx的这个人向诗乃酱的这个人连一条边，容量为1。
最后跑裸的最大流即可。

代码

#include<cstdio>
#include<algorithm>
#include<queue>
#include<iostream>
#include<cstring>
using namespace std;
typedef long long ll;
const int maxn=1e6+10;
struct node
{
    int to;
	int next;
	ll val;
}edge[maxn];
struct Node
{
	string s;
	ll life;
}p1[maxn],p2[maxn];
int head[maxn],tot;
void add(int u,int v,int val)
{
    edge[tot].to=v;
	edge[tot].val=val;
    edge[tot].next=head[u];
	head[u]=tot++;
}
int deep[maxn];
int cur[maxn];
ll n,m,sum1,sum2;
ll res;
int s;
int t;
int dfs(int pos,int flow)
{
    if(pos==t) return flow;
    int rest=flow,k,i;
    for(i=cur[pos];i!=-1&&rest;i=edge[i].next)
    {
    	int v=edge[i].to;
        int val=edge[i].val;
        if(deep[v]==deep[pos]+1&&val)
        {
            k=dfs(v,min(val,rest));
            if(!k) deep[v]=0;
            edge[i].val-=k;
            edge[i^1].val+=k;
            rest-=k; 
        }
	}
	cur[pos]=i;
    return flow-rest;
}
bool bfs()
{
	memset(deep,0,sizeof(deep));
    deep[s]=1;
    queue<int>q; 
    q.push(s);
    cur[s]=head[s];
    while(q.size())
    {
        int now=q.front();
		q.pop();
		for(int i=head[now];i!=-1;i=edge[i].next)
		{
			int v=edge[i].to;
            if(edge[i].val&&!deep[v])
            {
            	q.push(v);
            	cur[v]=head[v];
                deep[v]=deep[now]+1;
				if(v==t) return 1;
            }
		}
    }
    return 0;
}
void dinic()
{
	while(bfs())
        res+=dfs(s,0x3f3f3f3f);
}
int main()
{
	memset(head,-1,sizeof(head));
    cin>>n>>m;
    for(int i=1;i<=n;i++)
	{
		cin>>p1[i].s;
		if(p1[i].s[0]=='Y') sum1++;
	} 
    for(int i=1;i<=n;i++)
    {
    	cin>>p2[i].s;
    	if(p2[i].s[0]=='Y') sum2++;
	}
    for(int i=1;i<=n;i++)
    {
    	cin>>p1[i].life;
    	if(p1[i].s[0]=='J') p1[i].life+=sum1;
	}
    for(int i=1;i<=n;i++)
    {
    	cin>>p2[i].life;
    	if(p2[i].s[0]=='J') p2[i].life+=sum2;
	}
	s=2*n+1,t=2*n+2;
	for(int i=1;i<=n;i++)
	{
		add(s,i,p1[i].life);
		add(i,s,0);
	}
	for(int i=1;i<=n;i++)
	{
		add(i+n,t,p2[i].life);
		add(t,i+n,0);
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			if(p1[i].s[0]=='W'&&(p2[j].s[0]=='Y'||p2[j].s[0]=='E'))
				add(i,j+n,1),add(j+n,i,0);
			if(p1[i].s[0]=='J'&&(p2[j].s[0]=='H'||p2[j].s[0]=='W'))
				add(i,j+n,1),add(j+n,i,0);
			if(p1[i].s[0]=='E'&&(p2[j].s[0]=='Y'||p2[j].s[0]=='J'))
				add(i,j+n,1),add(j+n,i,0);
			if(p1[i].s[0]=='Y'&&(p2[j].s[0]=='H'||p2[j].s[0]=='J'))
				add(i,j+n,1),add(j+n,i,0);
			if(p1[i].s[0]=='H'&&(p2[j].s[0]=='W'||p2[j].s[0]=='E'))
				add(i,j+n,1),add(j+n,i,0);
		}
	}
    dinic();
    printf("%lld\n",min(res,m));
	return 0;
}