「模板」Splay
代码说明
对于一些变量进行说明：
| 变量名 | 说明 |
|---|---|
rt |
树根 |
ff[u] |
点 \(u\) 的父节点，特别地， ff[rt]=0 |
ch[u][0|1] |
点 \(u\) 的 左/右儿子 |
siz[u] |
点 \(u\) 及其子树大小 |
val[u] |
点 \(u\) 对应的值 |
recy[u] |
点 \(u\) 对应的 val[u] 出现的次数 |
代码
#include<cstdio>
#include<vector>
using namespace std;
#define rep(i,__l,__r) for(signed i=__l,i##_end_=__r;i<=i##_end_;++i)
#define fep(i,__l,__r) for(signed i=__l,i##_end_=__r;i>=i##_end_;--i)
#define writc(a,b) fwrit(a),putchar(b)
#define mp(a,b) make_pair(a,b)
#define ft first
#define sd second
#define LL long long
#define ull unsigned long long
#define uint unsigned int
#define pii pair< int,int >
#define Endl putchar('\n')
// #define FILEOI
// #define int long long
// #define int unsigned
#ifdef FILEOI
# define MAXBUFFERSIZE 500000
inline char fgetc(){
static char buf[MAXBUFFERSIZE+5],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,MAXBUFFERSIZE,stdin),p1==p2)?EOF:*p1++;
}
# undef MAXBUFFERSIZE
# define cg (c=fgetc())
#else
# define cg (c=getchar())
#endif
template<class T>inline void qread(T& x){
char c;bool f=0;
while(cg<'0'||'9'<c)f|=(c=='-');
for(x=(c^48);'0'<=cg&&c<='9';x=(x<<1)+(x<<3)+(c^48));
if(f)x=-x;
}
inline int qread(){
int x=0;char c;bool f=0;
while(cg<'0'||'9'<c)f|=(c=='-');
for(x=(c^48);'0'<=cg&&c<='9';x=(x<<1)+(x<<3)+(c^48));
return f?-x:x;
}
template<class T,class... Args>inline void qread(T& x,Args&... args){qread(x),qread(args...);}
template<class T>inline T Max(const T x,const T y){return x>y?x:y;}
template<class T>inline T Min(const T x,const T y){return x<y?x:y;}
template<class T>inline T fab(const T x){return x>0?x:-x;}
inline int gcd(const int a,const int b){return b?gcd(b,a%b):a;}
inline void getInv(int inv[],const int lim,const int MOD){
inv[0]=inv[1]=1;for(int i=2;i<=lim;++i)inv[i]=1ll*inv[MOD%i]*(MOD-MOD/i)%MOD;
}
template<class T>void fwrit(const T x){
if(x<0)return (void)(putchar('-'),fwrit(-x));
if(x>9)fwrit(x/10);putchar(x%10^48);
}
inline LL mulMod(const LL a,const LL b,const LL mod){//long long multiplie_mod
return ((a*b-(LL)((long double)a/mod*b+1e-8)*mod)%mod+mod)%mod;
}
const int MAXN=1e5;
const int INF=0x3f3f3f3f;
class Splay{
private:
int rt;
int ff[MAXN+5];
int ch[MAXN+5][2];
int siz[MAXN+5];
int recy[MAXN+5];
int val[MAXN+5];
vector<int>pool;
inline void release(const int x){
if(x)pool.push_back(x);
}
inline int generate(){
int ret=pool.back();
pool.pop_back();
return ret;
}
public:
Splay(){
rt=0;
rep(i,1,MAXN)pool.push_back(i);
}
inline void pushup(const int x){
siz[x]=siz[ch[x][0]]+siz[ch[x][1]]+recy[x];
}
inline void Rotate(const int x){
int y=ff[x],z=ff[y],k=(ch[y][1]==x);
ch[z][ch[z][1]==y]=x;
ff[x]=z;
ch[y][k]=ch[x][k^1];
ff[ch[x][k^1]]=y;
ch[x][k^1]=y;
ff[y]=x;
pushup(y),pushup(x);
//注意 上传 的先后顺序
}
/*
void Rotate(int x)//旋转
{
register int y=ff[x];
register int z=ff[y];
register int k=ch[y][1]==x;//x是y的左或右儿子
ch[z][ch[z][1]==y]=x;
ff[x]=z;
ch[y][k]=ch[x][k^1];
ff[ch[x][k^1]]=y;
ch[x][k^1]=y;
ff[y]=x;
pushup(y);pushup(x);
}
*/
inline void splay(const int x,const int goal=0){
int y,z;
while(ff[x]!=goal){
y=ff[x],z=ff[y];
if(z!=goal)
(ch[z][0]==y)^(ch[y][0]==x)?Rotate(x):Rotate(y);
Rotate(x);
}
if(goal==0)rt=x;
pushup(x);
}
inline void Insert(const int x){
int u=rt,fa=0;
while(u && val[u]!=x)fa=u,u=ch[u][val[u]<x];
if(u)++recy[u];
else{
u=generate();
if(fa)ch[fa][val[fa]<x]=u;
ff[u]=fa,val[u]=x;
ch[u][0]=ch[u][1]=0;
siz[u]=recy[u]=1;
}
splay(u);
}
inline void Find(const int x){
int u=rt;
if(!u)return;
while(x!=val[u] && ch[u][val[u]<x])u=ch[u][val[u]<x];
splay(u);
}
inline int Next(const int x,const int f){
Find(x);
int u=rt;
if((f && val[u]>x) || (!f && val[u]<x))return u;
for(u=ch[u][f];ch[u][f^1];u=ch[u][f^1]);
return u;
}
inline void Delete(const int x){
int pre=Next(x,0);
int suf=Next(x,1);
// printf("pre_val == %d, suf_val == %d\n",val[pre],val[suf]);
splay(pre);
splay(suf,pre);
int del=ch[suf][0];
if(recy[del]>1){
--recy[del];
splay(del);
}
else{
if(ch[suf][0])release(ch[suf][0]);
ch[suf][0]=0;
splay(suf);
//这里需要 splay 吗 ?
}
}
inline int Atrnk(int rnk){
int u=rt,y;
if(siz[u]<rnk)return -INF;
while(233){
y=ch[u][0];
if(rnk>siz[y]+recy[u]){
rnk-=siz[y]+recy[u];
u=ch[u][1];
}
else if(rnk<=siz[y])u=ch[u][0];
else return val[u];
}
}
inline int Getrnk(const int x){
Find(x);
return siz[ch[rt][0]]+1;
}
inline int Getpre(const int x){
int u=Next(x,0);
return val[u];
}
inline int Getsuf(const int x){
int u=Next(x,1);
return val[u];
}
void chk_tre(int u=-1){
if(u==-1)u=rt;
printf("Now u == %d, ff == %d, val == %d, siz == %d, lc == %d, rc == %d, recy == %d\n",u,ff[u],val[u],siz[u],ch[u][0],ch[u][1],recy[u]);
if(ch[u][0])chk_tre(ch[u][0]);
if(ch[u][1])chk_tre(ch[u][1]);
}
}Tre;
int n,opt,x;
signed main(){
#ifdef FILEOI
freopen("file.in","r",stdin);
freopen("file.out","w",stdout);
#endif
Tre.Insert(INF);
Tre.Insert(-INF);
// Tre.chk_tre();
// Endl;
qread(n);
while(n--){
qread(opt,x);
if(opt==1)Tre.Insert(x);
else if(opt==2)Tre.Delete(x);
else if(opt==3)writc(Tre.Getrnk(x)-1,'\n');
else if(opt==4)writc(Tre.Atrnk(x+1),'\n');
else if(opt==5)writc(Tre.Getpre(x),'\n');
else writc(Tre.Getsuf(x),'\n');
// Tre.chk_tre();
// Endl;
}
return 0;
}
扩展版
这个板子是针对有 pushdown() 与 pushup() 操作的 Splay ，这里的模板题目是 POJ3580 。
#include<cstdio>
#include<vector>
using namespace std;
#define rep(i,__l,__r) for(signed i=__l,i##_end_=__r;i<=i##_end_;++i)
#define fep(i,__l,__r) for(signed i=__l,i##_end_=__r;i>=i##_end_;--i)
#define writc(a,b) fwrit(a),putchar(b)
#define mp(a,b) make_pair(a,b)
#define ft first
#define sd second
#define LL long long
#define ull unsigned long long
#define uint unsigned int
#define pii pair< int,int >
#define Endl putchar('\n')
// #define FILEOI
// #define int long long
// #define int unsigned
#ifdef FILEOI
# define MAXBUFFERSIZE 500000
inline char fgetc(){
static char buf[MAXBUFFERSIZE+5],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,MAXBUFFERSIZE,stdin),p1==p2)?EOF:*p1++;
}
# undef MAXBUFFERSIZE
# define cg (c=fgetc())
#else
# define cg (c=getchar())
#endif
template<class T>inline void qread(T& x){
char c;bool f=0;
while(cg<'0'||'9'<c)f|=(c=='-');
for(x=(c^48);'0'<=cg&&c<='9';x=(x<<1)+(x<<3)+(c^48));
if(f)x=-x;
}
inline int qread(){
int x=0;char c;bool f=0;
while(cg<'0'||'9'<c)f|=(c=='-');
for(x=(c^48);'0'<=cg&&c<='9';x=(x<<1)+(x<<3)+(c^48));
return f?-x:x;
}
template<class T,class... Args>inline void qread(T& x,Args&... args){qread(x),qread(args...);}
template<class T>inline T Max(const T x,const T y){return x>y?x:y;}
template<class T>inline T Min(const T x,const T y){return x<y?x:y;}
template<class T>inline T fab(const T x){return x>0?x:-x;}
inline int gcd(const int a,const int b){return b?gcd(b,a%b):a;}
inline void getInv(int inv[],const int lim,const int MOD){
inv[0]=inv[1]=1;for(int i=2;i<=lim;++i)inv[i]=1ll*inv[MOD%i]*(MOD-MOD/i)%MOD;
}
template<class T>void fwrit(const T x){
if(x<0)return (void)(putchar('-'),fwrit(-x));
if(x>9)fwrit(x/10);putchar(x%10^48);
}
inline LL mulMod(const LL a,const LL b,const LL mod){//long long multiplie_mod
return ((a*b-(LL)((long double)a/mod*b+1e-8)*mod)%mod+mod)%mod;
}
const int MAXN=1e5;
const int INF=0x3f3f3f3f;
class Splay{
private:
struct node{
int ff,ch[2],siz,val,lag,delta,minn;
node(){ff=ch[0]=ch[1]=siz=val=lag=delta=minn=0;}
};
// node t[MAXN+5];
#define ch(x,i) t[x].ch[i]
vector<node>t;
int rt,ncnt,__siz;
vector<int>pol;
#define SUBBUFFERSIZ 100005
inline void release(const int u){
if(u)pol.push_back(u);
}
inline void addBuffer(){
fep(i,__siz+SUBBUFFERSIZ,__siz+1)pol.push_back(i);
__siz+=SUBBUFFERSIZ;
}
inline int generate(){
if(pol.empty())addBuffer();
int ret=pol.back();
pol.pop_back();
if(ret>ncnt)++ncnt,t.push_back(node());
return ret;
}
public:
Splay(){
rt=ncnt=__siz=0;
addBuffer();
t.push_back(node());
t[0].minn=INF;//pay attention !
}
inline void pushup(const int x){
t[x].siz=t[ch(x,0)].siz+t[ch(x,1)].siz+1;
t[x].minn=Min(t[x].val,Min(t[ch(x,0)].minn,t[ch(x,1)].minn));
}
inline void add(const int u,const int delta){
if(!u)return;
t[u].val+=delta;
t[u].delta+=delta;
t[u].minn+=delta;
}
inline void revers(const int u){
if(!u)return;
swap(ch(u,0),ch(u,1));
t[u].lag^=1;
}
inline void pushdown(const int x){
if(t[x].delta){
add(ch(x,0),t[x].delta);
add(ch(x,1),t[x].delta);
t[x].delta=0;
}
if(t[x].lag){
revers(ch(x,0));
revers(ch(x,1));
t[x].lag=0;
}
}
inline void rotate(const int x){
const int y=t[x].ff,z=t[y].ff,k=(ch(y,1)==x);
pushdown(y),pushdown(x);
if(z)ch(z,ch(z,1)==y)=x;
t[x].ff=z;
ch(y,k)=ch(x,k^1);
t[ch(x,k^1)].ff=y;
ch(x,k^1)=y;
t[y].ff=x;
pushup(y),pushup(x);
}
inline void splay(const int x,const int goal=0){
int y,z;
pushdown(x);
//是否需要 pushdown ?
while(t[x].ff!=goal){
y=t[x].ff,z=t[y].ff;
if(z!=goal)
(ch(z,0)==y)^(ch(y,0)==x)?rotate(x):rotate(y);
rotate(x);
}
if(goal==0)rt=x;
pushup(x);
}
inline int rnkFind(int rnk,const int f=0){
int u=rt,y;
if(!u || t[u].siz<rnk)return 0;
while(233){
pushdown(u);
y=ch(u,0);
if(rnk>t[y].siz+1){
rnk-=t[y].siz+1;
u=ch(u,1);
}
else if(rnk<=t[y].siz)u=ch(u,0);
else return f?splay(u),u:u;
}
}
inline void rnkInsert(const int rnk,const int x){
rnkFind(rnk,1);
int u=rt,ins=generate(),ext=ch(rt,1);
// printf("Now new node ins == %d, u == %d\n",ins,u);
if(u)ch(u,1)=ins;
t[ins].ff=u;
ch(ins,0)=0,ch(ins,1)=ext;
if(ext)t[ext].ff=ins;
t[ins].siz=1,t[ins].val=t[ins].minn=x;
t[ins].lag=t[ins].delta=0;
// printf("the new node's info :\n");
// printf("ins == %d, lc == %d, rc == %d, ff == %d, val == %d, siz == %d, delta == %d, lag == %d\n",ins,t[ins].ch[0],t[ins].ch[1],t[ins].ff,t[ins].val,t[ins].siz,t[ins].delta,t[ins].lag);
splay(ins);
}
inline void rnkDelete(const int rnk){
int pre=rnkFind(rnk-1);
int suf=rnkFind(rnk+1);
splay(pre);
splay(suf,pre);
release(ch(suf,0));
ch(suf,0)=0;
splay(suf);
}
inline void Reverse(const int lrnk,const int rrnk){
int l=lrnk-1,r=rrnk+1;
int ln=rnkFind(l);
int rn=rnkFind(r);
splay(ln,0);
splay(rn,ln);
revers(ch(rn,0));
return;
}
inline void modify(const int lrnk,const int rrnk,const int delta){
int l=lrnk-1,r=rrnk+1;
int ln=rnkFind(l);
int rn=rnkFind(r);
splay(ln);
// printf("After the first splay, the tree :\n");
// chk_tre();
// Endl;
splay(rn,ln);
// printf("After the second splay, the tree :\n");
// chk_tre();
// Endl;
// printf("modify : ln == %d, rn == %d\n",ln,rn);
// printf("add node %d, delta == %d\n",ch(rn,0),delta);
add(ch(rn,0),delta);
splay(rn);
return;
}
inline int Getminn(const int lrnk,const int rrnk){
int l=lrnk-1,r=rrnk+1;
int ln=rnkFind(l);
int rn=rnkFind(r);
splay(ln);
splay(rn,ln);
return t[ch(rn,0)].minn;
}
inline void Revolve(const int lrnk,const int rrnk,const int k){
Reverse(rrnk-k+1,rrnk);
Reverse(lrnk,rrnk-k);
Reverse(lrnk,rrnk);
}
inline int Atrnk(const int rnk){
// rnkFind(rnk,1);
int u=rnkFind(rnk);
return t[u].val;
}
inline void chk_tre(const int u=-1,int rnk=0){
// if(rnk>10)return;
if(u==-1)return chk_tre(rt);
if(ch(u,0))chk_tre(ch(u,0),rnk);
rnk+=t[ch(u,0)].siz;
printf("rnk == %d\n",++rnk);
printf("Now u == %d, lc == %d, rc == %d, ff == %d, val == %d, siz == %d, delta == %d, lag == %d\n",u,t[u].ch[0],t[u].ch[1],t[u].ff,t[u].val,t[u].siz,t[u].delta,t[u].lag);
Endl;
if(ch(u,1))chk_tre(ch(u,1),rnk);
}
inline void Writtre(const int u=-1){
if(u==-1)return Writtre(rt);
pushdown(u);
if(ch(u,0))Writtre(ch(u,0));
if(t[u].val!=-INF && t[u].val!=INF)
printf("%d ",t[u].val);
if(ch(u,1))Writtre(ch(u,1));
}
}tre;
int n,m;
signed main(){
#ifdef FILEOI
freopen("file.in","r",stdin);
freopen("file.out","w",stdout);
#endif
qread(n);
tre.rnkInsert(0,-INF);
rep(i,1,n)tre.rnkInsert(i,qread());
tre.rnkInsert(n+1,INF);
// tre.chk_tre();
// Endl;
qread(m);
char opt[20];
int x,y;
while(m--){
scanf("%s",opt+1);
x=qread();
if(opt[1]=='D')tre.rnkDelete(x+1);
else y=qread();
if(opt[1]=='A')tre.modify(x+1,y+1,qread());
else if(opt[1]=='R' && opt[4]=='E')tre.Reverse(x+1,y+1);
else if(opt[1]=='R')tre.Revolve(x+1,y+1,qread()%(y-x+1));
else if(opt[1]=='I')tre.rnkInsert(x+1,y);
else if(opt[1]=='M')writc(tre.Getminn(x+1,y+1),'\n');
// puts("After this option, the tre : ");
// tre.chk_tre();
// tre.Writtre();
// Endl;
}
return 0;
}
「模板」Splay的更多相关文章
- 「模板」 FHQ_Treap
「模板」 FHQ_Treap 我也是偶然发现我还没发过FHQ_Treap的板子. 那就发一波吧. 这个速度实在不算快,但是不用旋转,并且好写. 更重要的是,Splay 可以做的事情它都可以做!比如区间 ...
- 「模板」 线段树——区间乘 && 区间加 && 区间求和
「模板」 线段树--区间乘 && 区间加 && 区间求和 原来的代码太恶心了,重贴一遍. #include <cstdio> int n,m; long l ...
- 「模板」 FHQ_Treap 区间翻转
「模板」 FHQ_Treap 区间翻转 没有旋转的 Treap 实现区间操作的功能,很好理解,也很好写,只是速度不算太快. 对于要翻转的区间,把整棵 Treap(存有区间 \([1,n]\) 的信息) ...
- 「模板」 树链剖分 HLD
「模板」 树链剖分 HLD 不懂OOP的OIer乱用OOP出人命了. 谨此纪念人生第一次类套类. 以及第一次OI相关代码打过200行. #include <algorithm> #incl ...
- 「模板」「讲解」Treap名次树
Treap实现名次树 前言 学平衡树的过程可以说是相当艰难.浏览Blog的过程中看到大量指针版平衡树,不擅长指针操作的我已经接近崩溃.于是,我想着一定要写一篇非指针实现的Treap的Blog. 具体如 ...
- 「模板」AC自动机
目录 说明 普通版本 询问更改版 拓扑优化版本 说明 这篇博客只挂模板,具体分析请膜拜大佬 hyfhaha 大佬. 普通版本 题目传送门 #include<cstdio> #include ...
- 「模板」可持久化 HFQ-Treap
老师用的是静态数组的写法,开了很多数组- 其实个人更倾向于 struct 或者用 class 封装起来. 但是鉴于太难打 好吧,是我懒得打. 然后就借鉴了老师的模板,写出了属于自己的 压行 风格. 代 ...
- Solution -「LOJ #138」「模板」类欧几里得算法
\(\mathcal{Description}\) Link. \(T\) 组询问,每次给出 \(n,a,b,c,k_1,k_2\),求 \[\sum_{x=0}^nx^{k_1}\left\ ...
- Solution -「LOJ #141」回文子串 ||「模板」双向 PAM
\(\mathcal{Description}\) Link. 给定字符串 \(s\),处理 \(q\) 次操作: 在 \(s\) 前添加字符串: 在 \(s\) 后添加字符串: 求 \(s\ ...
随机推荐
- 3个N加上各种运算符号结果等于6（纯属娱乐）C#
网上的题目: 题有点难 但都有解 2 2 2 = 6 3 3 3 = 6 4 4 4 = 6 5 5 5 = 6 6 6 ...
- C# interact with Command prompt
using System.IO; using System.Diagnostics; static void Main(string[] args) { CmdDemo("dir" ...
- Java后端API调用身份验证的思考
在如今信息泛滥的数字时代中对产品安全性的要求越来越高了,就比如说今天要讨论的Java后端API调用的安全性,在你提供服务的接口中一定要保证调用方身份的有效性和合法性,不能让非法的用户进行调用,避免数据 ...
- Intel 8086 常用汇编指令表
一.数据传输指令 它们在存贮器和寄存器.寄存器和输入输出端口之间传送数据. 1. 通用数据传送指令. MOV 传送字或字节. MOVSX 先符号扩展,再传送. MOVZX 先零扩展,再传送. PUSH ...
- The Way to Home CodeForces - 910A
4个月前做的一道题,当时不知道为什么,写了一个bfs,直接就超时了. 现在再看这个题目,发现就是一个简单的贪心,每次走最远即可. #include <bits/stdc++.h> usin ...
- 【你不知道的javaScript 中卷 笔记2】javaScript中的类型转换
1.1 对象内部属性 [[Class]] 常见的原生函数: String() Number() Boolean() Array() Object() Function() RegExp() Date( ...
- Linux一些基本命令、inode定义、软硬链接
1.创建普通文件命令:touch 命令 2.创建目录文件命令:mkdir 命令 3.删除普通文件命令:rm 命令 4.删除目录文件命令:rmdir 命令 5.给普通文件写东西命令:vim 命令 6.查 ...
- 回文串--Manacher算法（模板）
用途:在O(n)时间内,求出以每一个点为中心的回文串长度. 首先,有一个非常巧妙的转化.由于回文串长度有可能为奇数也有可能为偶数,说明回文中心不一定在一个字符上.所以要将字符串做如下处理:在每两个字母 ...
- Docker最全教程——从理论到实战（二十一）
前言 MySQL是目前最流行的开源的关系型数据库,MySQL的容器化之前有朋友投稿并且写过此块,本篇仅从笔者角度进行总结和编写. 目录 镜像说明 运行MySQL容器镜像 1.运行MySQL容器 ...
- springboot~集成DataSource 与 Druid监控配置
介绍 Druid首先是一个数据库连接池.Druid是目前最好的数据库连接池,在功能.性能.扩展性方面,都超过其他数据库连接池,Druid已经在阿里巴巴部署了超过600个应用,经过一年多生产环境大规模部 ...