平衡树之非旋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最近被公司升任为营业部经理,他上 ...
随机推荐
- python上数据存储 .h5格式或者h5py
最近在做城市计算的项目,数据文件是以.h5的格式存储的,总结下其用法和特点 来自百度百科的简介: HDF(Hierarchical Data Format),可以存储不同类型的图像和数码数据的文件格式 ...
- ElasticSearch学习笔记(一)-- 查询索引分词
# 查看所有索引 GET _cat/indices # 创建一个索引 PUT /test_index # 插入一条数据(指定id)PUT /test_index/doc/ { "userna ...
- Spring MVC 使用 HttpServletResponseWrapper 修改返回结果
HttpServletResponseWrapper 是什么? ServletResponse 的包装类,相关设计模式 装饰者模式. 运行环境 jdk 1.7 spring boot 整合的web环境 ...
- 【Combinations】cpp
题目: Given two integers n and k, return all possible combinations of k numbers out of 1 ... n. For ex ...
- python批量爬取文档
最近项目需要将批量链接中的pdf文档爬下来处理,根据以下步骤完成了任务: 将批量下载链接copy到text中,每行1个链接: 再读txt文档构造url_list列表,利用readlines返回以行为单 ...
- hdu6212[区间dp] 2017青岛ACM-ICPC网络赛
原题: BZOJ1032 (原题数据有问题) /*hdu6212[区间dp] 2017青岛ACM-ICPC网络赛*/ #include <bits/stdc++.h> using name ...
- [LOJ#2326]「清华集训 2017」简单数据结构
[LOJ#2326]「清华集训 2017」简单数据结构 试题描述 参加完IOI2018之后就是姚班面试.而你,由于讨厌物理.并且想成为乔布斯一样的创业家,被成功踢回贵系. 转眼,时间的指针被指向201 ...
- 洛谷P3832 [NOI2017]蚯蚓排队 【链表 + 字符串hash】
题目链接 洛谷P3832 题解 字符串哈希然后丢到hash表里边查询即可 因为\(k \le 50\),1.2操作就暴力维护一下 经复杂度分析会发现直接这样暴力维护是对的 一开始自然溢出WA了,还以为 ...
- c#的字典序
//Dictionary System.Collections.DictionaryEntry dic=new System.Collections.DictionaryEntry("key ...
- ZJOI2017D1
假装我还活着. 去温州前沉迷各种奇怪的动画片..嗯补了不少高达.. 到温州以后继续看片..嗯ZG还是挺不错的..然后接着就FA♂现我什么都不会写..有点尴尬.. 因为宾馆离温州中学比较远就完全没去听课 ...