BZOJ3712 [PA2014]Fiolki

容易发现我们只需要求出每种反应的发生时间，然后排序做一遍即可

求反应的发生时间，就相当于一张图，开始每个点都是单独一个连通块，连通块会不断合并，求每条边的两个端点联通的时间

整体二分一下即可，每次判断现在是否联通吧边分到左右两边

用dfs版的话并查集回滚要按秩合并，多个log，所以写bfs版的就好

#include<iostream>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<cstdlib>
#include<cstdio>
#include<map>
#include<bitset>
#include<set>
#include<stack>
#include<vector>
#include<queue>
using namespace std;
#define MAXN 500010
#define MAXM 200010
#define ll long long
#define eps 1e-8
#define MOD 1000000007
#define INF 1000000000
struct data{
	int x;
	int y;
	int X;
	int Y;
	data(){
		
	}
	data(int _x,int _y,int _X=0,int _Y=0){
		x=_x;
		y=_y;
		X=_X;
		Y=_Y;
	}
};
data q[MAXN];
int hd,tl;
data th[MAXN];
int n,m,H;
int a[MAXN];
data e[MAXN];
data h[MAXN];
int f[MAXN];
ll ans;
int fa(int x){
	return f[x]==x?x:f[x]=fa(f[x]);
}
int main(){
	int i;
	scanf("%d%d%d",&n,&m,&H);
	for(i=1;i<=n;i++){
		scanf("%d",&a[i]);
	}
	for(i=1;i<=m;i++){
		scanf("%d%d",&e[i].x,&e[i].y);
	}
	for(i=1;i<=H;i++){
		scanf("%d%d",&h[i].x,&h[i].y);
		h[i].Y=i;
	}
	q[(tl%=MAXN)++]=data(1,m,1,H);
	int wzh=m+1;
	while(hd!=tl){
		data x=q[(hd%=MAXN)++];
		int l=x.x,r=x.y,L=x.X,R=x.Y;
		int mid=l+r>>1;
		if(wzh>mid){
			for(i=1;i<=n;i++){
				f[i]=i;
			}
			wzh=1;
		}
		while(wzh<=mid){
			f[fa(e[wzh].x)]=fa(e[wzh].y);
			wzh++;
		}
		int tl=L-1,tr=R+1;
		for(i=L;i<=R;i++){
			if(fa(h[i].x)==fa(h[i].y)){
				th[++tl]=h[i];
			}
		}
		for(i=R;i>=L;i--){
			if(fa(h[i].x)!=fa(h[i].y)){
				th[--tr]=h[i];
			}
		}
		memcpy(h+L,th+L,sizeof(data)*(R-L+1));
		if(l!=r){
			if(tl>=L){
				q[(::tl%=MAXN)++]=data(l,mid,L,tl);
			}
			if(tl<R){
				q[(::tl%=MAXN)++]=data(mid+1,r,tl+1,R);
			}
		}
	}
	for(i=1;i<=n;i++){
		f[i]=i;
	}
	for(i=1;i<=m;i++){
		f[fa(e[i].x)]=fa(e[i].y);
	}
	for(i=1;i<=H;i++){
		if(fa(h[i].x)==fa(h[i].y)){
			int t=min(a[h[i].x],a[h[i].y]);
			ans+=t*2;
			a[h[i].x]-=t;
			a[h[i].y]-=t;
		}
	}
	printf("%lld\n",ans);
	return 0;
}

/*
4 3 2
27 46 75 69 
4 2
3 2
2 1
1 3
3 4



*/