鲜花

一些鲜花放在前面,平衡树学了很久,但是每学一遍都忘,原因就在于我只能 70% 理解 + 30% 背板子,所以每次都忘。这次我采取了截然不同的策略,自己按照自己的理解打一遍,大获成功(?),大概打 20 min,调 10 min 结束,然后写下了这篇文章。

虽然但是,感觉 Treap 还是很强的,代码好写好调,而且可以解决很多问题(下面将会提到),好像说是常数大一点?无伤大雅吧……

Treap

首先,从 Treap 的定义开始。Treap 实际上是一种笛卡尔树(笛卡尔树可以看 这篇文章,每个点有两个信息 \((v_u,p_u)\),分别表示点 \(u\) 的权值与优先度。\(v_u\) 形成一个二叉搜索树(左儿子的值 \(\leq\) 自己的值 \(\leq\) 右儿子的值),\(p_u\) 形成一个大根堆(自己的值 \(\geq\) 左右儿子的值)。

你可以利用 Treap 做很多有用的操作,但是 Treap 有一个缺点,就是如果被卡成链怎么办?

FHQ Treap

这时候就需要用到 FHQ Treap,它基于 Treap 的基础上对每个点的 \(p\) 值都进行了随机化操作,这样就可以达到平衡的目的了。为什么?因为随机出来不平衡的概率很小,大多随机数都是无规律的,这样通过大根堆的性质就可以使它达到平衡,这样树高期望就是 \(O(\log n)\) 层的了,并且无法卡掉,因为随机数是程序生成的,并非数据所决定。

系统随机函数若以 srand(time(0)) 为随机种子,效果不佳,推荐使用 mt19937 rnd(233),生成随机数更加均匀。

FHQ Treap 又叫无旋 Treap,它只需要两个核心操作 merge()split() 即可完成所有的复杂操作,就像玩拼图一样把一棵树拆开再拼起来一样。下面就来具体介绍一下这两个函数到底是如何实现的。

关于 upd 函数的说明

(这个可以学完下面的再看)upd 函数,即刷新一个节点大小的函数。在 split() 中,在分裂完子树后最后应该刷新,不然可能大小信息已经被更改而不知道。同理,在 merge() 中,合并子树后也应当刷新,调不出来一定要检查一下是否因忘记刷新而错误。

split 函数

split 函数实现的功能是把一棵树按照权值大小拆分成两棵树,具体来说,split(int cur,int k,int &x,int &y) 表示将以 cur 为根的子树拆分成两棵树 \(x\) 和 \(y\),其中 \(x\) 里的权值都 \(\leq k\),\(y\) 里的权值都 \(>k\)。

怎么写呢?首先,为了方便返回,直接将要拆分成的两棵子树引用定义在函数参数中。先特判一下,如果要拆分的树为空,则拆分出来的两棵树也为空。不然就判断树根是否 \(\leq k\),如果是,则树左子树都 \(\leq k\),即属于 \(x\),那么只要分裂右子树即可,反之亦然。

代码
void split(int cur,int k,int &x,int &y)
{
if(!cur)
{
x=y=0;
return;
}
if(tree[cur].val<=k)
{
x=cur;
split(tree[x].rs,k,tree[x].rs,y);
}
else
{
y=cur;
split(tree[y].ls,k,x,tree[y].ls);
}
upd(cur);
}

merge 函数

merge 函数就是把两个树 \(x,y\) 合并,保证所有 \(x\) 中的 \(v\) 都 \(\leq\) 所有 \(y\) 中的 \(v\),具体来说,就是 int merge(int x,int y)

写法就是判断 \(x,y\) 中是否有空树,有的话直接返回另一个即可。不然比较两者根的有先值,如果第一个大,根据大根堆性质,显然应该合并 \(x\) 的右子树和 \(y\),反之亦然。

代码
int merge(int x,int y)
{
if(!x||!y)
return x+y;
if(tree[x].key>=tree[y].key)
{
tree[x].rs=merge(tree[x].rs,y);
upd(x);
return x;
}
else
{
tree[y].ls=merge(x,tree[y].ls);
upd(y);
return y;
}
}

其他操作的实现

下面来看看 FHQ Treap 能实现哪些操作吧,请看模板题 普通平衡树。(下面全部内容为笔者自己发挥,若有更好做法请在评论区留言或私信笔者)

插入 \(x\):将根以 \(x\) 值分裂,在中间建一个新的点合并回去即可,比较容易。

删除 \(x\):将根分别以 \(x-1,x\) 分裂,中间一个树就是权值为 \(x\) 的树,把这个树替换为它左右子树合并的结果即可(因为只要删一个,相当于把根节点给删了。

查询 \(x\) 的排名:这个非常容易,直接按 \(x-1\) 分裂并输出第一个子树大小 +1 即可。

查询排名为 \(x\) 的数:这个可以从根节点开始,不停地判断应该往左子树走还是右子树走,判断方式就是看当前点地左子树大小 +1 和 \(x\) 的大小关系,若相等则直接退出。

**查询 \(x\) 的前驱:$$笔者做法比较暴力,直接将数按 \(x-1\) 分裂,第一棵树的最后一个就是,最后一个求可以查询排名为数大小的数,用之前的函数可以求出。

**查询 \(x\) 的后继:$$同前驱做法,按 \(x\) 分裂,第二棵树的第一个就是,同样查询排名为 1 的树即可,使用之前的函数。

具体的还是看代码吧。

代码
#include <bits/stdc++.h>
#define TIME 1e3*clock()/CLOCKS_PER_SEC
using namespace std;
// stay organized
mt19937 rnd(233);
const int maxn=1e5+10;
struct Node
{
int ls,rs;
int siz;
int val,key;
}tree[maxn];
int rt=0,tot=0;
int newnode(int val)
{
tot++;
tree[tot].ls=tree[tot].rs=0;
tree[tot].siz=1;
tree[tot].val=val;
tree[tot].key=rnd();
return tot;
}
void upd(int x)
{
tree[x].siz=tree[tree[x].ls].siz+1+tree[tree[x].rs].siz;
}
void split(int cur,int k,int &x,int &y)
{
if(!cur)
{
x=y=0;
return;
}
if(tree[cur].val<=k)
{
x=cur;
split(tree[x].rs,k,tree[x].rs,y);
}
else
{
y=cur;
split(tree[y].ls,k,x,tree[y].ls);
}
upd(cur);
}
int merge(int x,int y)
{
if(!x||!y)
return x+y;
if(tree[x].key>=tree[y].key)
{
tree[x].rs=merge(tree[x].rs,y);
upd(x);
return x;
}
else
{
tree[y].ls=merge(x,tree[y].ls);
upd(y);
return y;
}
}
int x,y,z;
void ins(int val)
{
split(rt,val,x,y);
rt=merge(merge(x,newnode(val)),y);
}
void del(int val)
{
split(rt,val-1,x,y);
split(y,val,y,z);
y=merge(tree[y].ls,tree[y].rs);
rt=merge(merge(x,y),z);
}
int rnk(int val)
{
split(rt,val-1,x,y);
int ans=tree[x].siz+1;
rt=merge(x,y);
return ans;
}
int kth(int now,int k)
{
while(tree[tree[now].ls].siz+1!=k)
{
if(tree[tree[now].ls].siz+1>k)
now=tree[now].ls;
else
{
k-=tree[tree[now].ls].siz+1;
now=tree[now].rs;
}
}
return tree[now].val;
}
int pre(int val)
{
split(rt,val-1,x,y);
int ans=kth(x,tree[x].siz);
rt=merge(x,y);
return ans;
}
int suf(int val)
{
split(rt,val,x,y);
int ans=kth(y,1);
rt=merge(x,y);
return ans;
}
int main()
{
ios::sync_with_stdio(false);
cin.tie(0);
cout.tie(0);
int t;
cin>>t;
int opt,x;
while(t--)
{
cin>>opt>>x;
if(opt==1)
ins(x);
else if(opt==2)
del(x);
else if(opt==3)
cout<<rnk(x)<<'\n';
else if(opt==4)
cout<<kth(rt,x)<<'\n';
else if(opt==5)
cout<<pre(x)<<'\n';
else cout<<suf(x)<<'\n';
}
return 0;
// you should actually read the stuff at the bottom
} /* stuff you should look for
* int overflow, array bounds
* clear the arrays?
* special cases (n=1?),
* WRITE STUFF DOWN,
* DON'T GET STUCK ON ONE APPROACH
*/

完结撒花 owo~

FHQ Treap 详解的更多相关文章

  1. 【数据结构】FHQ Treap详解

    FHQ Treap是什么? FHQ Treap,又名无旋Treap,是一种不需要旋转的平衡树,是范浩强基于Treap发明的.FHQ Treap具有代码短,易理解,速度快的优点.(当然跟红黑树比一下就是 ...

  2. Treap详解

    今天一天怼了平衡树.深深地被她的魅力折服了.我算是领略到了高级数据结构的美妙.oi太神奇了. 今天初识平衡树,选择了Treap. Treap又叫树堆,是一个二叉搜索树.我们知道,它的节点插入是随机的, ...

  3. 非旋 treap 结构体数组版(无指针)详解,有图有真相

    非旋  $treap$ (FHQ treap)的简单入门 前置技能 建议在掌握普通 treap 以及 左偏堆(也就是可并堆)食用本blog 原理 以随机数维护平衡,使树高期望为logn级别, FHQ  ...

  4. Linq之旅:Linq入门详解(Linq to Objects)

    示例代码下载:Linq之旅:Linq入门详解(Linq to Objects) 本博文详细介绍 .NET 3.5 中引入的重要功能:Language Integrated Query(LINQ,语言集 ...

  5. 架构设计:远程调用服务架构设计及zookeeper技术详解(下篇)

    一.下篇开头的废话 终于开写下篇了,这也是我写远程调用框架的第三篇文章,前两篇都被博客园作为[编辑推荐]的文章,很兴奋哦,嘿嘿~~~~,本人是个很臭美的人,一定得要截图为证: 今天是2014年的第一天 ...

  6. EntityFramework Core 1.1 Add、Attach、Update、Remove方法如何高效使用详解

    前言 我比较喜欢安静,大概和我喜欢研究和琢磨技术原因相关吧,刚好到了元旦节,这几天可以好好学习下EF Core,同时在项目当中用到EF Core,借此机会给予比较深入的理解,这里我们只讲解和EF 6. ...

  7. Java 字符串格式化详解

    Java 字符串格式化详解 版权声明:本文为博主原创文章,未经博主允许不得转载. 微博:厉圣杰 文中如有纰漏,欢迎大家留言指出. 在 Java 的 String 类中,可以使用 format() 方法 ...

  8. Android Notification 详解(一)——基本操作

    Android Notification 详解(一)--基本操作 版权声明:本文为博主原创文章,未经博主允许不得转载. 微博:厉圣杰 源码:AndroidDemo/Notification 文中如有纰 ...

  9. Android Notification 详解——基本操作

    Android Notification 详解 版权声明:本文为博主原创文章,未经博主允许不得转载. 前几天项目中有用到 Android 通知相关的内容,索性把 Android Notificatio ...

随机推荐

  1. spring使用junit单元测试

    <dependency> <groupId>org.springframework</groupId> <artifactId>spring-test& ...

  2. 熔断器-Hystrix。。。之降级方法

    与Feign的Fallback降级方法不同,这个hystrix降级方法是写在被调用方的 需要依赖: <dependency> <groupId>org.springframew ...

  3. LOJ6062「2017 山东一轮集训 Day2」Pair(Hall定理,线段树)

    题面 给出一个长度为 n n n 的数列 { a i } \{a_i\} {ai​} 和一个长度为 m m m 的数列 { b i } \{b_i\} {bi​},求 { a i } \{a_i\} ...

  4. 【c语言学习】1 基础环境安装调试

    1-1下载 vs2019 vs2019下载链接https://visualstudio.microsoft.com/zh-hans/vs/community/ 1-2安装配置环境 记得勾选上c++开发 ...

  5. 第四十六篇:工程化概念以及什么是webpack

    好家伙,这波是概念补充 1.什么是工程化概念? 我的开发: 开个项目,想怎么改怎么改,拉个东西过来就用 实际的前端开发: (1) 模块化(js的模块化,css的模块化,资源的模块化) (2) 组件化( ...

  6. Java开发学习(二十九)----Maven依赖传递、可选依赖、排除依赖解析

    现在的项目一般是拆分成一个个独立的模块,当在其他项目中想要使用独立出来的这些模块,只需要在其pom.xml使用<dependency>标签来进行jar包的引入即可. <depende ...

  7. 【短道速滑九】仿halcon中gauss_filter小半径高斯模糊优化的实现

    通常,我们谈的高斯模糊,都知道其是可以行列分离的算法,现在也有着各种优化算法实现,而且其速度基本是和参数大小无关的.但是,在我们实际的应用中,我们可能会发现,有至少50%以上的场景中,我们并不需要大半 ...

  8. PostgreSQL 涉及复杂视图查询的优化案例

    一.前言 对于含有union , group by 等的视图,我们称之为复杂视图. 这类的视图会影响优化器对于视图的提升,也就是视图无法与父查询进行合并,从而影响访问路径.连接方法.连接顺序等.本文通 ...

  9. Java 监控直播流rtsp协议转rtmp、hls、httpflv协议返回浏览器

    Java 监控直播流rtsp协议转rtmp.hls.httpflv协议返回浏览器 目录 需求背景: 一:了解音视频流协议: 二:方案一 rtsp 转rtmp 1.下载nginx + nginx-rtm ...

  10. MySQL 不同隔离级别,都使用了什么锁?

    大家好,我是树哥. 在上篇文章,我们聊了「MySQL 啥时候会用表锁,啥时候用行锁」这个问题.在文章中,我们还留了一个问题,即:如果查询或更新时的数据特别多,是否从行锁会升级为表锁?此外,还有朋友留言 ...