模板——Splay
$Splay$
#include <bits/stdc++.h>
#define inf (int)1e9
using namespace std;
const int N=1e5+100;
int n,tot,root,val[N],sz[N],son[N][2];
int fa[N],sf[N],re[N];
int newnode()
{
return tot++;
}
void connect(int x,int y,int dir)
{
son[y][dir]=x;
fa[x]=y;
sf[x]=dir;
}
void pushup(int x)
{
sz[x]=sz[son[x][0]]+sz[son[x][1]]+re[x];
}
void clear()
{
val[0]=son[0][0]=son[0][1]=sz[0]=re[0]=fa[0]=sf[0]=0;
}
void rotate(int x)
{
int f,gf,xd,fd,s;
f=fa[x];gf=fa[f];
xd=sf[x];fd=sf[f];
s=son[x][xd^1];
connect(x,gf,fd);connect(f,x,xd^1);connect(s,f,xd);
clear();
pushup(f);pushup(x);
}
void splay(int x,int y)
{
while (fa[x]!=y)
{
if (fa[fa[x]]==y)
rotate(x);
else
if (sf[fa[x]]==sf[x])
{
rotate(fa[x]);
rotate(x);
}
else
{
rotate(x);
rotate(x);
}
}
if (y==0)
root=x;
}
void find(int v)
{
int cur=root;
while (son[cur][v>val[cur]] && val[cur]!=v)
cur=son[cur][v>val[cur]];
splay(cur,0);
}
int per(int v)
{
find(v);
if (val[root]<v) return root;
int cur=son[root][0];
while (son[cur][1]) cur=son[cur][1];
return cur;
}
int suc(int v)
{
find(v);
if (val[root]>v) return root;
int cur=son[root][1];
while (son[cur][0]) cur=son[cur][0];
return cur;
}
void insert(int v)
{ if (root==0)
{
int x=newnode();
val[x]=v;sz[x]=1;
root=x;
return;
}
int cur=root;
while (son[cur][v>val[cur]] && val[cur]!=v)
cur=son[cur][v>val[cur]];
if (val[cur]==v)
{
re[cur]++;sz[cur]++;
splay(cur,0);
return;
}
int x=newnode();
val[x]=v;re[x]=sz[x]=1;
connect(x,cur,v>val[cur]);
pushup(cur);
splay(x,0);
}
void del(int v)
{
int p,s;
p=per(v);s=suc(v);
splay(p,0);
splay(s,p);
re[son[s][0]]--;sz[son[s][0]]--;
if (re[son[s][0]]==0)
son[s][0]=0;
pushup(s);pushup(p);
}
int rk(int v)
{
find(v);
return sz[son[root][0]];
}
int kth(int x,int k)
{
if (k>sz[son[x][0]]+re[x])
return kth(son[x][1],k-sz[son[x][0]]-re[x]);
if (k<=sz[son[x][0]])
return kth(son[x][0],k);
return x;
}
int main()
{
tot=1;
insert(inf);insert(-inf);
scanf("%d",&n);
for (int i=1;i<=n;i++)
{
int op,x;
scanf("%d%d",&op,&x);
if (op==1) insert(x);
if (op==2) del(x);
if (op==3) printf("%d\n",rk(x));
if (op==4) printf("%d\n",val[kth(root,x+1)]);
if (op==5) printf("%d\n",val[per(x)]);
if (op==6) printf("%d\n",val[suc(x)]);
}
}
维护数列
#include <bits/stdc++.h>
#define inf (int)1e9
using namespace std;
const int N=500100;
int n,m,tot,root,a[N];
struct node
{
int sz,val,sum,res,lx,rx,tx,vc;
int son[2],fa,sf;
}sh[N+100];
queue <int> q;
int newnode()
{
int x;
if (tot>=N)
{
x=q.front();
q.pop();
}
else
x=tot++;
sh[x].sz=sh[x].val=sh[x].sum=sh[x].res=sh[x].lx=sh[x].rx=sh[x].tx=sh[x].vc=0;
sh[x].son[0]=sh[x].son[1]=sh[x].sf=sh[x].fa=0;
sh[x].vc=inf;
return x;
}
void clear(int x)
{
q.push(x);
if (sh[x].son[0]) clear(sh[x].son[0]);
if (sh[x].son[1]) clear(sh[x].son[1]);
}
void connect(int x,int y,int dir)//x->y
{
sh[y].son[dir]=x;
sh[x].fa=y;
sh[x].sf=dir;
}
void pushdown(int x)
{
int ls,rs;
ls=sh[x].son[0];rs=sh[x].son[1];
if (sh[x].res==1)
{
if (ls)
{
sh[sh[ls].son[0]].sf^=1;sh[sh[ls].son[1]].sf^=1;
swap(sh[ls].son[0],sh[ls].son[1]);
swap(sh[sh[ls].son[0]].lx,sh[sh[ls].son[0]].rx);
swap(sh[sh[ls].son[1]].lx,sh[sh[ls].son[1]].rx);
sh[ls].res^=1;
}
if (rs)
{
sh[sh[rs].son[0]].sf^=1;sh[sh[rs].son[1]].sf^=1;
swap(sh[rs].son[0],sh[rs].son[1]);
swap(sh[sh[rs].son[0]].lx,sh[sh[rs].son[0]].rx);
swap(sh[sh[rs].son[1]].lx,sh[sh[rs].son[1]].rx);
sh[rs].res^=1;
}
sh[x].res=0;
}
if (sh[x].vc!=inf)
{
if (ls)
{
sh[ls].val=sh[ls].vc=sh[x].vc;
sh[ls].sum=sh[ls].val*sh[ls].sz;
sh[ls].tx=max(sh[ls].val,sh[ls].val*sh[ls].sz);
sh[ls].lx=sh[ls].rx=max(0,sh[ls].val*sh[ls].sz);
}
if (rs)
{
sh[rs].val=sh[rs].vc=sh[x].vc;
sh[rs].sum=sh[rs].val*sh[rs].sz;
sh[rs].tx=max(sh[rs].val,sh[rs].val*sh[rs].sz);
sh[rs].lx=sh[rs].rx=max(0,sh[rs].val*sh[rs].sz);
}
sh[x].vc=inf;
}
}
void pushup(int x)
{
int ls,rs;
ls=sh[x].son[0];rs=sh[x].son[1];
sh[x].sz=sh[ls].sz+sh[rs].sz+1;
sh[x].sum=sh[ls].sum+sh[rs].sum+sh[x].val;
sh[x].lx=max(sh[ls].lx,sh[ls].sum+sh[rs].lx+sh[x].val);
sh[x].rx=max(sh[rs].rx,sh[rs].sum+sh[ls].rx+sh[x].val);
sh[x].tx=max(sh[ls].tx,max(sh[rs].tx,sh[ls].rx+sh[rs].lx+sh[x].val));
}
void rotate(int x)
{
int fa,gf,son,xd,fd;
fa=sh[x].fa;gf=sh[sh[x].fa].fa;
xd=sh[x].sf;fd=sh[fa].sf;
son=sh[x].son[xd^1];
connect(x,gf,fd);connect(fa,x,xd^1);connect(son,fa,xd);
sh[0].fa=sh[0].son[0]=sh[0].son[1]=sh[0].sf=0;
pushup(fa);pushup(x);
}
void splay(int x,int y)
{
while (sh[x].fa!=y)
{
if (sh[sh[x].fa].fa==y)
rotate(x);
else
if (sh[x].sf==sh[sh[x].fa].sf)
{
rotate(sh[x].fa);
rotate(x);
}
else
{
rotate(x);
rotate(x);
}
}
if (y==0)
root=x;
}
int find(int x,int k)
{
pushdown(x);
if (k>sh[sh[x].son[0]].sz+1)
return find(sh[x].son[1],k-sh[sh[x].son[0]].sz-1);
if (k<=sh[sh[x].son[0]].sz)
return find(sh[x].son[0],k);
return x;
}
int build(int l,int r,int father,int dir)
{
int mid=(l+r)>>1;
int x=newnode();
sh[x].val=a[mid];
sh[x].fa=father;sh[x].sf=dir;
if (l==r)
{
sh[x].sz=1;
sh[x].sum=sh[x].tx=a[mid];
sh[x].lx=sh[x].rx=max(0,a[mid]);
return x;
}
if (l<=mid-1)
sh[x].son[0]=build(l,mid-1,x,0);
if (r>=mid+1)
sh[x].son[1]=build(mid+1,r,x,1);
pushup(x);
return x;
}
int split(int l,int r)
{
int per,suc;
per=find(root,l);suc=find(root,r+2);
splay(per,0);
splay(suc,per);
return sh[suc].son[0];
}
int main()
{
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++)
scanf("%d",&a[i]);
tot=1;sh[0].tx=a[0]=a[n+1]=-inf;
root=build(0,n+1,0,0);
while (m--)
{
char ch[20];
scanf("%s",ch);
if (ch[0]=='I')
{
int pos,tot;
scanf("%d%d",&pos,&tot);
for (int i=1;i<=tot;i++)
scanf("%d",&a[i]);
int x=build(1,tot,0,0);
int A,B;
A=find(root,pos+1);B=find(root,pos+2);
splay(A,0);splay(B,A);
connect(x,B,0);
pushup(B);pushup(A);
}
if (ch[0]=='D')
{
int pos,tot;
scanf("%d%d",&pos,&tot);
int per,suc;
per=find(root,pos);suc=find(root,pos+tot+1);
splay(per,0);
splay(suc,per);
clear(sh[suc].son[0]);
sh[suc].son[0]=0;
pushup(suc);pushup(per);
}
if (ch[0]=='M' && ch[2]=='K')
{
int pos,tot,c;
scanf("%d%d%d",&pos,&tot,&c);
int x=split(pos,pos+tot-1);
sh[x].val=sh[x].vc=c;
sh[x].sum=c*sh[x].sz;
sh[x].tx=max(sh[x].val,c*sh[x].sz);
sh[x].lx=sh[x].rx=max(0,c*sh[x].sz);
pushup(sh[x].fa);pushup(sh[sh[x].fa].fa);
}
if (ch[0]=='R')
{
int pos,tot;
scanf("%d%d",&pos,&tot);
int x=split(pos,pos+tot-1);
sh[x].res^=1;
int ls,rs;
ls=sh[x].son[0];rs=sh[x].son[1];
sh[ls].sf^=1;sh[rs].sf^=1;
swap(sh[x].son[0],sh[x].son[1]);
swap(sh[ls].lx,sh[ls].rx);
swap(sh[rs].lx,sh[rs].rx);
pushup(x);pushup(sh[x].fa);pushup(sh[sh[x].fa].fa);
}
if (ch[0]=='G')
{
int pos,tot;
scanf("%d%d",&pos,&tot);
int x=split(pos,pos+tot-1);
printf("%d\n",sh[x].sum);
}
if (ch[0]=='M' && ch[2]=='X')
{
printf("%d\n",sh[root].tx);
}
}
}
模板——Splay的更多相关文章
- 算法模板——splay区间反转 2
实现功能:同splay区间反转 1(基于BZOJ3223 文艺平衡树) 这次改用了一个全新的模板(HansBug:琢磨了我大半天啊有木有),大大简化了程序,同时对于splay的功能也有所完善 这里面没 ...
- 算法模板——splay区间反转 1
实现的功能:将序列区间反转,并维护 详见BZOJ3223 var i,j,k,l,m,n,head,a1,a2:longint; s1:ansistring; a,b,c,d,fat,lef,rig: ...
- 【luogu P3369 普通平衡树（Treap/SBT）】 模板 Splay
题目链接:https://www.luogu.org/problemnew/show/P3369 #include <cstdio> #include <algorithm> ...
- 洛谷 P3391 模板Splay
#include<bits/stdc++.h> using namespace std; #define maxn 200000 int read() { ,w=; ;ch=getchar ...
- [模板] Splay
欠了好久的Splay,以后就它了. 默写真不容易,过几天估计就忘了.. 整个Splay真的精妙,不拖泥带水那种.. 前驱后继之所以不能用rk转到根,是因为这个数不一定存在.. kth中<=老忘记 ...
- 模板—splay
#include<iostream> #include<cstdio> #define cin(x) scanf("%d",&x) using na ...
- Splay 伸展树
废话不说,有篇论文可供参考:杨思雨:<伸展树的基本操作与应用> Splay的好处可以快速分裂和合并. ===============================14.07.26更新== ...
- [NOI2003][bzoj1507] 文本编辑器 editor [splay]
其实看明白了就是一道水题 毕竟模板 splay敲一发,插入一个串的时候先把它构建成一棵平衡树,再挂到原来的splay上面去即可 没别的了,上代码 #include<iostream> #i ...
- splay最终模板
来自wjmzbmr的splay模板 #include<cstdio> #include<iostream> #include<algorithm> using na ...
随机推荐
- [VBA原创源代码] excelhome 汇总多工作表花名册
生病了,一点一滴的积累,慢慢康复,今年十月,我就 2 周岁了. 以下代码完成了excelhome中留的作业 http://club.excelhome.net/forum.php?mod=viewth ...
- Java知识系统回顾整理01基础03变量02基本变量类型
一.变量类型分类 一个变量的类型,决定了该变量可以包含什么样的值. Java中有八种基本类型,都是Java语言预先定义好的,并且是关键字. 这八种基本类型分别是: 整型 (4种) 字符型 (1种) ...
- 【随笔---转载】xip.io
http://xip.io/ wildcard DNS for everyone 今天看到一个老外发布的DNS小工具,XIP.IO.功能十分简单,就是将foo.bar.10.0.0.1.xip.io ...
- python文档下载
网址记录:https://docs.python.org/3.6/
- Noip 2013 真题练习
Day1 T1 转圈游戏 Link 一句话题意: 让你求 \({x + m \times 10^k} \bmod n\) 的结果. 直接套上快速幂的板子. code #include<iostr ...
- IDEA中创建父子工程与maven打包Springboot聚合工程报错程序包不存在问题处理
公司新项目需使用java技术栈,便使用IDEA搭建了一个多SpringBoot项目的聚合工程,因为初次使用,遇到了很多问题,maven打包时各种报错,在网上查了好多终于解决了,为巩固记忆,特作此记录. ...
- redis协议规范
好多年前看过redis的代码,那个时候还是2.6的版本,集群和哨兵还没加入正式代码,这几年redis发展的好快.简略翻译一篇文章redis的https://redis.io/topics/protoc ...
- springboot2.2.2企业级项目整合redis与redis 工具类大全
1.springboot2.2.2整合redis教程很多,为此编写了比较完整的redis工具类,符合企业级开发使用的工具类 2.springboot与redis maven相关的依赖 <depe ...
- 制作西北地区地图数据并maskout
1.从全国地图数据中选中西北5省:打开bou2_4p.shp文件添加相应的图层(中国各省的行政区域),选中工具栏中的"通过矩形选择要素"工具,用鼠标点击选择要输出的图元,按住ctr ...
- react中 受控组件和 非受控组件 浅析
一 受控组件 顾名思义,受控 也就是能够被控制,简而言之也就是 该组件ui的显示或者内部state逻辑的变化依赖外部的 props的传入. 二 非受控组件 顾名思义,非受控,也就是内部的视图变化,st ...