动态树之Link-Cut Trees

    xiaoxiao2024-11-20  2

    关于动态树

    动态树( Dynamic Trees )问题, 即要求我们维护一个由若干棵子结点无序的有根树组成的森林。 要求这个数据结构支持对树的分割,合并,对某个点到它的根的路径的某些操作, 以及对某个点的子树进行的某些操作。(来源——QTREE解法的一些研究byYangZhe)

    解决动态树问题有很多种方法 ,这里介绍 Link - Cut Trees ,它不能够实现“对某个点的子树进行的某些操作”,但是对于大部分动态树问题来说还是够用了(其他动态树好像比较冷门)。

    动态树的主体思想和树链剖分很相似(树链剖分戳这里),非常频繁地运用到了 Splay Tree (Splay模板戳这里)


    无关内容: 还有一种动态树叫做Euler-Tour Trees,能够实现对子树的某些操作,然而我并没有找到有任何关于它的中文文献。。。 所以等到我哪天英文达到了那个水平再去学吧。。。。


    介绍

    Link - Cut Trees 是由 Sleator Tarjan 发明的,这里膜拜一下 Tarjan (Orz),LCT,LCA,LCP的发明者中都有他的名字。。。(还有很多很多算法。。。)

    定义一些神奇的称呼:

    access(x) :访问节点 x PreferredChild:如果结点 v 的子树中, 最后被访问的结点在子树w中, 这里 w v的儿子, 那么就称 w v PreferredChild PreferredEdge :每个点到它的 PreferredChild 的边称作 PreferredEdge PreferredPath :由 PreferredEdge 连接成的不可再延伸的路径称为 PreferredPath

    其中 access(x) Link - Cut Trees 最基本的操作,没有之一。就像 Splay 操作对于 Splay Tree 一样。

    容易得出整棵树就被划分成了若干条 PreferredPath 。对每条 PreferredPath , 用这条路上的点的深度作为关键字, 用一棵平衡树来维护它(一般使用Splay,理论上Treap也可以,可是我从没有看见有人这样做过)。然后这棵平衡树就被叫做 AuxiliaryTree 。我们把 AuciliaryTree 中深度最小的节点的父亲节点称为 PathParent

    Link - Cut Trees 就是将要维护的森林中的每棵树 T 表示为若干个 AuxiliaryTree , 并通过 PathParent 将这些 AuxiliaryTree 连接起来的数据结构

    access(x) 操作

    一旦我们调用 access(x) ,那么从点 x 到根结点的路径就成为一条新的 PreferredPath。 如果路径上的某个节点 u 不是它的父亲v PreferredChild , 那么我们要将 v PreferredChild变为 u , 原本包含v PreferredPath/AuxiliaryTree 将不再包含节点 v 及其之上的部分。

    具体操作: 首先,我们将节点v到它的 PreferredChild 的边删除,将节点 v Splay到它所属的 AuxiliaryTree 的根,然后将 v 与它的右子树断开,使它的右子树变为一棵新的AuxiliaryTree,并将它的 PathParent 设置为节点 v 。 然后,如果点v所属的 PreferredPath 并不包含根结点,设它的 PathParent u , 那么需要将u旋转到 u 所属的AuxiliaryTree的根,并用点 v 所属的AuxiliaryTree替换到点 u 所属的AuxiliaryTree中点 u 的右子树, 再将原来点u所属的 AuxiliaryTree 中点 u 的右子树的PathParent设置为 u

    重复以上操作,直到到达包含根结点的PreferredPath。 给出几张图片,方便大家更好地理解 access(x) 操作。

    有了 access(x) 剩下的操作都变得十分地简单。

    find _ root(x) 操作

    找到节点 x 所在树的根节点。

    具体操作: 先access(x),然后将 x Splay到所在 AuxiliaryTree 的根节点。找到这棵 AuxiliaryTree 最左端的点即可。

    cut(x) 操作

    断开 x 与其父亲节点的边。

    具体操作: 先access(x),然后将 x Splay到所在 AuxiliaryTree 的根节点。断开 x 和父亲节点的边。

    join(v,w)操作

    v 成为w的新的儿子。其中 v 是一棵树的根结点,并且v w 是不同的两棵树中的结点。

    具体操作: 先访问v,然后修改 v 所属的AuxiliaryTree PathParent w ,然后再次访问v

    具体代码

    struct LCT { int fa[MAXN], ch[MAXN][2], data[MAXN], pathparent[MAXN], pn; void init(int n) { pn = n; memset(fa, 0, sizeof(fa)); memset(ch, 0, sizeof(ch)); } void Rotate(int p, bool t) { int f = fa[p]; fa[ch[f][t^1] = ch[p][t]] = f; fa[ch[fa[f]][ch[fa[f]][1]==f] = p] = fa[f]; ch[fa[f] = p][t] = f; } void splay(int x) { int p; while(fa[x]) { p = fa[x]; if(!fa[p]) { Rotate(x, x==ch[p][0]); break; } bool f = x==ch[p][0], f1 = p==ch[fa[p]][0], f2 = p==ch[fa[p]][1]; Rotate(f?f1?p:x:f2?p:x, f); Rotate(x, f1); } } void access(int v) { int u = v; v = 0; while(u) { splay(u); pathparent[ch[u][1]] = u; ch[u][1] = v; pathparent[v] = 0; v = u; u = pathparent[u]; } } int find_root(int v) { access(v); splay(v); while(ch[v][0]) v = ch[c][0]; splay(v); return v; } void cut(int v) { access(v); ch[v][0] = 0; } void join(int v, int w) { access(v); pathparent[v] = w; access(v); } };

    我们发现,当一个节点u位于所在的 AuxiliaryTree 的根节点时,它的父亲为0,我们可以做相应的优化,不再需要 PathParent 数组,根节点的父亲设置为 PathParent ,然后增加一个 bool 型数组,表示一个节点是否为根节点,这样可以节省一部分空间。

    struct LCT { bool root[MAXN]; int fa[MAXN], ch[MAXN][2]; void init() { memset(fa, 0, sizeof(fa)); memset(ch, 0, sizeof(ch)); memset(root, true, sizeof(root)); } void Rotate(int p, bool t) { int f = fa[p]; fa[ch[f][t^1] = ch[p][t]] = f; if(!root[f]) ch[fa[f]][ch[fa[f]][1]==f] = p; else root[p] = !(root[f] = false); fa[p] = fa[f]; ch[fa[f] = p][t] = f; } void splay(int x) { while(!root[x]) { int p = fa[x]; if(root[p]) { Rotate(x, x==ch[p][0]); break; } bool f = x==ch[p][0], f1 = p==ch[fa[p]][0], f2 = p==ch[fa[p]][1]; Rotate(f?f1?p:x:f2?p:x, f); Rotate(x, f1); } } void access(int u) { int v = 0; while(u) { splay(u); root[ch[u][1]] = !(root[v] = false); ch[u][1] = v; u = fa[v = u]; } } int find_root(int v) { access(v); splay(v); while(ch[v][0]) v = ch[v][0]; splay(v); return v; } void cut(int v) { access(v); splay(v); root[ch[v][0]] = true; ch[v][0] = fa[ch[v][0]] = 0; } void join(int v, int w) { access(v); fa[v] = w; access(v); } };

    以下是裸题: HNOI2010 bounce 弹飞绵羊 题目戳这里 题解戳这里 SDOI2008 cave 洞穴勘测 题目戳这里 HDU4010 Query on The Trees 题目戳这里 SPOJ00913 Query on a tree II 题目戳这里

    转载请注明原文地址: https://ju.6miu.com/read-1293838.html
    最新回复(0)