[笔记] 平衡树合集(Treap,Splay,替罪羊,FHQ Treap)
今天翻了翻其他大佬的博客,发现自己有些。。。颓废。。。
有必要洗心革面,好好学习
序:正常的BST有可能退化,成为链,大大降低效率,所以有很多方法来保持左右size的平衡,本文将简单介绍Treap,Splay,替罪羊,FHQ Treap;
另:代码都是普通平衡树
1.Treap
树堆,在数据结构中也称Treap,是指有一个随机附加域满足堆的性质的二叉搜索树,其结构相当于以随机数据插入的二叉搜索树。其基本操作的期望时间复杂度为O(logn)。相对于其他的平衡二叉搜索树,Treap的特点是实现简单,且能基本实现随机平衡的结构。——百度百科
好的treap=tree+heap(为何不叫hee??)
首先易知堆是棵二叉树,BST也是棵二叉树,
又易知:当堆的中的数据是随机插入(即不是有序数据&&有序插入),堆的的高度是趋于log级别的
于是我们让BST中的节点满足堆性质,让BST中的每一个节点带上一个随机权值dat,作为他在这个满足堆性质的BST中的优先级;
然后为了让BST中的节点满足堆性质,我们要rotate(旋)他
易知以下两种是等价的
仍然满足BST的性质,但是改变了父子关系。
这就是如何在BST中维护堆性质:旋,改变父子关系,直到满足堆性质
PS:此处的旋好像叫单旋,只会改变父子关系,而Splay有一种操作较双旋(见下)
rotate(旋)在一类BST中我认为是最重要的操作
上代码:
ch[x][0/1]左右儿子,vl[x]权值,dat[x]在堆中的优先级,sz[x]子树大小,cnt[x]是vl[x]出现的次数
#include<cstdio> #include<iostream> #include<cstdlib> #define ls ch[x][0] #define rs ch[x][1] #define R register int using namespace std; const int N=100010,Inf=0x3f3f3f3f; inline int g() { R ret=0,fix=1; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-1:fix; do ret=ret*10+(ch^48); while(isdigit(ch=getchar())); return ret*fix; } int n,tot,rt; int sz[N],ch[N][2],vl[N],dat[N],cnt[N]; inline void upd(int x) {sz[x]=sz[ls]+sz[rs]+cnt[x];} inline int cre(int v) {R x=++tot; cnt[x]=1,vl[x]=v,dat[x]=rand(),upd(x); return tot;} inline void rot(int& x,int d) { R y=ch[x][d]; ch[x][d]=ch[y][d^1]; ch[y][d^1]=x; upd(x),upd(y); x=y; } inline void ins(int& x,int v) { if(!x) {x=cre(v); return ;} if(vl[x]==v) {++cnt[x]; upd(x); return ;} R d=vl[x]<v; ins(ch[x][d],v); upd(x); if(dat[ch[x][d]]<dat[x]) rot(x,d); } inline void del(int& x,int v) { if(!x) return ; if(vl[x]==v) { if(cnt[x]>1) --cnt[x]; else { if(!ls) x=rs; else if(!rs) x=ls; else {R d=dat[ls]>dat[rs]; rot(x,d); del(ch[x][d^1],v);}//看谁大就把谁旋上来,把根旋下去 } } else del(ch[x][vl[x]<v],v); upd(x); } inline void build() {srand(100023323); rt=cre(-Inf); ins(rt,Inf);} inline int getpre(int x,int v) { if(!x) return -Inf; if(vl[x]<v) return max(getpre(rs,v),vl[x]);//右边可能没有 else return getpre(ls,v); } inline int getnxt(int x,int v) { if(!x) return Inf; if(vl[x]>v) return min(getnxt(ls,v),vl[x]);//同上 else return getnxt(rs,v); } inline int getrk(int x,int v) { if(!x) return 0; if(vl[x]==v) return sz[ls]+1; else if(vl[x]>v) return getrk(ls,v); else return sz[ls]+cnt[x]+getrk(rs,v); } inline int getvl(int x,int rk) { if(!x||!rk) return 0; if(rk<=sz[ls]) return getvl(ls,rk); else if(rk<=sz[ls]+cnt[x]) return x; return getvl(rs,rk-sz[ls]-cnt[x]); } signed main() { //freopen("in.in","r",stdin);freopen("out.out","w",stdout); n=g(); for(R i=1;i<=n;++i) { R k=g(),x=g(); if(k==1) ins(rt,x); else if(k==2) del(rt,x); else if(k==3) printf("%d\n",getrk(rt,x)); else if(k==4) printf("%d\n",vl[getvl(rt,x)]); else if(k==5) printf("%d\n",getpre(rt,x)); else printf("%d\n",getnxt(rt,x)); } //while(1); }
2.Splay
伸展树(Splay)是一种平衡二叉树,即优化后的二叉查找树。伸展树可以自我调整,这就要依靠伸展操作Splay(x,S),使得提升效率。——洛谷日报
Splay,伸展树。。。维持左右子树平衡用到了另一种旋:双旋
设fa[x]=y,fa[y]=g
双旋,同时改变x,y与y,g之间的关系,它会使g变成x的孙子,y变为x的孩子,g变为y的孩子
分两种情况:
第一种如下
此时要先旋fa,再旋son(纯手绘不喜勿喷qwq);
第二种如下(不在一条链)
那么我们旋两次son
然后Splay的思路是:不管如何操作,将操作的点通过两种旋法,旋至根节点
至于为什么这么旋请找tarjan。。。至于时间复杂度请找tarjan
#include<cstdio> #include<iostream> #define R register int #define ls (ch[x][0]) #define rs (ch[x][1]) const int N=100005,Inf=0x3f3f3f3f; using namespace std; inline int g() { R ret=0,fix=1; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-1:fix; do ret=ret*10+(ch^48); while(isdigit(ch=getchar())); return ret*fix; } int n,tot; int fa[N],ch[N][2],sz[N],vl[N],cnt[N]; inline int cre(int v) {vl[++tot]=v,sz[tot]=cnt[tot]=1; return tot;} inline void upd(int x) {sz[x]=sz[ls]+cnt[x]+sz[rs];} inline void rot(int x) { R y=fa[x],d=ch[y][1]==x; if(fa[y]) ch[fa[y]][ch[fa[y]][1]==y]=x; fa[x]=fa[y]; fa[ch[y][d]=ch[x][d^1]]=y; fa[ch[x][d^1]=y]=x; upd(y); } int rt; inline void print(int x) { if(!x) return ; print(ls); printf("%d\n",vl[x]); print(rs); } inline void Splay(int x,int f) { while(fa[x]!=f) { R y=fa[x]; if(fa[y]!=f) rot((ch[y][1]==x)==(ch[fa[y]][1]==y)?y:x); //在不在一条链上 rot(x); } upd(x); if(!f) rt=x; } inline void ins(int v) { R x=rt; while(1) { if(vl[x]==v) {++cnt[x]; break;} if(!ch[x][vl[x]<v]) { fa[ch[x][vl[x]<v]=cre(v)]=x; x=tot; break; } x=ch[x][vl[x]<v]; } Splay(x,0); } inline void build() {rt=cre(Inf),ins(-Inf);} inline int getrk(int v) { R x=rt,ret=0; while(1) { if(vl[x]==v) {ret+=sz[ls]+1; Splay(x,0); return ret;} if(vl[x]<v) ret+=sz[ls]+cnt[x]; if(!ch[x][vl[x]<v]) {++ret; Splay(x,0); return ret;} x=ch[x][vl[x]<v]; } } inline int getpos(int x,int k) { if(!x) return 0; if(k<=sz[ls]) return getpos(ls,k); if(k<=sz[ls]+cnt[x]) return x; return getpos(rs,k-sz[ls]-cnt[x]); } inline int getvl(int rk) {R x=getpos(rt,rk); Splay(x,0); return vl[x];} inline int getmx(int x,int y) {if(!x||!y) return x|y; return vl[x]>vl[y]?x:y;} inline int ppos(int x,int v) { if(!x) return 0; if(vl[x]<v) return getmx(x,ppos(rs,v)); return ppos(ls,v); } inline int getpre(int v) {R x=ppos(rt,v); Splay(x,0); return vl[x];} inline int getmn(int x,int y) {if(!x||!y) return x|y; return vl[x]<vl[y]?x:y;} inline int npos(int x,int v) { if(!x) return 0; if(v<vl[x]) return getmn(x,npos(ls,v)); return npos(rs,v); } inline int getnxt(int v) {R x=npos(rt,v); Splay(x,0); return vl[x];} inline void del(int v) { Splay(ppos(rt,v),0),Splay(npos(rt,v),rt); R& x=ch[ch[rt][1]][0]; if(!(--cnt[x])) x=0; else Splay(x,0); } signed main() { //freopen("in.in","r",stdin); R n=g(); build(); while(n--) { R k=g(),x=g(); if(k==1) ins(x); else if(k==2) del(x); else if(k==3) printf("%d\n",getrk(x)-1); else if(k==4) printf("%d\n",getvl(x+1)); else if(k==5) printf("%d\n",getpre(x)); else printf("%d\n",getnxt(x)); } //system("pause"); while(1); }
3.替罪羊树
为何叫替罪羊。。。据说拍扁他是他儿子的锅。。。
替罪羊维护左右孩子平衡思路:当max(size(x.ls),size(x.rs))>size(x)*alpha(一个常量,一般0.7-0.8,看个人的喜好。。。),就暴力重构以x为根的子树。
具体地,就是把树拍扁,扔到数组中sort一遍,然后选mid,递归左子树和右子树;
然而替罪羊的删除是懒惰删除。。就是打一个tag。、所以用到儿子时要向下传递。。暴力重构时把删除的节点扔到内存池里去
所以, 特别地,当整个树实际存在的的节点数<整个树的节点数*B(另一个常量,合法范围是0.0-1.0,至于取多少看个人)
变量:sum总节点数=实际存在的节点数+删除的节点数; sz实际存在的节点数;del删除标记;mem存储删除的或没有使用的点;tmp存储拍扁重构的的点
#include<cstdio> #include<iostream> #define R register int using namespace std; const double A=0.72,RB=0.53; const int N=100010; inline int g() { R ret=0,fix=1; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-1:fix; do ret=ret*10+(ch^48); while(isdigit(ch=getchar())); return ret*fix; } struct node{ int ls,rs,vl,sz,sum,del; #define ls(x) t[x].ls #define rs(x) t[x].rs #define vl(x) t[x].vl #define sz(x) t[x].sz #define sum(x) t[x].sum #define del(x) t[x].del }t[N]; int n,rt; int mem[N],cm,tmp[N],ct; inline bool ck(int x) {return (double)sz(x)*A<=(double)max(sz(ls(x)),sz(rs(x)));} inline void dfs(int x) { if(!x) return ; dfs(ls(x)); if(!del(x)) tmp[++ct]=x; else mem[++cm]=x; dfs(rs(x)); } inline void build(int& x,int l,int r) { R md=l+r>>1; x=tmp[md]; if(l==r) { ls(x)=rs(x)=del(x)=0; sz(x)=sum(x)=1; return ; } if(l<md) build(ls(x),l,md-1); else ls(x)=0; build(rs(x),md+1,r); sz(x)=sz(ls(x))+sz(rs(x))+1, sum(x)=sum(ls(x))+sum(rs(x))+1; } inline void rebuild(int& x) { ct=0; dfs(x); if(ct) build(x,1,ct); else x=0; } inline void ins(int& x,int vl) { if(!x) { x=mem[cm--]; vl(x)=vl,ls(x)=rs(x)=del(x)=0; sz(x)=sum(x)=1; return ; } ++sz(x),++sum(x); if(vl(x)>=vl) ins(ls(x),vl); else ins(rs(x),vl); if(ck(x)) rebuild(x); } inline int getrk(int vl) { R x=rt; R ret=1; while(x) { if(vl(x)>=vl) x=ls(x); else {ret+=sz(ls(x))+(del(x)==0); x=rs(x);} } return ret; } inline int getvl(int rk) { R x=rt; while(x) { //cout<<x<<" "<<vl(x)<<endl; if(del(x)==0&&sz(ls(x))+1==rk) return vl(x); else { if(sz(ls(x))+1>rk) x=ls(x); else { rk-=sz(ls(x))+(del(x)==0); x=rs(x); } } } } inline void delrk(int& x,int rk) { if(del(x)==0&&sz(ls(x))+1==rk) {del(x)=1; --sz(x); return ;} --sz(x); if(sz(ls(x))+(del(x)==0)>=rk) delrk(ls(x),rk); else delrk(rs(x),rk-sz(ls(x))-(del(x)==0)); } inline void delvl(int vl) { R x=getrk(vl); //cerr<<x<<endl; delrk(rt,x); if(sum(rt)*RB>=sz(rt)) rebuild(rt); } signed main() { //freopen("in.in","r",stdin); freopen("out.out","w",stdout); n=g(); for(R i=100000;i>=1;--i) mem[++cm]=i; while(n--) { R k=g(),x=g(); if(k==1) ins(rt,x); else if(k==2) delvl(x); else if(k==3) printf("%d\n",getrk(x)); else if(k==4) printf("%d\n",getvl(x)); else if(k==5) printf("%d\n",getvl(getrk(x)-1)); else if(k==6) printf("%d\n",getvl(getrk(x+1))); } //while(1); }
4.FHQ Treap
首先%%%fhq%%%
其次就是这个Treap不用旋,很好理解。。。
两个操作:split和merge
1.split(int o,int v,int& x,int& y) 把以o为根的树分<=v和>v的两部分
也可按照rank去分成两部分
2.merge(int x,int y) 合并以x和y为根的子树,返回新的根的值
短小而精悍。。。
#include<cstdio> #include<iostream> #include<cstdlib> #define R register int #define ls(i) ch[i][0] #define rs(i) ch[i][1] using namespace std; const int N=500010; inline int g() { R ret=0,fix=1; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-1:fix; do ret=ret*10+(ch^48); while(isdigit(ch=getchar())); return ret*fix; } int tot; int ch[N][2],sz[N],vl[N],dat[N]; inline void upd(int x) {sz[x]=sz[ls(x)]+sz[rs(x)]+1;} inline int cre(int v) {R x=++tot; sz[x]=1,vl[x]=v,dat[x]=rand();return x;} inline int merge(int x,int y) { if(!x||!y) return x+y; if(dat[x]<dat[y]) {rs(x)=merge(rs(x),y); upd(x); return x;} else {ls(y)=merge(x,ls(y)); upd(y); return y;} } inline void split(int o,int v,int& x,int& y) { if(!o) {x=y=0; return ;} if(vl[o]<=v) x=o,split(rs(o),v,rs(o),y); else y=o,split(ls(o),v,x,ls(o)); upd(o); } inline int getvl(int x,int rk) { while(1) { if(rk<=sz[ls(x)]) x=ls(x); else if(rk==sz[ls(x)]+1) return x; else rk-=sz[ls(x)]+1,x=rs(x); } } signed main() { srand(100023323u); R n,x,y,z,rt=0; //freopen("1.in","r",stdin);freopen("out.out","w",stdout); n=g(); while(n--) { //cerr<<n<<" "<<rt<<" "<<endl; R k=g(),a=g(); if(k==1) split(rt,a,x,y),rt=merge(merge(x,cre(a)),y); else if(k==2) {split(rt,a,x,z),split(x,a-1,x,y); y=merge(ls(y),rs(y)),rt=merge(merge(x,y),z);} else if(k==3) split(rt,a-1,x,y),printf("%d\n",sz[x]+1),rt=merge(x,y); else if(k==4) printf("%d\n",vl[getvl(rt,a)]); else if(k==5) split(rt,a-1,x,y),printf("%d\n",vl[getvl(x,sz[x])]),rt=merge(x,y); else if(k==6) split(rt,a,x,y),printf("%d\n",vl[getvl(y,1)]),rt=merge(x,y); } while(1); }
2019.05.05&&2019.05.06