【PA2014】【BZOJ3712】Fiolki

Description

化学家吉丽想要配置一种神奇的药水来拯救世界。
吉丽有n种不同的液体物质，和n个药瓶（均从1到n编号）。初始时，第i个瓶内装着g[i]克的第i种物质。吉丽需要执行一定的步骤来配置药水，第i个步骤是将第a[i]个瓶子内的所有液体倒入第b[i]个瓶子，此后第a[i]个瓶子不会再被用到。瓶子的容量可以视作是无限的。
吉丽知道某几对液体物质在一起时会发生反应产生沉淀，具体反应是1克c[i]物质和1克d[i]物质生成2克沉淀，一直进行直到某一反应物耗尽。生成的沉淀不会和任何物质反应。当有多于一对可以发生反应的物质在一起时，吉丽知道它们的反应顺序。每次倾倒完后，吉丽会等到反应结束后再执行下一步骤。
吉丽想知道配置过程中总共产生多少沉淀。

Input

第一行三个整数n,m,k(0<=m< n<=200000,0<=k<=500000)，分别表示药瓶的个数（即物质的种数），操作步数，可以发生的反应数量。
第二行有n个整数g[1],g[2],…,g[n]（1<=g[i]<=10^9)，表示初始时每个瓶内物质的质量。
接下来m行，每行两个整数a[i],bi，表示第i个步骤。保证a[i]在以后的步骤中不再出现。
接下来k行，每行是一对可以发生反应的物质c[i],di，按照反应的优先顺序给出。同一个反应不会重复出现。

Output

Sample Input

3 2 1

2 3 4

1 2

3 2

2 3
Sample Output

6
HINT

Source

鸣谢Jcvb

刚开始做的时候就记得省队集训时候TA爷秒了…说是LCA balabala
抱着Delayyy的pdf看了好久好久(虽然里面只有一句话而且还是启发式合并的做法和LCA并无关系
最后发现一种做法
倾倒时候我们如果对每个倾倒新建一个容器(就是新建一个点),那么如果有两个试剂反应了,他们反应的位置就是LCA上.
于是我们可以把能反应的试剂按反应处点的深度排序然后做一遍累积答案

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 500010
#define GET (ch>='0'&&ch<='9')
#define LL long long
using namespace std;
int n,m,k,top,Top;
int g[MAXN],to[MAXN],du[MAXN];
int fa[MAXN][20],deep[MAXN];
LL ans;
struct edge
{
    int to;
    edge *next;
}e[MAXN<<1],*prev[MAXN];
struct Edge
{
    int u,v,w,id;
    bool operator <(const Edge& a)const {return w==a.w?id<a.id:w>a.w;}
}E[MAXN];
void insert(int u,int v)    {e[++top].to=v;e[top].next=prev[u];prev[u]=&e[top];du[v]++;}
void Insert(int u,int v,int w,int id)   {if (!w)    return;E[++Top]=(Edge){u,v,deep[w],id};}
void dfs(int x,int f)
{
    for (int i=1;i<=19;i++)
        if ((1<<i)>deep[x]) break;
        else    fa[x][i]=fa[fa[x][i-1]][i-1];
    for (edge *i=prev[x];i;i=i->next)
        if (i->to!=f)   fa[i->to][0]=x,deep[i->to]=deep[x]+1,dfs(i->to,x);
}
int lca(int a,int b)
{
    if (deep[a]<deep[b])    swap(a,b);
    int t=deep[a]-deep[b];
    for (int i=0;i<=19;i++) if (t&(1<<i))   a=fa[a][i];
    for (int i=19;i>=0;i--) if (fa[a][i]!=fa[b][i]) a=fa[a][i],b=fa[b][i];
    if (a==b)   return a;   return  fa[a][0];
}
void in(int &x)
{
    char ch=getchar();x=0;
    while (!GET)    ch=getchar();
    while (GET) x=x*10+ch-'0',ch=getchar();
}
int main()
{
    in(n);in(m);in(k);int u,v,w,z;
    for (int i=1;i<=n;i++)  in(g[i]),to[i]=i;
    for (int i=1;i<=m;i++)  in(u),in(v),insert(i+n,to[u]),insert(i+n,to[v]),to[v]=i+n;
    for (int i=1;i<=n+m;i++)    if (!du[i]) dfs(i,0);
    for (int i=1;i<=k;i++)  in(u),in(v),z=lca(u,v),Insert(u,v,z,i);
    sort(E+1,E+Top+1);
    for (int i=1;i<=Top;i++)    w=min(g[E[i].u],g[E[i].v]),g[E[i].u]-=w,g[E[i].v]-=w,ans+=w;
    ans<<=1;cout<<ans<<endl;
}