BZOJ 3712: [PA2014]Fiolki

我的心好len痛啊。。
蒟蒻并不会O（n） lca。。。
其实把它弄成树的样子就好了 注意可能是森林（简直就是嘛…）
反应就发生在lca 通过深度 合并顺序 反应顺序 三个东西排个序就好了。。
谁说sort是稳定排序的。。。蓝瘦 香菇

如果觉得代码丑 不要问我为什么QAQ。。。

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=200005;
inline int read(){
    char ch=getchar(); int x=0,f=1;
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0'; ch=getchar();}
    return x*f;
}
struct E{int y,next;}a[N]; int len,first[N];
void ins(int x,int y){a[++len]=(E){y,first[x]},first[x]=len;}
int g[N],f[N][20],dep[N],in[N],tf[N];
void dfs(int x){
    for(int i=1;i<=18;i++)
        if(f[x][i-1])f[x][i]=f[f[x][i-1]][i-1];
        else break;
    for(int k=first[x],u=0;k;k=a[k].next){
        int y=a[k].y; if(y==f[x][0])continue;
        dep[y]=dep[x]+1,f[y][0]=x,tf[y]=++u,dfs(y);
    }
}
struct node{int x,y,d,o,id;}p[N/2*5];
void lca(int x,int y,int k){
    if(dep[x]<dep[y])x^=y^=x^=y;
    p[k].x=x,p[k].y=y,p[k].id=k;
    for(int i=18;~i;i--)
        if(dep[x]-dep[y]>=(1<<i))x=f[x][i];
    if(x==y){
        x=p[k].x;
        for(int i=18;~i;i--)
            if(dep[x]-dep[y]>(1<<i))x=f[x][i];
        p[k].o=tf[x],p[k].d=dep[y];
        return;
    }
    for(int i=18;~i;i--)
        if(f[x][i]!=f[y][i])x=f[x][i],y=f[y][i];
    if(f[x][0]!=f[y][0] || !f[x][0])p[k].d=-1;
    else p[k].d=dep[f[x][0]];
    p[k].o=min(tf[x],tf[y]);
}
bool Cmp(node x1,node x2){
    if(x1.d!=x2.d)return x1.d>x2.d;
    if(x1.o!=x2.o)return x1.o>x2.o;
    return x1.id<x2.id;
}
int main()
{
    int n=read(),m=read(),k=read(),i,x,y;
    for(i=1;i<=n;i++)g[i]=read();
    for(i=1;i<=m;i++)x=read(),y=read(),ins(y,x),in[x]++;
    for(i=1;i<=n;i++)if(!in[i])dfs(i);
    for(i=1;i<=k;i++)
        x=read(),y=read(),lca(x,y,i);
    sort(p+1,p+1+k,Cmp);
    LL ans=0;
    for(i=1;i<=k;i++){
        if(p[i].d<0)break;
        int r=min(g[p[i].x],g[p[i].y]);
        g[p[i].x]-=r,g[p[i].y]-=r,ans+=2ll*r;
    }
    printf("%lld\n",ans);
    return 0;
}