Treap——大名鼎鼎的随机二叉查找树,以优异的性能和简单的实现在OIer们中广泛流传。

这篇blog介绍一种不需要旋转操作来维护的Treap,即无旋Treap,也称Fhq-Treap。

它的巧妙之处在于只需要分离和合并两种基本操作,就能实现任意的平衡树常用修改操作。

而不需要旋转的特性也使编写代码时不需要考虑多种情况和复杂的父亲儿子关系的更新,同时降低了时间复杂度。

此外,它还可以方便地支持可持久化,实在是功能强大。

接下来系统性地介绍一下无旋Treap的原理和实现,最后讲解一下应用和例题。

一、Treap是啥?

Treap是一棵二叉查找树,满足中序遍历始终是有序的。

Treap的每个节点除了保存信息外,还要保存一个值\(pri\)。

这个\(pri\)存储的是这个节点的优先级。

它是干嘛用的呢?

实际上,Treap除了是一个二叉查找树之外,它在另一个意义下还是一个堆。

每个节点的子节点,一定要满足子节点的\(pri\)比父节点的\(pri\)小。\(pri\)越大的点越往上放。

这有什么意义呢?为什么要这样定义?

事实上,\(pri\)值是程序随机指派的,每个点的\(pri\)值是与这个点的权值无关的,是随机的。

这样随机化就可以保证Treap的深度是\(O(log\;n)\),这就是随机化的力量。

而且可以保证这样建树一定不会出现矛盾,现在模拟一下建树的过程:

先把所有节点按照中序遍历排好,然后找到其中\(pri\)最大的,把它作为整个Treap的根,左边的节点形成左子树,右边的节点形成右子树。再对左右子树递归处理。

这样最终就能建好一棵Treap。

二、Treap的基本操作

讲完了Treap的定义,来看看Treap的两个基本操作:

分离(Split)和合并(Merge)。

分离:指的是将一棵Treap按照中序遍历的顺序,分割成左右两半,满足左右两半组成的Treap的所有值都不变。

合并:指的是将两棵Treap(一般是从原先的TreapSplit出来的)合并在一起,按照中序遍历的顺序,并且所有节点的值都不变。

Split操作比较简单,先讲讲如何实现:当然,Split之前要先指定一个值k,表示Split出这个Treap的中序遍历中的前k个数作为第一棵Split出的Treap。

从这个Treap的根开始,看它的左子树的大小是否大于等于k,如果是,那么说明右子树和根都在第二棵中,继续递归到左子树中。

如果不是,那么说明左子树和根都在第一棵Treap中,继续递归到右子树中,而且k要减去左子树的大小加一。

代码:

void Split(int rt,int k,int&rt1,int&rt2){
if(!rt) {rt1=rt2=0; return;}
if(k<=siz[ls[rt]]){
Split(ls[rt],k,rt1,rt2);
ls[rt]=rt2;
combine(rt);
rt2=rt;
}
else{
Split(rs[rt],k-siz[ls[rt]]-1,rt1,rt2);
rs[rt]=rt1;
combine(rt);
rt1=rt;
}
}

其中,rt是根,k是分出的第一棵子树的大小,rt1和rt2用来返回。

combine函数用来维护节点的大小。

接下来看Merge操作:

有两棵Treap,假设要把第二棵接到第一棵后面,那么应该怎么合并呢?

考虑两个根节点的\(pri\)值,因为第一棵在第二棵前面,所以要不然rt1(第一棵的根)在rt2(第二棵的根)的左子树,要不然rt2在rt1的右子树。

但是因为有了\(pri\)的影响,所以只能rt1和rt2中\(pri\)较大的那个作为根。

如果rt1为根,那么有rt1的右子树和rt2合并作为rt1的现在的右子树。

如果rt2为根,那么有rt2的左子树和rt1合并作为rt2的现在的左子树。

两种情况都递归进子树中即可。

代码如下:

int Merge(int rt1,int rt2){
if(!rt1) return rt2;
if(!rt2) return rt1;
if(pri[rt1]<pri[rt2]){
rs[rt1]=Merge(rs[rt1],rt2);
combine(rt1);
return rt1;
}
else{
ls[rt2]=Merge(rt1,ls[rt2]);
combine(rt2);
return rt2;
}
}

这就是两种Treap的基本操作,实现时要注意当树为空时,做特殊处理。

三、主要操作

看完了这两种操作,你肯定会问:这有个鬼用?

没法插入删除查询前驱后继排名位置的平衡树,只能分分合合的平衡树,我才不要!

但是,以上的所有操作,其实都能通过这两种基本操作实现:

查询小于等于val的数的个数:考察根节点的值和val,如果val小于根节点,递归进左子树,否则递归进右子树。

这个查询操作命名为Rank,代码:

int Rank(int rt,int v){
if(!rt) return 0;
if(v<val[rt]) return Rank(ls[rt],v);
else return siz[ls[rt]]+Rank(rs[rt],v)+1;
}

插入val:先查询Rank(val),然后按照rank把整个TreapSplit成两个,把val做成一个新节点,Merge到里面即可。

void Insert(int v){
val[++cnt]=v, pri[cnt]=ran(), siz[cnt]=1;
int rank=Rank(Root,v);
int rt1,rt2;
Split(Root,rank,rt1,rt2);
Root=Merge(Merge(rt1,cnt),rt2);
}

删除val:先查询Rank(val),然后按照rank把整个TreapSplit成三个,删除需要的点,最后Merge剩下两个。

void Delete(int v){
int rank=Rank(Root,v);
int rt1,rt2,rt3,tmp;
Split(Root,rank,rt1,rt2);
Split(rt1,rank-1,rt3,tmp);
Root=Merge(rt3,rt2);
// Memory Recycle?
}

查询第K个值:把整个TreapSplit成三个,输出需要的值,最后合并起来。

int Kth(int k){
int rt1,rt2,rt3,c;
Split(Root,k,rt1,rt2);
Split(rt1,k-1,rt3,c);
Root=Merge(rt3,Merge(c,rt2));
return val[c];
}

前驱:Kth(Rank(val-1))。

后继:Kth(Rank(val)+1)。

还有很多很多操作,供大家脑补。

四、区间操作

无旋Treap和旋转Treap的更重要区别是,无旋Treap可以很方便地支持区间的操作。

如何支持?你已经看到了,Split操作分离出的就是一个个区间啊!把一整棵Treap分离出来一段区间,在上面尽情地修改吧。不过要记得像线段树写好区间pushdown哦!

比如区间翻转,就是左右子树调换,并且打上标记。区间加减乘除更不用说。

记得在Split,Merge和Rank三个函数内部加上pushdown!

五、实战应用

有很多平衡树的题目,都能用Treap解决。

推荐几道题:

Luogu P3369 测试你的Treap普通操作的熟练程度。

Luogu P3391 区间操作的应用。

Luoge P3165 区间操作和其他技巧。

【算法学习】Fhq-Treap(无旋Treap)的更多相关文章

  1. fhq Treap(无旋Treap)

    先吹一波fhq dalao,竟然和我一个姓,我真是给他丢脸. 昨天treap就搞了一下午,感觉自己弱爆了.然后今天上午又看了一个上午的无旋treap再次懵逼,我太弱了,orzorz. 所以写个博客防止 ...

  2. 浅谈无旋treap(fhq_treap)

    一.简介 无旋Treap(fhq_treap),是一种不用旋转的treap,其代码复杂度不高,应用范围广(能代替普通treap和splay的所有功能),是一种极其强大的平衡树. 无旋Treap是一个叫 ...

  3. [转载]无旋treap:从单点到区间(例题 BZOJ1500&NOI2005 维护数列 )

    转自ZZH大佬,原文:http://www.cnblogs.com/LadyLex/p/7182631.html 1500: [NOI2005]维修数列 Time Limit: 10 Sec  Mem ...

  4. [转载]无旋treap:从好奇到入门(例题:bzoj3224 普通平衡树)

    转载自ZZH大佬,原文:http://www.cnblogs.com/LadyLex/p/7182491.html 今天我们来学习一种新的数据结构:无旋treap.它和splay一样支持区间操作,和t ...

  5. [您有新的未分配科技点]无旋treap:从好奇到入门(例题:bzoj3224 普通平衡树)

    今天我们来学习一种新的数据结构:无旋treap.它和splay一样支持区间操作,和treap一样简单易懂,同时还支持可持久化. 无旋treap的节点定义和treap一样,都要同时满足树性质和堆性质,我 ...

  6. Luogu 3369 / BZOJ 3224 - 普通平衡树 - [无旋Treap]

    题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id=3224 https://www.luogu.org/problemnew/show/P3 ...

  7. [您有新的未分配科技点] 无旋treap:从单点到区间(例题 BZOJ1500&NOI2005 维护数列 )

    1500: [NOI2005]维修数列 Time Limit: 10 Sec  Memory Limit: 64 MB Description Input 输入的第1 行包含两个数N 和M(M ≤20 ...

  8. HNOI2012 永无乡 无旋Treap

    题目描述 永无乡包含 nnn 座岛,编号从 111 到 nnn ,每座岛都有自己的独一无二的重要度,按照重要度可以将这 nnn 座岛排名,名次用 111 到 nnn 来表示.某些岛之间由巨大的桥连接, ...

  9. [Bzoj3224][Tyvj1728] 普通平衡树(splay/无旋Treap)

    题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=3224 平衡树入门题,学习学习. splay(学习yyb巨佬) #include<b ...

随机推荐

  1. 多进程编程之system()函数

    1.system函数: 使用函数system,在程序中执行一个shell命令字符串很方便.它是一个和操作系统紧密相关的函数,用户可以使用它在自己的程序中调用系统提供的各种命令,执行系统的命令行,其实也 ...

  2. c语言宏定义#define

    1. 利用define来定义 数值宏常量 #define 宏定义是个演技非常高超的替身演员,但也会经常耍大牌的,所以我们用它要慎之又慎.它可以出现在代码的任何地方,从本行宏定义开始,以后的代码就就都认 ...

  3. Java并发编程:线程池

    一.为什么使用线程池 使用线程的时候直接就去创建一个线程,这样实现起来非常简便,但是就会有一个问题: 如果并发的线程数量很多,并且每个线程都是执行一个时间很短的任务就结束了,这样频繁创建线程就会大大降 ...

  4. http的无状态无连接

    搞爬虫的核心:http协议. 在理解http中的无状态和无连接时,有一些困惑,下文可以解决. 转自:http://www.cnblogs.com/bellkosmos/p/5237146.html h ...

  5. [提升性选讲] 树形DP进阶:一类非线性的树形DP问题(例题 BZOJ4403 BZOJ3167)

    转载请注明原文地址:http://www.cnblogs.com/LadyLex/p/7337179.html 树形DP是一种在树上进行的DP相对比较难的DP题型.由于状态的定义多种多样,因此解法也五 ...

  6. ORA-01410: 无效的 ROWID

    视图查询单表是有这个东西的,但是视图的SQL涉及多表关联,就没这个rowid了,要么自己写个,要么不用这个ROWID做啥动作

  7. MT【130】Heilbronn问题

    (清华THUSSAT,多选题) 平面上 4 个不同点 \(P_1,P_2,P_3,P_4\),在每两个点之间连接线段得到 6 条线段. 记 \[L=\max_{1\leq i<j\leq 4}| ...

  8. 51nod 1218 最长递增子序列 | 思维题

    51nod 1218 最长递增子序列 题面 给出一个序列,求哪些元素可能在某条最长上升子序列中,哪些元素一定在所有最长上升子序列中. 题解 YJY大嫂教导我们,如果以一个元素结尾的LIS长度 + 以它 ...

  9. Alpha 冲刺 —— 十分之四

    队名 火箭少男100 组长博客 林燊大哥 作业博客 Alpha 冲鸭鸭鸭鸭! 成员冲刺阶段情况 林燊(组长) 过去两天完成了哪些任务 协调各成员之间的工作 协助前后端接口的开发 测试项目运行的服务器环 ...

  10. bzoj 4521: [Cqoi2016]手机号码

    感觉get到了一种数位dp的新姿势,加一位表示当前要填的数有没有限制(感觉以前的写法都太蠢了). 这么写有两个地方要注意: 1.每dp到一位时需要f[i][初始状态]++,相当于这位前都是前导零(这道 ...