无旋treap的简单思想以及模板
因为学了treap,不想弃坑去学splay,终于理解了无旋treap...
好像普通treap没卵用。。。(再次大雾)
简单说一下思想免得以后忘记。普通treap因为带旋转操作似乎没卵用,而无旋treap可以不旋转。
经典地不能再经典的例题
题目描述
您需要写一种数据结构(可参考题目标题),来维护一些数,其中需要提供以下操作:
-
插入x数
-
删除x数(若有多个相同的数,因只删除一个)
-
查询x数的排名(若有多个相同的数,因输出最小的排名)
-
查询排名为x的数
-
求x的前驱(前驱定义为小于x,且最大的数)
- 求x的后继(后继定义为大于x,且最小的数)
输入输出格式
输入格式:
第一行为n,表示操作的个数,下面n行每行有两个数opt和x,opt表示操作的序号(1<=opt<=6)
输出格式:
对于操作3,4,5,6每行输出一个数,表示对应答案
输入输出样例
10 1 106465 4 1 1 317721 1 460929 1 644985 1 84185 1 89851 6 81968 1 492737 5 493598
106465 84185 492737
说明
时空限制:1000ms,128M
1.n的数据范围:n<=100000
2.每个数的数据范围:[-1e7,1e7]
无旋treap最基本的操作就是merge和split。修改都得靠这个。首先说一下merge和split的操作方法。
merge:给你两个平衡树a和b,让你合并它们,但b中所有权值都得大于a中所有权值以维护treap性质。如果是tree只需要建个根节点然后连两条边,但无法维护heap性质。其实操作也很简单,如果a的rand值小于b的rand值,则merge(a.rs,b),否则merge(a,b.ls)。
split:将平衡树a分成两个,一个是1~K名,一个是K+1~size名,返回两个树的根,我存在了一个pair中。记左子树size为s,如果k==s则两根就是左儿子和根,如果k==s+1则两根就是根和右儿子。若都不是,继续向下递归,那个就参照以下查询排名的操作。注意split要将儿子设为0
那么基本操作怎么实现???
查询排名、查询排名为x的数、查询前驱后继代码基本不变,只是要注意这里需要合并分离所以不可以记录出现次数要新开一个数,否则后果自行脑补。
插入操作:先查询这个数x的排名K,然后把treap分成两个treap,即执行split(root,K),得到两根a,b。再合并两次,即执行root=merge(a,点的新编号),root=merge(root,b)(注意要重设root)
删除操作:也差不多。先查询排名K,然后把treap分成三个treap,分别为a=1~K-1,b=K,c=K+1~n,要删去K,只需再合并a,c子树即可
不旋转应该就可以维护一些别的好东西了
然而并不会
然后附上普通平衡树的AC代码:
1 // It is made by XZZ 2 #include<cstdio> 3 #include<algorithm> 4 using namespace std; 5 #define rep(a,b,c) for(rg int a=b;a<=c;a++) 6 #define drep(a,b,c) for(rg int a=b;a>=c;a--) 7 #define erep(a,b) for(rg int a=fir[b];a;a=nxt[a]) 8 #define il inline 9 #define rg register 10 #define vd void 11 #define mp make_pair 12 typedef long long ll; 13 il int gi(){ 14 rg int x=0,f=1;rg char ch=getchar(); 15 while(ch<'0'||ch>'9')f=ch=='-'?-1:f,ch=getchar(); 16 while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar(); 17 return x*f; 18 } 19 #define Now tree[now] 20 struct node{int ls,rs,value,rand,size;}tree[100010]; 21 int root,siz; 22 int seed=151806+150605+151127;//mou shen ben xue hao+mou da lao xue hao+mou ju ruo xue hao 23 il int Rand(){return seed=seed*48271%2147483647;} 24 il vd reset(int now){Now.size=tree[Now.ls].size+tree[Now.rs].size+1;} 25 il int merge(int a,int b){ 26 if(!a||!b)return a|b; 27 if(tree[a].rand<tree[b].rand){tree[a].rs=merge(tree[a].rs,b),reset(a);return a;} 28 else {tree[b].ls=merge(a,tree[b].ls),reset(b);return b;} 29 } 30 il pair<int,int>split(int now,int num){ 31 if(!now)return mp(0,0); 32 int ls=Now.ls,rs=Now.rs; 33 if(num==tree[ls].size){Now.ls=0,reset(now);return mp(ls,now);} 34 if(num==tree[ls].size+1){Now.rs=0,reset(now);return mp(now,rs);} 35 if(num<tree[ls].size){ 36 pair<int,int>T=split(ls,num); 37 Now.ls=T.second,reset(now); 38 return mp(T.first,now); 39 }else{ 40 pair<int,int>T=split(rs,num-tree[ls].size-1); 41 Now.rs=T.first,reset(now); 42 return mp(now,T.second); 43 } 44 } 45 il int getrank(int now,int num){ 46 int ret=0,t=1e9; 47 while(now){ 48 if(num==Now.value)t=min(t,ret+tree[Now.ls].size+1); 49 if(num<=Now.value)now=Now.ls; 50 else ret+=tree[Now.ls].size+1,now=Now.rs; 51 } 52 return t==1e9?ret:t; 53 } 54 il int getnum(int now,int num){ 55 while(1){ 56 if(tree[Now.ls].size==num-1)return Now.value; 57 if(tree[Now.ls].size>num-1)now=Now.ls; 58 else num-=tree[Now.ls].size+1,now=Now.rs; 59 } 60 } 61 il int lower(int now,int num){ 62 int ret; 63 while(now)if(tree[now].value<num)ret=tree[now].value,now=Now.rs; 64 else now=Now.ls; 65 return ret; 66 } 67 il int upper(int now,int num){ 68 int ret; 69 while(now)if(tree[now].value>num)ret=tree[now].value,now=Now.ls; 70 else now=Now.rs; 71 return ret; 72 } 73 il vd ins(int num){ 74 int Rank=getrank(root,num),now; 75 pair<int,int>tmp=split(root,Rank); 76 now=++siz; 77 Now.value=num,Now.rand=Rand(),Now.size=1; 78 root=merge(tmp.first,siz); 79 root=merge(root,tmp.second); 80 } 81 il vd del(int num){ 82 int Rank=getrank(root,num); 83 pair<int,int>t1=split(root,Rank),t2=split(t1.first,Rank-1); 84 root=merge(t2.first,t1.second); 85 } 86 int main(){ 87 int n=gi(),opt,x; 88 while(n--){ 89 opt=gi(),x=gi(); 90 switch(opt){ 91 case 1:ins(x);break; 92 case 2:del(x);break; 93 case 3:printf("%d\n",getrank(root,x));break; 94 case 4:printf("%d\n",getnum(root,x));break; 95 case 5:printf("%d\n",lower(root,x));break; 96 case 6:printf("%d\n",upper(root,x));break; 97 } 98 } 99 return 0; 100 }
还有一篇博客代码写的不错,这是链接
本博客中博文均为原创,未经博主允许请勿随意转载,谢谢。