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
; int len,first
;
void ins(int x,int y){a[++len]=(E){y,first[x]},first[x]=len;}
int g
,f
[20],dep
,in
,tf
;
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;
}