平衡树之非旋Treap
平衡树(二叉树)
线段树不支持插入or删除一个数于是平衡树产生了
常见平衡树:treap(比sbt慢,好写吧),SBT(快,比较好写,有些功能不支持),splay(特别慢,复杂度当做根号n来用,功能强大,不好写),rbt(红黑树,特别快),//替罪羊树,朝鲜树
晚上要讲的不旋转平衡树:
平衡树:
节点的左儿子中的每一个一定比他小,右儿子中的每一个一定比他大
那么它的中序遍历是有序的
用下标建树,那么区间询问的话就是求一棵子数和子树根和领一棵子数的一部分
treap:
tree+heap,平衡树和heap的性质是矛盾的,所以每个节点存一个key和value
key值满足heap性质,value满足平衡树的性质,这样的树叫做treap?
插入:
插入的新节点的key值随机,调用rand函数(这样保证树的深度一定是logn的)改变树的形态使它重新满足hea与平衡树性质
操作1.merge:
merge(P1,P2):把以p1为根的treap和以p2为根的treap合并成一个treap(p1中的所有制小于
操作2.split:
把以p为根的treap中拿出k小的数,组成一个新treap
保证原先树中的所有数>新树中所有数
可持久化treap :
插:
建一个只有一个点的树(要插得数)例如(2.33)把(1,2)splay出来,再把新树(2.33)和(1,2)merge起来,再把(1,2,2.33)和(4,5)merge 一下
删除一个:
如删除(2.33),先把split(treap,3),此时把splay把(1,2)与(2.33,4,5)分离在split(treep2,1),此时(2.33)与(4,5)分离
在merge(treap1,treap3)合并即把(1,2),(4,5)合并,那么2.33就没了
实际操作
merge时,找key值最大的作为新treap的根,不是p1就是p2
1要是p1.p>=p2.p此时p1作为新根,那么p1的左儿子不会变换,右子树就是p1的右子树和p2 merge 一下,即 merge(p1.r,p2);
2要是p2.p>p1.p此时p2作为新根,那么p2的右儿子不会变换,左儿子就是p2的
左子树 和 p1 mege 一下 即 merge(p2.l,p1);
split(p,k)几点记录value,key,l,r,size
p.L<-p->p.r;
1.要是k<=p.l.size 说明k小的点全在左子数,递归split(p.l,k);构成新树的时候直接把split后剩下的左子树接到P根上就好了
2.k=p.l.size+1;,返回两棵树(p.l-p,p.r)
3.k>p.l.siz+1,左边已经全不要,那么就split(p.r,k-p.l.size-1);
返回两棵树(p.l-p-p.r,剩余p.r)
merge:
int merge(int x,int y) {
if(!x) return y;//left son is empty
else if(!y)return x;//right son is empty
if(key[x]<key[y]) {
tree[x][1]=merge(tree[x][1],y);//union y_tree to x's right son
update(x);
return x;
}
else {
tree[y][0]=merge(x,tree[y][0]);//union x_tree to y'left son
update(y);
return y;
}
}
split:
void split(int now,int k,int &x,int &y) {//根据插入数的val值,把树分为小于等于k的,和大于k的
if (!now) x=y=0;//指针为空
else {
if (val[now]<=k)
x=now,split(tree[now][1],k,tree[now][1],y);//分割标准大于当前节点val访问右子树
else
y=now,split(tree[now][0],k,x,tree[now][0]);//否则访问左子树
update(now);
}
}
query_min:
查询那些数比x数小,当找到一个根节点比x小时,那么该节点的所有子树都比他小,那么就把子树size+1加到答案里-->删除一个数的时候时用来确定split的k(比要删除的数小的)值
query_kth
查询第k大的数
int kth(int now,int k) {
while(1) {
if (k<=siz[tree[now][0]])//左子树大小比k大说明第k大在左子树中
now=tree[now][0];
else
if (k==siz[tree[now][0]]+1)
return now;
else
k-=siz[tree[now][0]]+1,now=tree[now][1];//在右子树中查询第k-左子树子树个数 大的数
}
}
查询数a的rank
按照权值进行split,split出的x数(所有节点权值小于a)的大小+1,就是a的rank
split(root,a-1,x,y);
printf("%d\n",siz[x]+1);
root=merge(x,y);
插入一个新节点
int make_new(int v) {
siz[++sz]=1;
val[sz]=v;
key[sz]=rand();
return sz;
}
delite a number
删除一个数,把这个数split出来,然后将另外两颗split出来的树merge起来
split(root,a,x,z);
split(x,a-1,x,y);
y=merge(tree[y][0],tree[y][1]);
root=merge(merge(x,y),z);
查询前驱
查询a的前驱
查询按照a作为标准split出的x子树中的最大值
split(root,a-1,x,y);
printf("%d\n",val[kth(x,siz[x])]);
root=merge(x,y);
后继
查询a的后继
查询按照a作为标准split出的y子树中的最小值
split(root,a,x,y);
printf("%d\n",val[kth(y,1)]);
root=merge(x,y);
区间操作
总结,利用split与merge完成各种操作
平衡树之非旋Treap的更多相关文章
- [模板] 平衡树: Splay, 非旋Treap, 替罪羊树
简介 二叉搜索树, 可以维护一个集合/序列, 同时维护节点的 \(size\), 因此可以支持 insert(v), delete(v), kth(p,k), rank(v)等操作. 另外, prev ...
- 2018.08.05 bzoj3223: Tyvj 1729 文艺平衡树(非旋treap)
传送门 经典的平衡树问题,之前已经用splay写过一次了,今天我突发奇想,写了一发非旋treap的版本,发现挺好写的(虽然跑不过splay). 代码: #include<bits/stdc++. ...
- bzoj 3224: Tyvj 1728 普通平衡树【非旋treap】
就是非旋treap的板子 #include<iostream> #include<cstdio> #include<cstdlib> using namespace ...
- 平衡树简单教程及模板(splay, 替罪羊树, 非旋treap)
原文链接https://www.cnblogs.com/zhouzhendong/p/Balanced-Binary-Tree.html 注意是简单教程,不是入门教程. splay 1. 旋转: 假设 ...
- 非旋 treap 结构体数组版(无指针)详解,有图有真相
非旋 $treap$ (FHQ treap)的简单入门 前置技能 建议在掌握普通 treap 以及 左偏堆(也就是可并堆)食用本blog 原理 以随机数维护平衡,使树高期望为logn级别, FHQ ...
- 2827: 千山鸟飞绝 非旋treap
国际惯例的题面:看起来很不可做的样子,我们先来整理一下题意吧.就是,维护每个点曾经拥有过的最大的两个属性值,支持把点的位置移动.我们用map对每个位置进行离散化,对每个位置建立一个平衡树.为了方便分离 ...
- 2018.08.06 bzoj1500: [NOI2005]维修数列(非旋treap)
传送门 平衡树好题. 我仍然是用的fhqtreap,感觉速度还行. 维护也比线段树splay什么的写起来简单. %%%非旋treap大法好. 代码: #include<bits/stdc++.h ...
- 2018.07.24 loj#107. 维护全序集(非旋treap)
传送门 就是普通平衡树,可以拿来练非旋treap" role="presentation" style="position: relative;"&g ...
- 2018.07.06 BZOJ 1588: HNOI2002营业额统计(非旋treap)
1588: [HNOI2002]营业额统计 Time Limit: 5 Sec Memory Limit: 162 MB Description 营业额统计 Tiger最近被公司升任为营业部经理,他上 ...
随机推荐
- 将FragmentManger事务添加到返回栈中
FragmentManger事务添加或替换的 Fragment 后,这时点击 Back 键,程序并不会返回添加之前的状态. 我们可以使用 Transaction 对象的 addToBackStack( ...
- leetcode 【 Set Matrix Zeroes 】python 实现
题目: Given a m x n matrix, if an element is 0, set its entire row and column to 0. Do it in place. cl ...
- 微服务学习笔记——Spring Boot特性
1. 创建独立的Spring应用程序 2. 嵌入的Tomcat,无需部署WAR文件 3. 简化Maven配置 4. 自动配置Spring 5. 提供生产就绪型功能,如指标,健康检查和外部配置 6. 开 ...
- 图说不为人知的IT传奇故事-4-王安用一生来跟IBM抗衡
此系列文章为“图说不为人知的IT传奇故事”,各位大忙人可以在一分钟甚至几秒内了解把握整个内容,真可谓“大忙人的福利”呀!!希望各位IT界的朋友在钻研技术的同时,也能在文学.历史上有所把握.了解这些故事 ...
- IOS开发学习笔记035-UIScrollView-自动滚动
让图片自动滚动的话,需要使使用定时器,循环计算当前页的页码.并且在拖动图片时停止计时器,停止拖动时启动计时器. 定时器 方法1: performSelector [self performSelect ...
- PAT1026
要获得一个C语言程序的运行时间,常用的方法是调用头文件time.h,其中提供了clock()函数,可以捕捉从程序开始运行到clock()被调用时所耗费的时间.这个时间单位是clock tick,即“时 ...
- quagga源码学习--BGP协议的初始化
quagga支持BGP-4,BGP-4+协议,支持多协议(mpls,isis,ospf等等)以及单播,组播路由的导入和分发. 具体的协议,这里就不附录了,网络上有很多资料,或者RFC. 协议源码的学习 ...
- HDU 2036 求任意多边形面积向量叉乘
三角形的面积可以使用向量的叉积来求: 对于 三角形的面积 等于: [(x2 - x1)*(y3 - y1)- ( y2 - y1 ) * ( x3 - x1 ) ] / 2.0 但是面积是有方向的, ...
- Codeforces 785D Anton and School - 2 (组合数相关公式+逆元)
D. Anton and School - 2 time limit per test 2 seconds memory limit per test 256 megabytes input stan ...
- Java EE 学习(5):IDEA + maven + spring 搭建 web(1)
参考:http://www.cnblogs.com/lonelyxmas/p/5397422.html http://www.ctolib.com/docs-IntelliJ-IDEA-c--1590 ...