平衡树之替罪羊树

平衡树之替罪羊树

–yangkai


身为平衡树却不做任何形式的旋转,替罪羊树可以称得上是最暴力的平衡树了。

替罪羊树(SGT)保留有二叉搜索树的基本性质,即对于任意一个节点t,左儿子的所有节点比它小,右儿子的所有节点比它大。但是既然不基于翻转,它怎样维护平衡树的优秀复杂度呢?

SDT基于一个叫做“重构”的操作,听起来很是优美暴力,那么我们要如何重构?

首先,如果要时刻维持平衡树的平衡,即 a b s ( s i z [ l ] s i z [ r ] ) <= 1 ,显然我们平均每两次插入、删除操作就需要一次重构,这样的时间效率接近 O ( n 2 ) ,所以我们可以取一个阈值[0.5,1.0],当且仅当左右子树超过阈值的限制重构子树。在这里我们可以发现,当阈值限制越紧[0.5,0.75],对修改少询问多有利;当阈值限定越松[0.75,1.0],对修改多查询少有利。

根据阈值,我们就可以对子树进行重构了。

但是问题又来了,在条根到叶子的链上可能有许许多多的“不平衡点”,我们要怎样确定最优的一个,实际上,我们只需要选择最靠近根节点的一个就好了,因为如果修改的层数太深,在总体上不会对树的结构产生多么大的影响,所以选最靠近根的一个相对来说比较优秀。

那么维护平衡树性质的部分我们解决了,但是插入以及删除呢?

首先,对于插入,直接参见动态开点大法,在向下插入的时候动态开点就好了(此处应该注意回收内存)

其次,对于删除,我们既不能把它转到叶子结点(treap),也不能直接裂开丢掉再合并(非旋转treap),所以我们给它打上一个real标记,real标记存在说明这个节点信息为’真‘,否则为‘假’(已经删除),然后在重构子树的时候直接丢掉然后回收内存就好了

其他操作rank,get_kth,pre,nxt之类的操作直接套用平衡树常见套路就好了


例题:BZOJ3224 Tyvj 1728 普通平衡树

Description

您需要写一种数据结构(可参考题目标题),来维护一些数,其中需要提供以下操作:
\1. 插入x数
\2. 删除x数(若有多个相同的数,因只删除一个)
\3. 查询x数的排名(若有多个相同的数,因输出最小的排名)
\4. 查询排名为x的数
\5. 求x的前驱(前驱定义为小于x,且最大的数)
\6. 求x的后继(后继定义为大于x,且最小的数)

Input

第一行为n,表示操作的个数,下面n行每行有两个数opt和x,opt表示操作的序号(1<=opt<=6)

Output

对于操作3,4,5,6每行输出一个数,表示对应答案

Sample Input

10
1 106465
4 1
1 317721
1 460929
1 644985
1 84185
1 89851
6 81968
1 492737
5 493598

Sample Output

106465
84185
492737

HINT

1.n的数据范围:n<=100000

2.每个数的数据范围:[-2e9,2e9]


上板子

//yangkai
#include<bits/stdc++.h>
using namespace std;
inline int read(){
    int data=0,w=1;char ch=getchar();
    while(!isdigit(ch)&&ch!='-')ch=getchar();
    if(ch=='-')w=-1,ch=getchar();
    while(isdigit(ch))data=(data<<1)+(data<<3)+(ch-'0'),ch=getchar();
    return data*w;
}
struct SGT{
    static const int N=100010;
    static const double alpha=0.75;
    //内存池 
    int st[N],top,tot;
    int get_place(){return top?st[top--]:++tot;}
    void del_place(int t){st[++top]=t;}
    void clean_st(){top=tot=0;}
    //SGT
    int root;
    int ls[N],rs[N];
    int val[N],siz[N],all[N];
    bool real[N];
    void pushup(int t){
        siz[t]=siz[ls[t]]+siz[rs[t]]+real[t];
        all[t]=all[ls[t]]+all[rs[t]]+1;
    }
    bool check(int t){//返回1说明需要重构 
        if(all[ls[t]]>=alpha*all[t])return 1;
        if(all[rs[t]]>=alpha*all[t])return 1;
        return 0;
    }
    int newnode(int vl=0){
        int t=get_place();
        ls[t]=rs[t]=0;
        val[t]=vl;
        siz[t]=all[t]=1;
        real[t]=1;
        return t;
    }
    //收集需要重构的下标 
    void collect(int t,vector<int> &v){
        if(!t)return;
        collect(ls[t],v);
        if(real[t])v.push_back(t);
        else del_place(t);
        collect(rs[t],v); 
    }
    //重构划分 
    int divide(int l,int r,vector<int> v){
        if(l>=r)return 0;
        int mid=(l+r)>>1;
        int t=v[mid];
        ls[t]=divide(l,mid,v);
        rs[t]=divide(mid+1,r,v);
        pushup(t);
        return t;
    }
    //重新建树 
    void rebuild(int &t){
        static vector<int> v;v.clear();
        collect(t,v);
        //创造临时节点 
        int newt=divide(0,v.size(),v);
        t=newt;
    }
    //查询vl的排名 
    int rank(int vl){
        int t=root,ans=1;
        while(t){
            if(vl<=val[t])t=ls[t];
            else{
                ans+=siz[ls[t]]+real[t];
                t=rs[t];
            }
        }
        return ans;
    }
    //查询排名为k的数
    int get_kth(int k){
        int t=root;
        while(t){
            if(siz[ls[t]]+1==k&&real[t])return val[t];
            if(siz[ls[t]]>=k)t=ls[t];
            else k-=siz[ls[t]]+real[t],t=rs[t];
        }
    } 
    //插入vl值 
    int insert(int &t,int vl){//返回深度最浅的不合法点 
        if(!t){
            t=newnode(vl);
            return 0;
        }
        siz[t]++;all[t]++;
        int res;
        if(vl<=val[t])res=insert(ls[t],vl);
        else res=insert(rs[t],vl);
        pushup(t);
        //检查是否需要重构 
        if(check(t))res=t;
        return res;
    }
    void insert(int vl){
        int t=insert(root,vl);
        if(!t)rebuild(t);
    }
    //删除排名为k的数 
    void erase(int t,int k){
        siz[t]--;
        if(real[t]&&k==siz[ls[t]]+real[t]){
            real[t]=0;
            return;
        }
        if(k<=siz[ls[t]])erase(ls[t],k);
        else erase(rs[t],k-siz[ls[t]]-real[t]);
    }
    //删除数vl 
    void erase(int vl){
        erase(root,rank(vl));
        if(siz[root]<alpha*all[root])rebuild(root);
    }
    //查询vl的前驱和后继 
    int pre(int vl){return get_kth(rank(vl)-1);} 
    int nxt(int vl){return get_kth(rank(vl+1));}
    //输出调试
    void out_put(int t){
        cout<<t<<"||"<<endl;
        if(ls[t])out_put(ls[t]);
        if(real[t])cout<<val[t]<<" L:"<<val[ls[t]]<<" R:"<<val[rs[t]]<<endl;
        if(rs[t])out_put(rs[t]);
    } 
}sgt;
int main(){
    int n=read();
    while(n--){
        int op=read(),t=read();
        switch(op){
            case 1:sgt.insert(t);break;
            case 2:sgt.erase(t);break;
            case 3:printf("%d\n",sgt.rank(t));break;
            case 4:printf("%d\n",sgt.get_kth(t));break;
            case 5:printf("%d\n",sgt.pre(t));break;
            case 6:printf("%d\n",sgt.nxt(t));break;
        }
    }
    return 0;
} 

猜你喜欢

转载自blog.csdn.net/dream_maker_yk/article/details/80556381