[HDU 3712] Fiolki (带边权并查集+启发式合并)

[HDU 3712] Fiolki (带边权并查集+启发式合并)

题面

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

\(n,m \leq 2 \times 10^5,k \leq 5 \times 10^5\)

分析

注意到初始状态下第i个瓶子里有物质i,也就是说每种物质恰好只在一个瓶子里。那么混合的过程中，每种反应至多发生一次。对于一个反应\((a,b)\),因为开始反应前只有1个瓶子里有a,1个瓶子里有b.而当a,b相遇时，会一直进行直到完全反应。

那么，我们只要知道第i个反应发生的时间，然后按时间给反应排序(时间相同时按优先级排序).然后一个个反应按顺序模拟，更新反应物的质量和沉淀质量。

如何求某个反应(a,b)发生的时间呢?。我们把处在同一个烧杯里的物质看成一个联通块，(a,b)发生的时间就是a和b最早连通的时间。用并查集维护连通性，每个点x还要另外记录tim[x],表示x什么时间与父亲相连。答案就是u到v路径上的点tim的最大值。具体参考[BZOJ 4668]冷战(并查集+启发式合并)

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 200000
#define maxk 500000
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
struct disjoint_set{
    int fa[maxn+5];
    int tim[maxn+5];
    int sz[maxn+5];
    int find(int x){
        while(fa[x]!=x) x=fa[x];
        return x;
    }
    int get_deep(int x){
        int ans=0;
        while(fa[x]!=x){
            ans++;
            x=fa[x];
        }
        return ans;
    }
    void merge(int x,int y,int t){
        int fx=find(x);
        int fy=find(y);
        if(sz[fx]>sz[fy]) swap(fx,fy);
        fa[fx]=fy;
        tim[fx]=t;
        sz[fy]+=sz[fx];
    }
    int query(int x,int y){
        if(find(x)!=find(y)) return INF;
        int ans=0;
        int dx=get_deep(x),dy=get_deep(y);
        if(dx<dy){
            swap(x,y);
            swap(dx,dy);
        }
        while(dx>dy){
            ans=max(ans,tim[x]);
            x=fa[x];
            dx--;
        }
        if(x==y) return ans;
        while(x!=y){
            ans=max(ans,max(tim[x],tim[y]));
            x=fa[x];
            y=fa[y];
        }
        return ans;
    }
    void ini(int n){
        for(int i=1;i<=n;i++){
            fa[i]=i;
            sz[i]=1;
        }
    }
}S;
int n,m,k;
int g[maxn+5];
struct rec{
    int x;
    int y;
    int tim;
    int id;
    friend bool operator < (rec p,rec q){
        if(p.tim==q.tim) return p.id<q.id;
        else return p.tim<q.tim;
    }
}q[maxk+5];
int main(){
    int u,v;
    scanf("%d %d %d",&n,&m,&k);
    for(int i=1;i<=n;i++){
        scanf("%d",&g[i]);
    }
    S.ini(n);
    for(int i=1;i<=m;i++){
        scanf("%d %d",&u,&v);
        S.merge(u,v,i);
    }
    for(int i=1;i<=k;i++){
        scanf("%d %d",&q[i].x,&q[i].y);
        q[i].tim=S.query(q[i].x,q[i].y);
        q[i].id=i;
    }
    sort(q+1,q+1+k);
    ll ans=0;
    for(int i=1;i<=k;i++){
        if(q[i].tim==INF) continue;
        int x=q[i].x,y=q[i].y;
        int sum=min(g[x],g[y]);
        g[x]-=sum;
        g[y]-=sum;
        ans+=sum*2;
    }
    printf("%lld\n",ans);
}

原文地址：https://www.cnblogs.com/birchtree/p/11517009.html