平衡树学习(2)——Splay

2024-02-04 12:08
文章标签 学习 平衡 splay

本文主要是介绍平衡树学习(2)——Splay,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

Splay是一种平衡树,它的代码复杂度和时间复杂度稍弱于Treap,但由于其可以支持区间操作,所以在实战中还是有许多用处。
我们先来看看Splay的定义和基本思路。

伸展树(Splay Tree),也叫分裂树,它能在O(log n)内完成插入、查找和删除操作。在伸展树上的一般操作都基于伸展操作:假设想要对一个二叉查找树执行一系列的查找操作,为了使整个查找时间更小,被查频率高的那些条目就应当经常处于靠近树根的位置。于是想到设计一个简单方法, 在每次查找之后对树进行重构,把被查找的条目搬移到离树根近一些的地方。伸展树应运而生。伸展树是一种自调整形式的二叉查找树,它会沿着从某个节点到树根之间的路径,通过一系列的旋转把这个节点搬移到树根去。
它的优势在于不需要记录用于平衡树的冗余信息。

Ps:我们介绍的Splay是双旋操作Splay。


数组定义:
struct node{int fa,siz,ct,val,son[2];
}T[MAXN];
//fa记录父亲节点,siz记录子树大小,val记录当前节点值,ct表示与当
//前节点值相同的点个数,son[2]记录其左儿子和右儿子。

预先操作:

1.clear:清空点k里的数据

void clear(int k){T[k].siz=T[k].ct=T[k].fa=T[k].val=T[k].son[0]=T[k].son[1]=0;
} 

2.up:统计k的子树大小

void up(int k){T[k].siz=T[L].siz+T[R].siz+T[k].ct;
} 

3.get:查询k是父亲的左儿子还是右儿子

int get(int k){return T[T[k].fa].son[1]==k;
}

旋转:

对于Splay来说,最重要的就是rotate旋转和splay伸展操作。
对于rotate操作,我们可以类比于Treap的旋转。

void rotate(int &k){int fa=T[k].fa,gran=T[fa].fa;int d1=get(k),d2=get(fa);T[T[k].son[d1^1]].fa=fa;T[fa].son[d1]=T[k].son[d1^1];T[k].son[d1^1]=fa;T[k].fa=gran;T[fa].fa=k;if(gran) T[gran].son[d2]=k;;up(fa);up(k);
}

splay操作是Rotate的升级版,该函数将子节点旋转到根,来保持Splay的复杂度。
每次的splay操作,我们都将其旋转到根节点。对于Splay的伸展操作,我们需要进行讨论:
1.当旋转的点中有根节点时,可以直接进行旋转。
这里写图片描述

2.当点x和其父亲,祖父三点共线,则先旋转x的父亲,然后旋转x。
这里写图片描述
3.如果x和其父亲,祖父三点不共线,旋转两次x。

void splay(int k){for(int fa;fa=T[k].fa;rotate(k))if(T[fa].fa) rotate(get(fa)==get(k)?fa:k);root=k; 
}

删除:

删除操作需要进行讨论,如果对于值x存在多个,那直接将x所在节点k的size和ct减1即可。如果不是这样,我们继续讨论,如果k没有儿子,直接clear节点k即可。如果k只有一个儿子,就将它的儿子接上来即可。如果k有两个儿子,就将它的前驱接到它父亲上,然后将k的右儿子接到前驱上,这样就删去了k节点。

int find(int k,int val){   //先求出要删除点的位置,并进行splay伸展if(!k) return 0;if(val==T[k].val){splay(k);return k;}return find(T[k].son[val>T[k].val],val);
}
void delet(int x){int pl=find(root,x);if(!pl) return;if(T[root].ct>1){T[root].ct--,T[root].siz--;return;}if(!T[root].son[1]&&!T[root].son[0]){clear(root),root=0;return;}if(!T[root].son[1]||!T[root].son[0]){int rt=root;root=T[root].son[0]+T[root].son[1];T[root].fa=0;clear(rt);return;}int rt=root;pre(root,x);splay(dist);T[T[rt].son[1]].fa=root;T[root].son[1]=T[rt].son[1];clear(rt);up(root);T[root].fa=0;return;
}

其它操作:

如求取前驱,后继等操作,类似于Treap

void insert(int &k,int val,int pos){   //插入操作if(!k){k=++sz;T[k].ct=T[k].siz=1;T[k].fa=pos;T[k].son[0]=T[k].son[1]=0;T[k].val=val;dist=k;return; }if(val==T[k].val){dist=k;T[k].siz++;T[k].ct++;return;}if(val<T[k].val) insert(T[k].son[0],val,k);if(val>T[k].val) insert(T[k].son[1],val,k);up(k);
}
int searchRANK(int k,int x){  //查询x排在第几if(!k) return 0;if(x<T[k].val) return searchRANK(L,x);if(x==T[k].val) return T[L].siz+1; if(x>T[k].val) return T[L].siz+T[k].ct+searchRANK(R,x);
}
int searchPLACE(int k,int x){    //查询排在x的数是几if(!k) return 0;if(x<=T[L].siz) return searchPLACE(L,x);x-=T[L].siz;if(x<=T[k].ct) return T[k].val;x-=T[k].ct;return searchPLACE(R,x);
}
void pre(int k,int val){   //前驱,dist记录位置if(!k) return;if(val<=T[k].val) pre(L,val);else dist=k,pre(R,val);
}
void ahe(int k,int val){   //后继if(!k) return;if(val>=T[k].val) ahe(R,val);else dist=k,ahe(L,val);
}

合并:

①当S1中的所有元素小于S2(比如S1和S2是刚刚分裂出来的)时,只需要把S1最大的点伸展到根,然后连边即可。
②当S1和S2大小任意时,启发式合并,把小的往大的身上挪。
分裂:(以k为界限,左边小于或等于k,右边大于或等于k)


BZOJ3224模板:
#include<bits/stdc++.h>
#define L T[k].son[0]
#define R T[k].son[1]
#define MAXN 100005
using namespace std;
int read(){char c;int x=0,y=1;while(c=getchar(),(c<'0'||c>'9')&&c!='-');if(c=='-') y=-1;else x=c-'0';while(c=getchar(),c>='0'&&c<='9') x=x*10+c-'0';return x*y;
}
int root,dist,sz,n;
struct node{int fa,siz,ct,val,son[2];
}T[MAXN];
void clear(int k){T[k].siz=T[k].ct=T[k].fa=T[k].val=T[k].son[0]=T[k].son[1]=0;
} 
void up(int k){T[k].siz=T[L].siz+T[R].siz+T[k].ct;
} 
int get(int k){return T[T[k].fa].son[1]==k;
}
void rotate(int &k){int fa=T[k].fa,gran=T[fa].fa;int d1=get(k),d2=get(fa);T[T[k].son[d1^1]].fa=fa;T[fa].son[d1]=T[k].son[d1^1];T[k].son[d1^1]=fa;T[k].fa=gran;T[fa].fa=k;if(gran) T[gran].son[d2]=k;;up(fa);up(k);
}
void splay(int k){for(int fa;fa=T[k].fa;rotate(k))if(T[fa].fa) rotate(get(fa)==get(k)?fa:k);root=k; 
}
void insert(int &k,int val,int pos){if(!k){k=++sz;T[k].ct=T[k].siz=1;T[k].fa=pos;T[k].son[0]=T[k].son[1]=0;T[k].val=val;dist=k;return; }if(val==T[k].val){dist=k;T[k].siz++;T[k].ct++;return;}if(val<T[k].val) insert(T[k].son[0],val,k);if(val>T[k].val) insert(T[k].son[1],val,k);up(k);
}
int searchRANK(int k,int x){if(!k) return 0;if(x<T[k].val) return searchRANK(L,x);if(x==T[k].val) return T[L].siz+1; if(x>T[k].val) return T[L].siz+T[k].ct+searchRANK(R,x);
}
int searchPLACE(int k,int x){if(!k) return 0;if(x<=T[L].siz) return searchPLACE(L,x);x-=T[L].siz;if(x<=T[k].ct) return T[k].val;x-=T[k].ct;return searchPLACE(R,x);
}
void pre(int k,int val){if(!k) return;if(val<=T[k].val) pre(L,val);else dist=k,pre(R,val);
}
void ahe(int k,int val){if(!k) return;if(val>=T[k].val) ahe(R,val);else dist=k,ahe(L,val);
}
int find(int k,int val){if(!k) return 0;if(val==T[k].val){splay(k);return k;}return find(T[k].son[val>T[k].val],val);
}
void delet(int x){int pl=find(root,x);if(!pl) return;if(T[root].ct>1){T[root].ct--,T[root].siz--;return;}if(!T[root].son[1]&&!T[root].son[0]){clear(root),root=0;return;}if(!T[root].son[1]||!T[root].son[0]){int rt=root;root=T[root].son[0]+T[root].son[1];T[root].fa=0;clear(rt);return;}int rt=root;pre(root,x);splay(dist);T[T[rt].son[1]].fa=root;T[root].son[1]=T[rt].son[1];clear(rt);up(root);T[root].fa=0;return;
}
int main()
{n=read();for(int i=1;i<=n;i++){int x=read(),y=read();if(x==1) insert(root,y,0),splay(dist);if(x==2) delet(y);if(x==3) printf("%d\n",searchRANK(root,y));if(x==4) printf("%d\n",searchPLACE(root,y));if(x==5) dist=0,pre(root,y),printf("%d\n",T[dist].val);if(x==6) dist=0,ahe(root,y),printf("%d\n",T[dist].val); }return 0;
}

文献资料:图片转自博客ZigZagK,感谢ZZK大佬

这篇关于平衡树学习(2)——Splay的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/677436

相关文章

Java深度学习库DJL实现Python的NumPy方式

《Java深度学习库DJL实现Python的NumPy方式》本文介绍了DJL库的背景和基本功能,包括NDArray的创建、数学运算、数据获取和设置等,同时,还展示了如何使用NDArray进行数据预处理... 目录1 NDArray 的背景介绍1.1 架构2 JavaDJL使用2.1 安装DJL2.2 基本操

HarmonyOS学习(七)——UI(五)常用布局总结

自适应布局 1.1、线性布局(LinearLayout) 通过线性容器Row和Column实现线性布局。Column容器内的子组件按照垂直方向排列,Row组件中的子组件按照水平方向排列。 属性说明space通过space参数设置主轴上子组件的间距,达到各子组件在排列上的等间距效果alignItems设置子组件在交叉轴上的对齐方式,且在各类尺寸屏幕上表现一致,其中交叉轴为垂直时,取值为Vert

Ilya-AI分享的他在OpenAI学习到的15个提示工程技巧

Ilya(不是本人,claude AI)在社交媒体上分享了他在OpenAI学习到的15个Prompt撰写技巧。 以下是详细的内容: 提示精确化:在编写提示时,力求表达清晰准确。清楚地阐述任务需求和概念定义至关重要。例:不用"分析文本",而用"判断这段话的情感倾向:积极、消极还是中性"。 快速迭代:善于快速连续调整提示。熟练的提示工程师能够灵活地进行多轮优化。例:从"总结文章"到"用

【前端学习】AntV G6-08 深入图形与图形分组、自定义节点、节点动画(下)

【课程链接】 AntV G6:深入图形与图形分组、自定义节点、节点动画(下)_哔哩哔哩_bilibili 本章十吾老师讲解了一个复杂的自定义节点中,应该怎样去计算和绘制图形,如何给一个图形制作不间断的动画,以及在鼠标事件之后产生动画。(有点难,需要好好理解) <!DOCTYPE html><html><head><meta charset="UTF-8"><title>06

学习hash总结

2014/1/29/   最近刚开始学hash,名字很陌生,但是hash的思想却很熟悉,以前早就做过此类的题,但是不知道这就是hash思想而已,说白了hash就是一个映射,往往灵活利用数组的下标来实现算法,hash的作用:1、判重;2、统计次数;

零基础学习Redis(10) -- zset类型命令使用

zset是有序集合,内部除了存储元素外,还会存储一个score,存储在zset中的元素会按照score的大小升序排列,不同元素的score可以重复,score相同的元素会按照元素的字典序排列。 1. zset常用命令 1.1 zadd  zadd key [NX | XX] [GT | LT]   [CH] [INCR] score member [score member ...]

【机器学习】高斯过程的基本概念和应用领域以及在python中的实例

引言 高斯过程(Gaussian Process,简称GP)是一种概率模型,用于描述一组随机变量的联合概率分布,其中任何一个有限维度的子集都具有高斯分布 文章目录 引言一、高斯过程1.1 基本定义1.1.1 随机过程1.1.2 高斯分布 1.2 高斯过程的特性1.2.1 联合高斯性1.2.2 均值函数1.2.3 协方差函数(或核函数) 1.3 核函数1.4 高斯过程回归(Gauss

【学习笔记】 陈强-机器学习-Python-Ch15 人工神经网络(1)sklearn

系列文章目录 监督学习:参数方法 【学习笔记】 陈强-机器学习-Python-Ch4 线性回归 【学习笔记】 陈强-机器学习-Python-Ch5 逻辑回归 【课后题练习】 陈强-机器学习-Python-Ch5 逻辑回归(SAheart.csv) 【学习笔记】 陈强-机器学习-Python-Ch6 多项逻辑回归 【学习笔记 及 课后题练习】 陈强-机器学习-Python-Ch7 判别分析 【学

系统架构师考试学习笔记第三篇——架构设计高级知识(20)通信系统架构设计理论与实践

本章知识考点:         第20课时主要学习通信系统架构设计的理论和工作中的实践。根据新版考试大纲,本课时知识点会涉及案例分析题(25分),而在历年考试中,案例题对该部分内容的考查并不多,虽在综合知识选择题目中经常考查,但分值也不高。本课时内容侧重于对知识点的记忆和理解,按照以往的出题规律,通信系统架构设计基础知识点多来源于教材内的基础网络设备、网络架构和教材外最新时事热点技术。本课时知识

线性代数|机器学习-P36在图中找聚类

文章目录 1. 常见图结构2. 谱聚类 感觉后面几节课的内容跨越太大,需要补充太多的知识点,教授讲得内容跨越较大,一般一节课的内容是书本上的一章节内容,所以看视频比较吃力,需要先预习课本内容后才能够很好的理解教授讲解的知识点。 1. 常见图结构 假设我们有如下图结构: Adjacency Matrix:行和列表示的是节点的位置,A[i,j]表示的第 i 个节点和第 j 个