算法提高——高级数据结构 (带权并查集)

---

一、例题

食物链
动物王国中有三类动物 A,B,C，这三类动物的食物链构成了有趣的环形。

A 吃 B，B 吃 C，C 吃 A。

现有 N 个动物，以 1∼N 编号。

每个动物都是 A,B,C 中的一种，但是我们并不知道它到底是哪一种。

有人用两种说法对这 N 个动物所构成的食物链关系进行描述：

第一种说法是 1 X Y，表示 X 和 Y 是同类。

第二种说法是 2 X Y，表示 X 吃 Y。

此人对 N 个动物，用上述两种说法，一句接一句地说出 K 句话，这 K 句话有的是真的，有的是假的。

当一句话满足下列三条之一时，这句话就是假话，否则就是真话。

当前的话与前面的某些真的话冲突，就是假话；
当前的话中 X 或 Y 比 N 大，就是假话；
当前的话表示 X 吃 X，就是假话。
你的任务是根据给定的 N 和 K 句话，输出假话的总数。

输入格式
第一行是两个整数 N 和 K，以一个空格分隔。

以下 K 行每行是三个正整数 D，X，Y，两数之间用一个空格隔开，其中 D 表示说法的种类。

若 D=1，则表示 X 和 Y 是同类。

若 D=2，则表示 X 吃 Y。

输出格式
只有一个整数，表示假话的数目。

数据范围
1≤N≤50000,
0≤K≤100000
输入样例：
100 7
1 101 1
2 1 2
2 2 3
2 3 3
1 1 3
2 3 1
1 5 5
输出样例：
3

题目来源：算法竞赛进阶指南/ NOI2001 / Kuangbin专题 /
题目链接：食物链

二、思路分析：

根据分析题目内容，理解题意大概是在只有三只动物的环形食物链，来分析n只生物分别属于那种动物。食物链大概是这样（A->B->C->A），我们来分析题目，这里我们来分析生物属于那种动物可以用并查集来做，但是每种动物之间的关系只有用带权并查集来维护。（如果第一次接触这里很难想到，要不然本题也不会成为经典）我们的权值可以来如下定义：1表示能吃根节点，2表示可以被根结点吃，0则表示和根节点是同类。如果画在图上这个关系也符合三种动物环环相扣的状态，但是我们需要维护的动物远远不知三只，这时我们可以将权值来对3取mod，这样一来我们就解决了动物和数据存储之间的关系。

三、理解带权并查集：

观察代码，首先我们要记住两点，第一并查集其实是从底往上查找的过程，第二并查集的权值维护只能自上而下来维护，要理解这点，我们来细想一下，如果我们是从下往上维护的，那当我们要更新x结点权值时，那我们就需要他的父节点已经被更新过的状态，因此如果从下往上查找的话显然没办法做到 d[x] += d[p[x]]，因此我们只有自上而下来更新权值。

好了，既然已经了解了运作顺序，那我们就可以来看下代码，每当我们遇到一个新结点x，并在下面时，在调用一次Find函数后，就会开始路径压缩，结点x会先从底往上找，这个过程不更新任何数据，指到找到根结点(p[x]==x)，并返回根节点,此时开始回溯，来逐步压缩路径并更新权值，更新权值时，我们其实直接加上一个结点权值即可，因为上一个结点必然已经被更新过了，所以我们直接d[x]+=d[p[x]]即可，但是注意，此时我们是需要让这个新结点也指向根结点，因此在回溯过程我们其实需要用一个变量u来保存根结点(u=Find(p[x]))，因此我们才能在后边p[x]=u,而此时的x结点父节点也是根节点，因此直接返回p[x],让下一个结点也能取到u即可。

四、如何合并？

因为结点是带权的，所以合并可不能含糊合并。假设两个颗树想要合并，那我们要考虑好fx结点到fy结点的路径值,在这里我们设为d。当D1时，如果两个结点要合并，因为x结点和y结点为同类，根据我们的定义，同类结点的权值为0，因此我们不难推出(d[x]+d-d[y])%3 == 0, 因此fx的权值就为d=d[y]-d[x]，同理我们可以推出D2的情况。

最后附上代码

#include <iostream>

using namespace std;

const int N = 50010;

int p[N],d[N];

int Find(int x)
{
    if(p[x] != x)
    {
        int u = Find(p[x]);
        d[x] += d[p[x]];
        p[x] = u;
    }
    return p[x];
}
int main()
{
    int n,k,res = 0;
    cin >> n >> k;
    for (int i = 1; i <= n; i ++ ) p[i] = i;
    while(k -- )
    {
        int a,x,y;
        cin >> a >> x >> y;
        if(x > n || y > n) res ++;
        else
        {
            if(a == 1)
            {
                int fx = Find(x),fy = Find(y);
                if(fx == fy && (d[x] - d[y]) % 3) res ++;
                else if(fx != fy)
                {
                    p[fx] = fy;
                    d[fx] = d[y] - d[x];
                }
            }
            else if(a == 2)
            {
                int fx = Find(x),fy = Find(y);
                if(fx == fy && (d[x] - d[y] - 1) % 3) res ++;
                else if(fx != fy)
                {
                    p[fx] = fy;
                    d[fx] = d[y] + 1 - d[x];
                }
            }
        }
    }
    
    cout<<res<<endl;
    
    return 0;
}

总结