浅谈无旋treap(fhq_treap)
一、简介
无旋Treap(fhq_treap),是一种不用旋转的treap,其代码复杂度不高,应用范围广(能代替普通treap和splay的所有功能),是一种极其强大的平衡树。
无旋Treap是一个叫做范浩强的大佬发明的(快%啊!)
在我们一起学习无旋Treap之前,本蒟蒻有几句活想说:
1.无旋Treap我个人认为是最容易理解的一种平衡树,而且编程复杂度不高,功能还那么强大。
我一开始学平衡树的时候是先从普通的带旋转的Treap开始学的。那种Treap,我现在都没搞懂什么左旋右旋究竟是怎么一回事。
我当时真的以为我这辈子都不指望学会平衡树了,直到djq大佬(快%!)告诉了我这种很美妙的数据结构。
我真就不明白了,为什么那么多人愿意去学普通的Treap而不学无旋Treap。无旋Treap更容易理解,最主要的是他与普通的Treap相比能可持久化!(尽管我不会)
所以,在大家学习完无旋Treap以后,本蒟蒻请诸君不妨推广一下无旋Treap,造福更多的Oier。
2.关于无旋Treap和其他平衡树的比较:(这个建议大家学完无旋Treap再来看,可能感触会更深刻一些)
与AVL相比:旋转操作真的很浪费时间,最坏情况下复杂度为O(log n),而且AVL树难写无比,不适合运用于算法竞赛。
与普通Treap相比:参见第一条
与splay相比:基本上可以代替splay的所有功能,但是在处理LCT问题上没有splay优秀
与rbt(红黑树相比):红黑树特别难写是众所周知的。
与sbt相比:sbt是我认为的最强大的平衡树。但是无旋Treap中的merge、split操作的应用的广泛(可持久化Treap维护Hash之类的)是sbt做不到的。
二、Treap
什么是Treap?
Treap=Tree+heap
相信大家都知道二叉堆吧。父节点的权值比子节点都要大(或小)
而Treap,则是在BST(二叉查找树)的基础上,添加二叉堆中的这个元素。
Treap与heap的区别是,heap是完全二叉树,而Treap不是。
下面的
三、核心操作
我在前面说过,普通的Treap最烦人的地方便是旋转。
而无旋Treap是如何做到无旋的呢?
关键就在两个操作:merge和split
1.split
split,顾名思义,就是把一个平衡树分成两棵树。
split有两种:一种是按照权值split,一种是按照size来split。
如果按照权值split,那么分出来两棵树的第一棵树上的每一个数的大小都小于(或小于等于,视具体情况而定)x;
如果按照size split,那么分出来两棵树的第一棵树恰好有x个节点。
我们可以结合具体代码讲解:
inline void split(int k,int& l,int& r,int x){//理解时,我们可以把l当做是答案的第一个,r当做答案的第二个,这个函数的意义是:将以k为根的树按照val分为以l为根的树和以r为根的树。注意引用的作用
if(!k){
l=r=;//分到底了,返回
return;
}
if(tree[k].val<x){//如果比它小
l=k;//那么x肯定在k的右子树里,先将k贴到第一个答案上
split(tree[l].r,tree[l].r,r,x);//把第一个答案的右子树按x分开,得到答案(这里自己理解一下,不难懂)
}else{//反之亦然
r=k;
split(tree[r].l,l,tree[r].l,x);
}
push_up(k);
}
而按size split的道理是一样的:
inline void split(int k,int& l,int& r,int x){
if(!k){
l=r=;
return;
}
if(tree[tree[k].l].size+<=x){
l=k;
split(tree[l].r,tree[l].r,r,x-tree[tree[k].l].size-);//注意这里有些变化
}else{
r=k;
split(tree[r].l,l,tree[r].l,x);
}
push_up(k);
}
这样就把一棵树分开了。
2.merge
merge就是把两颗原本分开的树合并在一起。
我们仍然结合具体代码讲解
inline void merge(int& k,int l,int r){//函数名的意义是:把以l为根的树和以r为根的树合并为以k为根的树
if(!l||!r){//合并到底,返回
k=l+r;
return;
}
if(tree[l].p>tree[r].p){//默认大根堆
k=l;//先把l挂到k上
merge(tree[k].r,tree[k].r,r);//注意要满足BST性质
}else{//反之亦然
k=r;
merge(tree[k].l,l,tree[k].l);
}
push_up(k);
}
有了merge和split,其他平衡树的基本操作就好做多了
三、其他操作
inline void insert(int val){
int x,y;
split(root,x,y,val-);
merge(x,x,New(val));
merge(root,x,y);
}
inline void Delete(int val){
int x,y,z;
split(root,x,y,val);
split(x,x,z,val-);
merge(z,tree[z].l,tree[z].r);
merge(x,x,z);
merge(root,x,y);
}
inline int rnk(int val){
int x,y;
split(root,x,y,val-);
int ans=tree[x].size+;
merge(root,x,y);
return ans;
}
inline int kth(int k){
int x=root;
while(true){
if(k==tree[tree[x].l].size+) return tree[x].val;
if(k<=tree[tree[x].l]) x=tree[x].l;
else k-=tree[tree[x].l].size+,x=tree[x].r;
}
}
inline int pre(int val){
int x,y;
split(root,x,y,val-);
int ans=tree[x].val;
merge(root,x,y);
return ans;
}
inline int suf(int val){
int x,y;
split(root,x,y,val);
int ans=kth(tree[x].size+);
merge(root,x,y);
return ans;
}
最后给大家放上一道模板题P3369 【模板】普通平衡树
如有不足请指正,谢谢
浅谈无旋treap(fhq_treap)的更多相关文章
- [转载]无旋treap:从单点到区间(例题 BZOJ1500&NOI2005 维护数列 )
转自ZZH大佬,原文:http://www.cnblogs.com/LadyLex/p/7182631.html 1500: [NOI2005]维修数列 Time Limit: 10 Sec Mem ...
- [转载]无旋treap:从好奇到入门(例题:bzoj3224 普通平衡树)
转载自ZZH大佬,原文:http://www.cnblogs.com/LadyLex/p/7182491.html 今天我们来学习一种新的数据结构:无旋treap.它和splay一样支持区间操作,和t ...
- [您有新的未分配科技点]无旋treap:从好奇到入门(例题:bzoj3224 普通平衡树)
今天我们来学习一种新的数据结构:无旋treap.它和splay一样支持区间操作,和treap一样简单易懂,同时还支持可持久化. 无旋treap的节点定义和treap一样,都要同时满足树性质和堆性质,我 ...
- Luogu 3369 / BZOJ 3224 - 普通平衡树 - [无旋Treap]
题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id=3224 https://www.luogu.org/problemnew/show/P3 ...
- 【算法学习】Fhq-Treap(无旋Treap)
Treap——大名鼎鼎的随机二叉查找树,以优异的性能和简单的实现在OIer们中广泛流传. 这篇blog介绍一种不需要旋转操作来维护的Treap,即无旋Treap,也称Fhq-Treap. 它的巧妙之处 ...
- 无旋treap的区间操作实现
最近真的不爽...一道维修数列就做了我1上午+下午1h+1晚上+晚上1h+上午2h... 一道不错的自虐题... 由于这一片主要讲思想,代码我放这里了 不会无旋treap的童鞋可以进这里 呵呵... ...
- 无旋treap的简单思想以及模板
因为学了treap,不想弃坑去学splay,终于理解了无旋treap... 好像普通treap没卵用...(再次大雾) 简单说一下思想免得以后忘记.普通treap因为带旋转操作似乎没卵用,而无旋tre ...
- [BZOJ3223]文艺平衡树 无旋Treap
3223: Tyvj 1729 文艺平衡树 Time Limit: 10 Sec Memory Limit: 128 MB Description 您需要写一种数据结构(可参考题目标题),来维护一个 ...
- [您有新的未分配科技点] 无旋treap:从单点到区间(例题 BZOJ1500&NOI2005 维护数列 )
1500: [NOI2005]维修数列 Time Limit: 10 Sec Memory Limit: 64 MB Description Input 输入的第1 行包含两个数N 和M(M ≤20 ...
随机推荐
- Bootstrap4 入门
http://www.runoob.com/bootstrap4/bootstrap4-navs.html 共五个部分 1 <!DOCTYPE html> <html lang=&q ...
- 外汇MT4编程手册
1.为了最大的方便用户,交易中断的一些变量可以从智能系统输入. AccountNumber-账号(同义词:AccNum) Ask –卖价(买方出价) Balance – 交易账户的余额值 Bars – ...
- javascript原型原型链 学习随笔
理解原型和原型链.需从构造函数.__proto__属性(IE11以下这个属性是undefined,请使用chrome调试).prototype属性入手. JS内置的好多函数,这些函数又被叫做构造函数. ...
- 时间都去哪儿了?开源一个统计iPhone上App运行时间和打开次数的小工具【iOS8已失效】
如今,大家每天都有大量时间花在手机上,但是,大家有没有想过自己的时间都花在哪些App上了呢?相信很多人都有这样的需求,不过iOS系统本身并不能显示每个App的运行时间和次数,因此,本人写了这样一个小工 ...
- Vue路由规则中定义参数
Vue使用routerLinke定义参数的时候 路由规则中不需要更改任何属性. 路由其实就是我们在html中定义的锚点,点击这个连接跳转一个锚点.vue中的路由也是这个原理, 前提是路由必须创建在实 ...
- MNIST手写数据集在运行中出现问题解决方案
今天在运行手写数据集的过程中,出现一个问题,代码没有问题,但是运行的时候一直报错,错误如下: urllib.error.URLError: <urlopen error [SSL: CERTIF ...
- JS在浏览器中输出各种三角形
直角三角形 <script type="text/javascript"> for(var i=1;i<=8;i++){ for(var j=1;j<=i; ...
- BZOJ 2321 星器
星器 思路: 势能分析法. 假设每颗星星的势能为\(x^2+y^2\) 那么对于一行的两颗星星\((i, j), (i, k), j < k\) 它转移到\((i, j+1), (i, k-1) ...
- 利用random模块做一个抢红包功能
我们都知道random模块是一个生成随机数的模块,用它来做抢红包的功能很合适. 抢红包,抢到的金额是随机的,但怎么让每个人抢到的随机金额公平合理呢 比如:我想发一个100元的红包让10个人抢,我可以把 ...
- certbot更新错误
自动更新老是提示这个错误. root@vultr:~/certbot# ./certbot-auto Upgrading certbot-auto 0.29.1 to 0.34.2... Couldn ...