时间限制: 3 Sec 内存限制: 128 MB
题目描述
化学家吉丽想要配置一种神奇的药水来拯救世界。
吉丽有n种不同的液体物质，和n个药瓶（均从1到n编号）。初始时，第i个瓶内装着g[i]克的第i种物质。吉丽需要执行一定的步骤来配置药水，第i个步骤是将第a[i]个瓶子内的所有液体倒入第b[i]个瓶子，此后第a[i]个瓶子不会再被用到。瓶子的容量可以视作是无限的。
吉丽知道某几对液体物质在一起时会发生反应产生沉淀，具体反应是1克c[i]物质和1克d[i]物质生成2克沉淀，一直进行直到某一反应物耗尽。生成的沉淀不会和任何物质反应。当有多于一对可以发生反应的物质在一起时，吉丽知道它们的反应顺序。每次倾倒完后，吉丽会等到反应结束后再执行下一步骤。
吉丽想知道配置过程中总共产生多少沉淀。
输入
第一行三个整数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],b[i](1<=a[i],b[i]<=n,a[i]≠b[i])，表示第i个步骤。保证a[i]在以后的步骤中不再出现。
接下来k行，每行是一对可以发生反应的物质c[i],d[i](1<=c[i],d[i]<=n,c[i]≠d[i])，按照反应的优先顺序给出。同一个反应不会重复出现。
输出
配置过程中总共产生多少沉淀。
样例输入
3 2 1
2 3 4
1 2
3 2
2 3
样例输出
6
题解
题目并不复杂，而且有很多可以利用的性质，比如说每个反应只会进行一次，所有物质在遇到反应物之前并没有什么用之类的，但是还是只想到了n^2的并查集+链表做法，水了30。通过这道题学会了一种很神奇的东西叫重构树，即每倾倒一次就新建一个节点作为这两个节点的父亲，然后以此类推，可以想见每一个反应一定是在LCA上发生的，然后我们就可以欢快地按优先级处理每个点发生的反应，只要倒着加边反应就会按优先级从大到小进行。重构树作为一种新思路让这道题一下子简明了起来，不过我实现得有点繁冗，因为用的是离线LCA，再加上把每个反应倒序加到LCA上一共用了三个邻接表两个并查集。
#include<iostream> #include<cstdio> #include<cstring> using namespace std; const int sj1=200005; const int sj2=500005; const int sj3=sj1+sj2; int n,m,k,f[sj3],g[sj1],yz[sj2],yy[sj2],h[sj1],zx[sj3]; int tp,d[sj3],e1,e2,lca[sj2],fa[sj2],size,w[sj3]; long long ans; bool r[sj3]; int read() { int jg=0,jk=getchar()-'0'; if(jk<=9&&jk>0) jg=jk; jk=getchar()-'0'; while(jk<=9&&jk>=0) jg=jg*10+jk,jk=getchar()-'0'; return jg; } int find(int x) { if(f[x]==-1) return x; f[x]=find(f[x]); return f[x]; } void hb(int x,int y) { x=find(x),y=find(y); if(x!=y) f[x]=y; } int fi(int x) { if(w[x]==-1) return x; w[x]=fi(w[x]); return w[x]; } void hbb(int x,int y) { x=fi(x),y=fi(y); if(x!=y) w[x]=y; } struct tree { int v,ne,num; }t[sj2<<1]; struct B { int v,ne; }b[sj2<<2]; void ad1(int x,int y,int z) { t[e1].v=y,t[e1].ne=h[x]; t[e1].num=z,h[x]=e1++; } void ad2(int x,int y) { b[e2].v=y,b[e2].ne=d[x]; d[x]=e2++; } void init() { memset(d,-1,sizeof(d)); memset(h,-1,sizeof(h)); memset(f,-1,sizeof(f)); memset(w,-1,sizeof(w)); n=read(),m=read(),k=read(); size=n; for(int i=1;i<=n;i++) g[i]=read(); for(int i=1;i<=m;i++) { tp=read(),e1=read(); size++; ad2(find(tp),size),ad2(find(e1),size); ad2(size,find(tp)),ad2(size,find(e1)); hb(tp,size),hb(e1,size); } for(int i=1;i<=k;i++) yz[i]=read(),yy[i]=read(); e1=0; for(int i=k;i>=1;i--) ad1(yz[i],yy[i],i),ad1(yy[i],yz[i],i); } void tarjan(int x) { r[x]=1,zx[x]=x; for(int i=d[x];i!=-1;i=b[i].ne) { if(!r[b[i].v]) { tarjan(b[i].v); hb(x,b[i].v); zx[find(b[i].v)]=x; } hbb(x,b[i].v); } if(x<=n) for(int i=h[x];i!=-1;i=t[i].ne) if(r[t[i].v]&&fi(x)==fi(t[i].v)) fa[t[i].num]=zx[find(t[i].v)]; } int main() { init(); memset(f,-1,sizeof(f)); for(int i=n+m;i>=1;i--) if(!r[i]) tarjan(i); memset(b,0,sizeof(b)); memset(d,-1,sizeof(d)); for(int i=k;i>=1;i--) if(fa[i]) ad2(fa[i],i); for(int i=n+1;i<=n+m;i++) for(int j=d[i];j!=-1;j=b[j].ne) { tp=b[j].v; if(g[yy[tp]]==0||g[yz[tp]]==0) continue; if(g[yy[tp]]<=g[yz[tp]]) { ans+=g[yy[tp]]; g[yz[tp]]-=g[yy[tp]],g[yy[tp]]=0; } else { ans+=g[yz[tp]]; g[yy[tp]]-=g[yz[tp]],g[yz[tp]]=0; } } printf("%lld",ans<<1); return 0; }