树套树模板 之 二逼平衡树

13 篇文章 0 订阅
9 篇文章 0 订阅

题目描述:

描述?描述个屁!
就是让你写个数据结构,资瓷一堆恶心的操作.

题目分析:

如果没有区间限制,这就是道裸的不能再裸的Splay!
然而加入区间限制要求我们可持久化
什么BIT套Seg啊 主席树套BIT啊
Seg套Treap啊 还有分块大佬!
然后我就写了Seg套Splay
(刚开始想搞普通线段树套个权值来着,无奈不会搞前驱后继,如果会了在搞吧qwq
其实树套树并不难理解
这道题目的一维就是一颗Seg 用来维护 区间查询操作
而不同的是,Seg上不是维护着一个个信息,而是直接挂了一颗Splay
这颗Splay存有所有属于本Seg节点区间所有值
然后我们在Splay里来进行维护即可!
第一次打这么长的代码

题目链接:

Luogu 3380
BZOJ 3196

Ac 代码:

#include <cstdio>
#include <iostream>
#include <algorithm>
#define il inline 
const int maxm=4000010;
const int inf=2147483647;
int ans;
int rt[maxm],siz[maxm],cnt[maxm],ch[maxm][2],fa[maxm],data[maxm];
int sz,MX,n,q;
int a[maxm];
il void splayclear(int x)
{
    siz[x]=cnt[x]=ch[x][0]=ch[x][1]=fa[x]=data[x]=0;
}
il void update(int x)
{
    siz[x]=cnt[x];
    if(ch[x][0]) siz[x]+=siz[ch[x][0]];
    if(ch[x][1]) siz[x]+=siz[ch[x][1]];
}
il bool get(int x)
{
    return ch[fa[x]][1]==x;
}
il void rotate(int x)
{
    int fax=fa[x],ffa=fa[fax],which=get(x);
    bool fx=get(fax);
    ch[fax][which]=ch[x][which^1]; 
    fa[ch[fax][which]]=fax;  
    ch[x][which^1]=fax;
    fa[fax]=x;  
    fa[x]=ffa;

    if(ffa)  
     ch[ffa][fx]=x;  
    update(fax),update(x);  
}
inline void splay(int i,int x,int top)
{
     for(int fax;(fax=fa[x])!=0;rotate(x))
       if(fa[fax]!=0)
         rotate(get(fax)==get(x)?fax:x);
    rt[i]=x;
}  
il void splayinsert(int i,int val)
{
    if(!rt[i])
    {
        sz++;
        fa[sz]=ch[sz][0]=ch[sz][1]=0;
        data[sz]=val,cnt[sz]=1,siz[sz]=1;
        rt[i]=sz;
        return;
    }
    int now=rt[i],fax=0; 
    while(1)
    {
        if(data[now]==val)
        {
            cnt[now]++,update(now),update(fax),splay(i,now,0);
            return;
        }
        fax=now;
        now=ch[fax][data[now]<val];
        if(!now)
        {
            sz++;
            ch[sz][0]=ch[sz][1]=0;
            data[sz]=val,cnt[sz]=1,siz[sz]=1;
            fa[sz]=fax,ch[fax][data[fax]<val]=sz;
            update(fax); 
            splay(i,sz,0);
            return;
        }
    }
}
il int splayrank(int i,int val)
{
    int x=rt[i],ret=0;
    while (x)
    {
        if (data[x]==val) return ret+((ch[x][0])?siz[ch[x][0]]:0);
        if (data[x]<val)
        {
            ret+=((ch[x][0])?siz[ch[x][0]]:0)+cnt[x];
            x=ch[x][1];
        } else x=ch[x][0];
    }
    return ret;
}
il int splayfind(int i,int val)
{
    int x=rt[i];
    while(x)
    {
        if(data[x]==val)
        {
            splay(i,x,0);
            return x;
        }
        x=ch[x][val>data[x]];
    }
} 
il int splaypre(int i)
{
    int x=ch[rt[i]][0];
    while(ch[x][1]) x=ch[x][1];
    return x;
}
il int splaynext(int i)
{
    int x=ch[rt[i]][1];
    while(ch[x][0]) x=ch[x][0];
    return x;
}
il void splaydel(int i,int val)
{
    int x=splayfind(i,val);
    if(cnt[x]>1)
    {
        cnt[x]--,update(x);
        return;
    }
    if(!ch[x][0]&&!ch[x][1])
    {
        splayclear(rt[i]);
        rt[i]=0;
        return;
    }
    if (!ch[x][0])
    {
        int y=ch[x][1];
        rt[i]=y,fa[y]=0;
        return;
    }
    if (!ch[x][1])
    {
        int y=ch[x][0]; 
        rt[i]=y,fa[y]=0;
        return;
    }
    int p=splaypre(i);
    int oldrt=rt[i];
    splay(i,p,0);
    ch[rt[i]][1]=ch[oldrt][1];
    fa[ch[oldrt][1]]=rt[i];
    splayclear(oldrt);
    update(rt[i]);
}
il int splaygetpre(int i,int val)
{
    int ret=-inf;
    int x=rt[i];
    while(x)
    {
        if (data[x]<val)
        {
            if(ret<data[x]) ret=data[x];
            x=ch[x][1];
        } else x=ch[x][0];
    } 
    return ret;
}
il int splaygetnext(int i, int v)
{
    int ret=inf;
    int x=rt[i];
    while(x)
    {
        if (data[x]>v)
        {
            if(ret>data[x]) ret=data[x];
            x=ch[x][0];
        } else x=ch[x][1];
    } 
    return ret;
}
il void seginsert(int o,int l,int r,int ind,int num)
{
    splayinsert(o,num);
    if(l>=r) return;
    int mid=(l+r)>>1;
    if(ind<=mid) seginsert((o<<1),l,mid,ind,num);
    else seginsert((o<<1)|1,mid+1,r,ind,num);
}
il void segrank(int o,int l,int r,int ql,int qr,int num)
{
    if(ql<=l&&r<=qr)
    {
        ans+=splayrank(o,num);
        return;
    }
    int mid=(l+r)>>1;
    if(ql<=mid) segrank((o<<1),l,mid,ql,qr,num);
    if(qr>mid) segrank((o<<1)|1,mid+1,r,ql,qr,num);
}
il void segchange(int o,int l,int r,int ind,int ynum,int cnum)
{
    splaydel(o,ynum),splayinsert(o,cnum);
    if(l>=r) return;
    int mid=(l+r)>>1;
    if(ind<=mid) segchange((o<<1),l,mid,ind,ynum,cnum);
    else segchange((o<<1)|1,mid+1,r,ind,ynum,cnum);
}
il void segpre(int o,int l,int r,int ql,int qr,int num)
{
    if(ql<=l&&r<=qr)
    {
        ans=std::max(ans,splaygetpre(o,num));
        return;
    }
    int mid=(l+r)>>1;
    if(ql<=mid) segpre((o<<1),l,mid,ql,qr,num);
    if(qr>mid) segpre((o<<1)|1,mid+1,r,ql,qr,num);
}
il void segnext(int o,int l,int r,int ql,int qr,int num)
{
    if(ql<=l&&r<=qr)
    {
        ans=std::min(ans,splaygetnext(o,num));
        return;
    }
    int mid=(l+r)>>1;
    if(ql<=mid) segnext((o<<1),l,mid,ql,qr,num);
    if(qr>mid) segnext((o<<1)|1,mid+1,r,ql,qr,num);
}
il int getkth(int ql,int qr,int k)
{
    int l=0,r=MX+1;
    while(l<r)
    {
        int mid=(l+r)>>1;
        ans=0;
        segrank(1,1,n,ql,qr,mid);
        if (ans<k) l=mid+1;
        else r=mid;
    }
    return l-1;
}
int main()
{
    //freopen("wc.in","r",stdin);
    //freopen("he.out","w",stdout);
    scanf("%d%d",&n,&q);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        MX=std::max(MX,a[i]);
        seginsert(1,1,n,i,a[i]);
    }
    for(int i=1,l,r,k,opt;i<=q;i++)
    {
        scanf("%d",&opt);
        if(opt==3)
        {
            scanf("%d%d",&l,&k);
            segchange(1,1,n,l,a[l],k);
            a[l]=k;
        }
        else
        {
            scanf("%d%d%d",&l,&r,&k);
            if(opt==1)
            {
                ans=0;
                segrank(1,1,n,l,r,k);
                printf("%d\n",ans+1);
            }
            if(opt==2)
             printf("%d\n",getkth(l,r,k));
            if(opt==4)
            {
                ans=-inf;
                segpre(1,1,n,l,r,k);
                printf("%d\n",ans);
            }
            if(opt==5)
            {
                ans=inf;
                segnext(1,1,n,l,r,k);
                printf("%d\n",ans);
            }
        }
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
平衡树也称为自平衡叉搜索,实现起来比较复杂,需要用到旋转操作来保持的平衡,如左旋、右旋、双旋等。 下面是一份平衡树模板代码,主要包括插入、删除、查找操作,以及各种旋转操作实现。 ```c #include <stdio.h> #include <stdlib.h> // 平衡树节点结构体 typedef struct TreeNode { int val, height; // val:节点权值,height:节点高度 struct TreeNode *left, *right; // 左右子节点指针 } TreeNode; // 获取节点高度 int getHeight(TreeNode *node) { return node ? node->height : 0; } // 获取节点平衡因子 int getBalanceFactor(TreeNode *node) { return getHeight(node->left) - getHeight(node->right); } // 更新节点高度 void updateHeight(TreeNode *node) { node->height = 1 + fmax(getHeight(node->left), getHeight(node->right)); } // 左旋操作 TreeNode *leftRotate(TreeNode *node) { TreeNode *right = node->right; node->right = right->left; right->left = node; updateHeight(node); updateHeight(right); return right; } // 右旋操作 TreeNode *rightRotate(TreeNode *node) { TreeNode *left = node->left; node->left = left->right; left->right = node; updateHeight(node); updateHeight(left); return left; } // 平衡操作(使以node为根的平衡,返回新根节点) TreeNode *balance(TreeNode *node) { // 更新节点高度 updateHeight(node); // 计算平衡因子 int balanceFactor = getBalanceFactor(node); // 如果平衡因子为2,说明左子比右子高 if (balanceFactor > 1) { if (getBalanceFactor(node->left) < 0) { // 如果左子的右子比左子高 node->left = leftRotate(node->left); // 先左旋左子节点 } return rightRotate(node); } // 如果平衡因子为-2,说明右子比左子高 else if (balanceFactor < -1) { if (getBalanceFactor(node->right) > 0) { // 如果右子的左子比右子高 node->right = rightRotate(node->right); // 先右旋右子节点 } return leftRotate(node); } // 否则说明当前子已经平衡,直接返回 return node; } // 插入操作(p为当前子的根节点,val为要插入的节点权值) TreeNode *insert(TreeNode *p, int val) { if (!p) { // 如果是空节点,说明找到了要插入的位置,新建一个节点 p = (TreeNode *) malloc(sizeof(TreeNode)); p->val = val; p->height = 1; p->left = p->right = NULL; } else if (val < p->val) { // 如果要插入的值小于节点值,向左子递归 p->left = insert(p->left, val); } else { // 否则向右子递归 p->right = insert(p->right, val); } return balance(p); // 插入成功后,再进行平衡操作 } // 寻找以p为根的中权值为val的节点 TreeNode *find(TreeNode *p, int val) { if (!p) { // 如果搜索到空节点,说明未找到 return NULL; } else if (val < p->val) { // 如果要查找的值小于节点值,向左子递归 return find(p->left, val); } else if (val > p->val) { // 如果要查找的值大于节点值,向右子递归 return find(p->right, val); } else { // 否则就是要查找的节点 return p; } } // 查找以p为根的中的最小节点 TreeNode *findMin(TreeNode *p) { while (p->left) { // 不停向左子节点遍历,直到为空 p = p->left; } return p; } // 删除操作(p为当前子的根节点,val为要删除的节点权值) TreeNode *delete(TreeNode *p, int val) { if (!p) { // 如果搜索到空节点,说明未找到,直接返回 return NULL; } else if (val < p->val) { // 如果要删除的值小于节点值,向左子递归 p->left = delete(p->left, val); } else if (val > p->val) { // 如果要删除的值大于节点值,向右子递归 p->right = delete(p->right, val); } else { // 否则就是要删除的节点 if (!p->left && !p->right) { // 如果是叶子节点,直接删除 free(p); return NULL; } else if (p->left && !p->right) { // 如果只有左子节点(或没有子节点),指向左子节点 TreeNode *left = p->left; free(p); return left; } else if (!p->left && p->right) { // 如果只有右子节点,指向右子节点 TreeNode *right = p->right; free(p); return right; } else { // 如果有左、右子节点,用右子节点的最小值代替当前节点,然后再删除右子节点最小值 TreeNode *minRight = findMin(p->right); p->val = minRight->val; p->right = delete(p->right, p->val); } } return balance(p); // 删除成功后,再进行平衡操作 } // 深度优先遍历(先序遍历) void dfs(TreeNode *p) { if (!p) return; printf("%d ", p->val); dfs(p->left); dfs(p->right); } // 主函数 int main() { TreeNode *root = NULL; root = insert(root, 10); root = insert(root, 20); root = insert(root, 30); root = insert(root, 40); root = insert(root, 50); root = delete(root, 30); root = delete(root, 20); root = insert(root, 25); dfs(root); return 0; } ``` 上述代码给出了平衡树节点结构体、获取节点高度、获取节点平衡因子、更新节点高度、左旋、右旋、平衡操作、插入、查找、查找最小节点和删除等常用函数实现。同时,还给出了一个简单的测试代码,用于对平衡树进行测试和验证。 使用平衡树可以大大提高数据插入、删除、查找等操作的效率,所以是实际应用中很常用的数据结构之一。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值