区间树Splay——[NOI2005]维护数列

时间:2022-01-13 16:19:58

无指针Splay超详细讲解

区间树这玩意真TM玄学。

学这东西你必须要拥有的

1.通过【模板】文艺平衡树(Splay)【模板】普通平衡树GSS3 - Can you answer these queries III

2.学会Splay,学会求最大子段和并知道怎么维护信息和下传标记,及会有区间修的最大子段和

3.多年的编程技巧,以及一颗写数据结构的良好心态

4.攒够两个月的肝,这很重要!

如果你不会上面东西的解决方法

1.看以下博客Splay入门解析文艺平衡树Splay题解GSS系列题解——最大子段和系列

2.看上面

3.别管,瞎逼的

4.好好养生,如果不够肝的话千万别写这道题

那么现在就可以开始了

既然你已经会了上面的前置技能,那么我们就可以开始分步解决这道题了。

先给出我们需要存的全部信息:

struct kkk{
int ch[2]; //左右儿子
int size; //子树大小
int fa; //父亲
int tag; //赋值标记
int val; //权值
int rev; //翻转标记
int sum; //区间权值和
int left; //左区间,指区间最大前缀和
int right; //右区间,指区间最大后缀和
int middle; //中区间,指区间最大子段和
void clear(){ch[0]=ch[1]=fa=rev=0;tag=TAGNONE;} //清空节点信息
}tree[maxn];

存的东西很多,大家务必要理解清楚每一个信息所表达的含义。

区间树Splay介绍

做过“普通平衡树”的都知道,在“普通平衡树”里,Splay是按照权值来排序的,所以能维护数的关系。那么现在到了维护区间上的操作了,也就不能按权值来排序了。

区间树,我们按照的是序列中的编号来排序。

我们可以发现,序列中的第k个点,在Splay中也是第k大的。(按编号排序嘛

所以我们想要查找序列中第k个位置,就直接找Splay中的第k大就可以了。

所以“普通平衡树”里的Splay操作,rotate操作和kth操作都是可以直接照搬的(一样的,只是维护编号而已

那么我们怎么在Splay中找到一个区间[x,y]呢?

我们可以考虑Splay的性质,将 x Splay上根,再将 y Splay上到x的右节点,那么我们得出的 y 的左子树就是我们要的[x,y]区间。

之后我们想对这个区间做什么就可以直接对那颗子树做了。

上面就是区间树的一些介绍

代码中的一些宏定义

#define TAGNONE 10000001				//没有赋值tag的标志
#define L(node) (tree[node].ch[0]) //替左儿子
#define R(node) (tree[node].ch[1]) //替右儿子
#define F(node) (tree[node].fa) //替父亲
#define V(node) (tree[node].val) //替权值
#define S(node) (tree[node].size) //替子树大小
#define compare(node,x) (tree[node].val<x) //比较node是权值x的左儿子还是右儿子

操作剖析

1.基本操作 Splay,rotate,kth

这个就不用怎么说了吧,大家在做平衡树Splay都写过的啦!

2.将指定区间找出来 split操作

和上面讲的区间树一样,先找到区间[l,r]的kth,计l的kth为xr的kth为y

然后Splay(x,0);Splay(y,x); (直接上代码解释)

最后返回y左儿子就是指定区间

代码:

int split(int k,int len){	//找到那个区间的位置
int x=kth(k),y=kth(k+len+1);
Splay(x,0);Splay(y,x);
return L(y);
}

3.建一颗平衡的Splay,build操作

一开始我们要构造一颗有初始信息的Splay,一个一个insert显然很慢,所以我们写一个build,可以将一段序列建成一颗平衡的Splay的操作。

其实写起来和线段树差不多,注意是以编号排序来建树。

void New(int node,int x){                       //新建节点
tree[node].middle=tree[node].sum=x; //赋值信息
tree[node].tag=TAGNONE;tree[node].rev=0; //标记初始化
tree[node].left=tree[node].right=max(x,0); //区间赋值
tree[node].size=1; //大小赋值
}
void build(int begin,int end,int fa){ //建树
int mid=(begin+end)>>1;int node=id[mid],pre=id[fa];
if(begin==end) //到达底部
New(node,a[begin]); //新建一个节点
if(begin<mid)build(begin,mid-1,mid); //建左子树
if(mid<end)build(mid+1,end,mid); //建右子树
tree[node].val=a[mid];tree[node].fa=pre;tree[node].tag=TAGNONE; //基本信息赋值
pushup(node); //维护信息
tree[pre].ch[mid>=fa]=node;
}

4.插入操作 insert

这里题目要求的是在x位置后插入一段长为len的序列

如果我们还是一个一个插入,仍然很慢,所以我们可以直接把插入的序列build成一颗平衡的子树,最后直接在x后插入建成的子树就可以了。

void insert(int k,int len){			//插入区间
for(int i=1;i<=len;i++)scanf("%d",&a[i]); //输入区间
for(int i=1;i<=len;i++)
id[i]=rublish(); //从垃圾桶里找一个编号
build(1,len,0); //将输入的区间建成一个完全二叉树
int z=id[(1+len)>>1];
int x=kth(k+1),y=kth(k+2); //找到要插入的位置
Splay(x,0);Splay(y,x);
tree[z].fa=y; tree[y].ch[0]=z; //将新建的子树插入树中
pushup(y);pushup(x); //维护信息
}

5.删除操作 eraser

这个就更简单了,直接找到那个区间,然后让那个子树的父亲将左儿子清为0就可以了。

但是,为了节省空间,我们加入了一个垃圾回收的操作,就是将删除的节点重新利用起来,以节省空间

所以我们还要遍历一遍子树将那颗子树的节点扔进垃圾桶里

int rublish(){				//垃圾回收
if(top==0)return ++cnt;
int node=rub[top--];
return node;
}
void remove(int node){ //将一个子树清空
if(L(node))remove(L(node)); //继续清空左子树
if(R(node))remove(R(node)); //继续清空右子树
rub[++top]=node; tree[node].clear(); //清空并仍进垃圾桶,定义里有
}
void eraser(int x,int len){ //删除区间
int node=split(x,len),y=F(node);//找到该区间
remove(node);tree[y].ch[0]=0; //删除该区间,子树清空
pushup(y);pushup(F(y)); //维护信息
}

6.修改操作 update

一样的,先找到指定区间的子树,然后直接修改信息,打上赋值标记

void change_val(int node,int val){          //更新点值
if(!node)return ; //空节点返回
tree[node].tag=tree[node].val=val; //打赋值标记,更新权值
tree[node].sum=val*tree[node].size; //更新区间权值和
tree[node].left=tree[node].right=max(tree[node].sum,0); //左右区间更新
tree[node].middle=max(tree[node].sum,val); //最大子段和更新
}
void update(int x,int tot,int val){ //更新区间的指
int node=split(x,tot),y=F(node); //找到该区间
change_val(node,val); //更新该区间
pushup(y);pushup(F(y)); //维护信息
}

7.翻转操作 reverse

一样的,先找到指定的区间的子树,然后直接翻转,打上翻转标记

void change_rev(int node){                  //更新翻转
swap(tree[node].ch[0],tree[node].ch[1]);//交换左右儿子
swap(tree[node].left,tree[node].right); //交换左右区间
tree[node].rev^=1; //打翻转标记
}
void reverse(int x,int len){ //翻转区间
int node=split(x,len),y=F(node);//找到该区间
if(tree[node].tag!=TAGNONE)return ; //如果已经有赋值标记就不用管了
change_rev(node); //翻转该区间
pushup(y);pushup(F(y)); //维护信息
}

8.求和操作 query

这就更简单了,找到指定区间的子树,然后直接输出那颗子树的sum就OK了

void query(int x,int len){	//查询区间权值和
int node=split(x,len); //找到该区间
printf("%d\n",tree[node].sum); //输出答案
}

9.求最大子段和

直接输出root的middle最大子段和

printf("%d\n",tree[root].middle);

难点

10.维护信息和下传标记

这玩意是真毒瘤,不过主要还是难在最大子段和上面,只要我们能理解“GSS3”中的求法,其实也很简单。

维护信息,所以除了这个最大子段和之外,好像还挺简单的。最大子段和那几个更新方法这里就不讲了,不知道可以看上面的博客。

void pushup(int node){					//维护信息
kkk &x=tree[L(node)],&y=tree[R(node)];int val=tree[node].val; //实质是将左右儿子合并,x代替左儿子,y代替右儿子
kkk &res=tree[node]; //res代替tree[node]
res.sum=x.sum+y.sum+val; res.size=x.size+y.size+1; //权值和更新,子树大小更新
res.middle=max(max(x.middle,y.middle),x.right+y.left+val); //最大子段和更新
res.left=max(x.left,x.sum+y.left+val); //区间最大前缀和更新
res.right=max(y.right,y.sum+x.right+val); //区间最大后缀和更新
}

下传标记,这本来是比较得毒瘤,我们要先更新赋值操作,左右儿子有很多信息需要更新,其中就有tag,sum,left,right和middle,更新起来十分的繁琐。但是在之前的赋值操作update中,我们引入了一个叫change_val的函数,所以这里,我们可以直接调用那个函数。于是代码就被减短了很多。

最后将tag标记为TAGNONE就OK了

然后要更新翻转操作,一样的,在之前翻转操作revrese中,我们引入了一个叫change_rev的函数,所以这里,我们还是可以直接调用。于是代码又被减了……

最后将rev标记为0就OK了。

代码:

void pushdown(int node){				//标记下传
if(tree[node].tag!=TAGNONE){ //判断有没有赋值标记
change_val(L(node),tree[node].tag); //更新左儿子
change_val(R(node),tree[node].tag); //更新右儿子
tree[node].tag=TAGNONE; //除去标记
}
if(tree[node].rev){ //判断有没有翻转标记
change_rev(L(node)); //更新左儿子
change_rev(R(node)); //更新右儿子
tree[node].rev=0; //除去标记
}
}

看,多简短

11.主函数

注意边界!注意边界!注意边界! 主要的事情说三遍!

其他就没什么了,都是输入嘛。

总代码

#include<bits/stdc++.h>
#define TAGNONE 10000001
#define maxn 1000010
#define inf 100000001
#define L(node) (tree[node].ch[0]) //替左儿子
#define R(node) (tree[node].ch[1]) //替右儿子
#define F(node) (tree[node].fa) //替父亲
#define V(node) (tree[node].val) //替权值
#define S(node) (tree[node].size) //替子树大小
#define compare(node,x) (tree[node].val<x) //比较node是权值x的左儿子还是右儿子
using namespace std;
int root,cnt,a[maxn],id[maxn],rub[maxn],top,n,m;
struct kkk{
int ch[2]; //左右儿子
int size; //子树大小
int fa; //父亲
int tag; //赋值标记
int val; //权值
int rev; //翻转标记
int sum; //区间权值和
int left; //左区间,指区间最大前缀和
int right; //右区间,指区间最大后缀和
int middle; //中区间,指区间最大子段和
void clear(){ch[0]=ch[1]=fa=rev=0;tag=TAGNONE;} //清空节点信息
}tree[maxn];
int rublish(){ //垃圾回收
if(top==0)return ++cnt;
int node=rub[top--];
return node;
}
void change_val(int node,int val){ //更新点值
if(!node)return ; //空节点返回
tree[node].tag=tree[node].val=val; //打赋值标记,更新权值
tree[node].sum=val*tree[node].size; //更新区间权值和
tree[node].left=tree[node].right=max(tree[node].sum,0); //左右区间更新
tree[node].middle=max(tree[node].sum,val); //最大子段和更新
}
void change_rev(int node){ //更新翻转
swap(tree[node].ch[0],tree[node].ch[1]);//交换左右儿子
swap(tree[node].left,tree[node].right); //交换左右区间
tree[node].rev^=1; //打翻转标记
}
void pushup(int node){ //维护信息
kkk &x=tree[L(node)],&y=tree[R(node)];int val=tree[node].val; //实质是将左右儿子合并,x代替左儿子,y代替右儿子
kkk &res=tree[node]; //res代替tree[node]
res.sum=x.sum+y.sum+val;res.size=x.size+y.size+1; //权值和更新,子树大小更新
res.middle=max(max(x.middle,y.middle),x.right+y.left+val); //最大子段和更新
res.left=max(x.left,x.sum+y.left+val); //区间最大前缀和更新
res.right=max(y.right,y.sum+x.right+val); //区间最大后缀和更新
}
void pushdown(int node){ //标记下传
if(tree[node].tag!=TAGNONE){ //判断有没有赋值标记
change_val(L(node),tree[node].tag); //更新左儿子
change_val(R(node),tree[node].tag); //更新右儿子
tree[node].tag=TAGNONE; //除去标记
}
if(tree[node].rev){ //判断有没有翻转标记
change_rev(L(node)); //更新左儿子
change_rev(R(node)); //更新右儿子
tree[node].rev=0; //除去标记
}
}
void rotate(int node){ //rotate 模板
int fa=F(node);
int gfa=F(fa);
int z=tree[fa].ch[1]==node;
tree[gfa].ch[tree[gfa].ch[1]==fa]=node; tree[node].fa=gfa;
tree[fa].ch[z]=tree[node].ch[z^1];tree[tree[node].ch[z^1]].fa=fa;
tree[node].ch[z^1]=fa;tree[fa].fa=node;
pushup(fa); pushup(node);
}
void Splay(int node,int goal){ //Splay 模板
while(tree[node].fa!=goal){
int fa=F(node);
int gfa=F(fa);
if(gfa!=goal)
(compare(fa,tree[node].val))!=(compare(gfa,tree[fa].val))
?rotate(node) : rotate(fa);
rotate(node);
}
if(!goal)root=node;
}
void New(int node,int x){ //新建节点
tree[node].middle=tree[node].sum=x; //赋值信息
tree[node].tag=TAGNONE;tree[node].rev=0; //标记初始化
tree[node].left=tree[node].right=max(x,0); //区间赋值
tree[node].size=1; //大小赋值
}
void build(int begin,int end,int fa){ //建树
int mid=(begin+end)>>1;int node=id[mid],pre=id[fa];
if(begin==end) //到达底部
New(node,a[begin]); //新建一个节点
if(begin<mid)build(begin,mid-1,mid); //建左子树
if(mid<end)build(mid+1,end,mid); //建右子树
tree[node].val=a[mid];tree[node].fa=pre;tree[node].tag=TAGNONE; //基本信息赋值
pushup(node); //维护信息
tree[pre].ch[mid>=fa]=node;
}
int kth(int x){ //kth模板
int node=root;
while(1){
pushdown(node);
if(tree[L(node)].size>=x)node=L(node);
else
if(tree[L(node)].size+1==x)return node;
else x-=tree[L(node)].size+1,node=R(node);
}
}
void remove(int node){ //将一个子树清空
if(L(node))remove(L(node)); //继续清空左子树
if(R(node))remove(R(node)); //继续清空右子树
rub[++top]=node; tree[node].clear(); //清空并仍进垃圾桶
}
int split(int k,int len){ //找到那个区间的位置
int x=kth(k),y=kth(k+len+1);
Splay(x,0);Splay(y,x);
return L(y);
}
void query(int x,int len){ //查询区间权值和
int node=split(x,len); //找到该区间
printf("%d\n",tree[node].sum); //输出答案
}
void update(int x,int tot,int val){ //更新区间的指
int node=split(x,tot),y=F(node); //找到该区间
change_val(node,val); //更新该区间
pushup(y);pushup(F(y)); //维护信息
}
void rever(int x,int len){ //翻转区间
int node=split(x,len),y=F(node);//找到该区间
if(tree[node].tag!=TAGNONE)return ; //如果已经有赋值标记就不用管了
change_rev(node); //翻转该区间
pushup(y);pushup(F(y)); //维护信息
}
void eraser(int x,int len){ //删除区间
int node=split(x,len),y=F(node);//找到该区间
remove(node);tree[y].ch[0]=0; //删除该区间,子树清空
pushup(y);pushup(F(y)); //维护信息
}
void insert(int k,int len){ //插入区间
for(int i=1;i<=len;i++)scanf("%d",&a[i]); //输入区间
for(int i=1;i<=len;i++)
id[i]=rublish();
build(1,len,0); //将输入的区间建成一个完全二叉树
int z=id[(1+len)>>1];
int x=kth(k+1),y=kth(k+2); //找到要插入的位置
Splay(x,0);Splay(y,x);
tree[z].fa=y; tree[y].ch[0]=z; //将新建的子树插入树中
pushup(y);pushup(x); //维护信息
}
int main(){
scanf("%d%d",&n,&m);
tree[0].middle=a[1]=a[n+2]=-inf; //边界
for(int i=1;i<=n;i++)scanf("%d",&a[i+1]); //输入
for(int i=1;i<=n+2;i++)id[i]=i;
build(1,n+2,0); //建成一颗Splay
root=(n+3)>>1;cnt=n+2; //指根,更新点数
for(int i=1;i<=m;i++){
string s; int x,len,y;
cin>>s;
if(s!="MAX-SUM")scanf("%d%d",&x,&len);
else printf("%d\n",tree[root].middle);
if(s=="INSERT")insert(x,len);
if(s=="DELETE")eraser(x,len);
if(s=="MAKE-SAME")
scanf("%d",&y),update(x,len,y);
if(s=="REVERSE")rever(x,len);
if(s=="GET-SUM")query(x,len);
}
}

后记

学习时有参考I_AM_HelloWord大佬的题解。所以有的地方和他的代码很像。

希望大家都能掌握区间树QwQ。

谢谢观赏