「模板」可持久化 HFQ-Treap
老师用的是静态数组的写法,开了很多数组…
其实个人更倾向于 struct 或者用 class 封装起来。
但是鉴于太难打 好吧,是我懒得打。
然后就借鉴了老师的模板,写出了属于自己的 压行 风格。
代码见下:(注:题目是洛谷P3835)
#include<cstdio>
#include<cstdlib>
#define rep(i,__l,__r) for(int i=__l,i##_end_=__r;i<=i##_end_;++i)
#define fep(i,__l,__r) for(int 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 pii pair<int,int>
#define Endl putchar('\n')
// #define FILEOI
#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=5e7;
const int INF=(1ll<<31)-1;
int ncnt;
int son[MAXN+5][2];
int keyof[MAXN+5];
int rd[MAXN+5];
int siz[MAXN+5];
int rt[MAXN+5];
int rcnt;
inline void pushup(const int p){
siz[p]=siz[son[p][0]]+siz[son[p][1]]+1;
}
inline int newnode(const int key){
++ncnt;
son[ncnt][0]=son[ncnt][1]=0;
keyof[ncnt]=key;
rd[ncnt]=rand();
siz[ncnt]=1;
return ncnt;
}
struct Pair{
int ele[2];
Pair(){ele[0]=ele[1]=0;}
Pair(const int A,const int B){ele[0]=A,ele[1]=B;}
};
inline int cpyNode(const int p){
++ncnt;
son[ncnt][0]=son[p][0];
son[ncnt][1]=son[p][1];
keyof[ncnt]=keyof[p];
rd[ncnt]=rd[p];
siz[ncnt]=siz[p];
return ncnt;
}
Pair split(int p,const int key,const bool mode,const bool flg=0){
//mode=0 时为按 key 分裂
//mode=1 时为按 rank 分裂
Pair ret=Pair(0,0);
if(p){
int d;
if(!flg)p=cpyNode(p);
if(mode)d=(siz[son[p][0]]+1<=key);
else d=(keyof[p]<=key);
if(mode && d)ret=split(son[p][d],key-siz[son[p][0]]-1,mode,flg);
else ret=split(son[p][d],key,mode,flg);
son[p][d]=ret.ele[d^1];
pushup(p);
ret.ele[d^1]=p;
}
return ret;
}
int unite(int r1,int r2,const bool flg=0){
if(r1*r2==0)return r1|r2;
if(rd[r1]<rd[r2]){
if(!flg)r1=cpyNode(r1);
son[r1][1]=unite(son[r1][1],r2,flg);
return pushup(r1),r1;
}
else{
if(!flg)r2=cpyNode(r2);
son[r2][0]=unite(r1,son[r2][0],flg);
return pushup(r2),r2;
}
}
int inser(const int p,const int key){
Pair tmp=split(p,key,0);
return unite(unite(tmp.ele[0],newnode(key)),tmp.ele[1]);
}
int delet(const int p,const int key){
Pair tmp0=split(p,key,0);
Pair tmp1=split(tmp0.ele[0],key-1,0);
return unite(tmp1.ele[0],tmp0.ele[1]);
}
int getRank(const int p,const int key){
if(!p)return 1;
if(keyof[p]>key)return getRank(son[p][0],key);
if(keyof[p]<key)return siz[son[p][0]]+1+getRank(son[p][1],key);
return siz[son[p][0]]+1;
}
int atRank(const int p,const int rk){
if(!p)return -INF;
if(siz[son[p][0]]>=rk)return atRank(son[p][0],rk);
else if(siz[son[p][0]]+1<rk)return atRank(son[p][1],rk-siz[son[p][0]]-1);
else return keyof[p];
}
int getPre(const int p,const int key){
Pair tmp1=split(p,key-1,0,1);
int tsiz=siz[tmp1.ele[0]]-1;
Pair tmp0=split(tmp1.ele[0],tsiz,1,1);
int ret=keyof[tmp0.ele[1]];
tmp1.ele[0]=unite(tmp0.ele[0],tmp0.ele[1],1);
unite(tmp1.ele[0],tmp1.ele[1],1);
return ret;
}
int getSuf(const int p,const int key){
Pair tmp0=split(p,key,0,1);
Pair tmp1=split(tmp0.ele[1],1,1,1);
int ret=keyof[tmp1.ele[0]];
tmp0.ele[1]=unite(tmp1.ele[0],tmp1.ele[1],1);
unite(tmp0.ele[0],tmp0.ele[1],1);
return ret;
}
int N,v,opt,x;
signed main(){
#ifdef FILEOI
freopen("file.in","r",stdin);
freopen("file.out","w",stdout);
#endif
srand(233);
qread(N);
rt[0]=inser(rt[0],-INF);
rt[0]=inser(rt[0],INF);
while(N--){
qread(v,opt,x);
rt[++rcnt]=rt[v];
if(opt==1)rt[rcnt]=inser(rt[rcnt],x);
else if(opt==2)rt[rcnt]=delet(rt[rcnt],x);
else if(opt==3)writc(getRank(rt[rcnt],x)-1,'\n');
else if(opt==4)writc(atRank(rt[rcnt],x+1),'\n');
else if(opt==5)writc(getPre(rt[rcnt],x),'\n');
else writc(getSuf(rt[rcnt],x),'\n');
}
return 0;
}
「模板」可持久化 HFQ-Treap的更多相关文章
- 「洛谷 P3834」「模板」可持久化线段树 题解报告
题目描述 给定n个整数构成的序列,将对于指定的闭区间查询其区间内的第k小值. 输入输出格式 输入格式 第一行包含两个正整数n,m,分别表示序列的长度和查询的个数. 第二行包含n个整数,表示这个序列各项 ...
- 「模板」 FHQ_Treap
「模板」 FHQ_Treap 我也是偶然发现我还没发过FHQ_Treap的板子. 那就发一波吧. 这个速度实在不算快,但是不用旋转,并且好写. 更重要的是,Splay 可以做的事情它都可以做!比如区间 ...
- 「模板」 FHQ_Treap 区间翻转
「模板」 FHQ_Treap 区间翻转 没有旋转的 Treap 实现区间操作的功能,很好理解,也很好写,只是速度不算太快. 对于要翻转的区间,把整棵 Treap(存有区间 \([1,n]\) 的信息) ...
- 「模板」 线段树——区间乘 && 区间加 && 区间求和
「模板」 线段树--区间乘 && 区间加 && 区间求和 原来的代码太恶心了,重贴一遍. #include <cstdio> int n,m; long l ...
- 「模板」 树链剖分 HLD
「模板」 树链剖分 HLD 不懂OOP的OIer乱用OOP出人命了. 谨此纪念人生第一次类套类. 以及第一次OI相关代码打过200行. #include <algorithm> #incl ...
- 「模板」「讲解」Treap名次树
Treap实现名次树 前言 学平衡树的过程可以说是相当艰难.浏览Blog的过程中看到大量指针版平衡树,不擅长指针操作的我已经接近崩溃.于是,我想着一定要写一篇非指针实现的Treap的Blog. 具体如 ...
- 「模板」Splay
代码说明 对于一些变量进行说明: 变量名 说明 rt 树根 ff[u] 点 \(u\) 的父节点,特别地, ff[rt]=0 ch[u][0|1] 点 \(u\) 的 左/右儿子 siz[u] 点 \ ...
- 「模板」AC自动机
目录 说明 普通版本 询问更改版 拓扑优化版本 说明 这篇博客只挂模板,具体分析请膜拜大佬 hyfhaha 大佬. 普通版本 题目传送门 #include<cstdio> #include ...
- Solution -「LOJ #138」「模板」类欧几里得算法
\(\mathcal{Description}\) Link. \(T\) 组询问,每次给出 \(n,a,b,c,k_1,k_2\),求 \[\sum_{x=0}^nx^{k_1}\left\ ...
随机推荐
- date命令的帮助信息,使用date命令输出数字做为命名标题则不会有重复标题
date命令的帮助信息,如下图 原文来自 https://blog.csdn.net/yz18931904/article/details/80985345 [root@localhost sourc ...
- tmux的基本使用
tmux的基本使用 tmux 最近发现了一个linux终端非常好用的工具,可以快速分屏 使用方法 参考
- vue Cli 按需引入Element UI 和全局引用Element UI
全局引用: 一.安装 Element UI npm i element-ui -S 二.在main.js 中引入 element UI import ElementUI from 'element-u ...
- 线段树+扫描线 NAIPC 2019 Intersecting Rectangles
你看看你有多菜,一点线段树的小小的运用,就不会写了: 题意:如果矩阵有交集,输出1,否则输出0(不包含内嵌): 思路:本题求交集,还得不包括内嵌的情况: 做过一道是求面积的题.跟这道类似,但在这里定义 ...
- 【游戏体验】Colour My World(让我的世界充满色彩)
这是一款玩法简单的游戏 但是它给你的感觉不一样 推荐试玩 个人测评 游戏性 7/10 音乐 9/10 剧情 6/10 总评 22/30
- Hello 2020D(多重集)
如果有一对时间对在某一场馆有时间重合而这一对时间对在另一场馆没有时间重合,则输出NO,否则输出YES. #define HAVE_STRUCT_TIMESPEC #include<bits/st ...
- JavaScript 推箱子游戏
推箱子游戏的 逻辑非常简单,但是如果不动手的话,还是不太清楚.我在这里讲一下自己的思路. 制作推箱子,首先要有自己的设计素材.如下我也是网上找的素材 第二步,理清游戏的规则. 游戏规则: 1.小人将箱 ...
- Django报错 The serializer field might be named incorrectly and not match any Got AttributeError when attempting to get a value for field `author_for` on serializer `KnownledgeBaseListSerializer`
1.问题描述,在设置,model部分字段的serialier时,出现如下报错 字段如下: # 知识库List class KnownledgeBaseListSerializer(serializer ...
- Nginx开始(自己使用,简单开始)
Nginx 1.网址 2.安装 如果需要资料的话可以直接来下载进行下载 前提 sudo yum install yum-utils 云主机重启之后需要执行 nginx -c /etc/nginx/ng ...
- MySQL双机热备环境搭建
一. 前期准备 准备两台服务器(电脑),接入到同一局域网中,能够使双方可以ping通: 安装MySQL数据库,具体安装方法网上很全面,但是安装的版本需保持一致: 服务器IP地址设置. l A服 ...