「模板」可持久化 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\ ...
随机推荐
- maven私服搭建&使用
Maven私服搭建教程 一.nexus安装 1,解压安装包 安装包下载地址 2,以管理员身份打开cmd 3,进入到nexus的bin目录 (1) 安装:nexus install (2) 启动:nex ...
- C# 委托实例实现的多种类型
using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.T ...
- nginx配置访问密码,输入用户名和密码才能访问
1. 安装 htpasswd 工具 yum install httpd-tools -y 设置用户名和密码,并把用户名.密码保存到指定文件中: [sandu@bogon conf]$ sudo mkd ...
- Analog power pin UPF defination
在一个analog macro端口上,有些pin是always on的,有些是shut down的,如何描述这些pin的power属性?这是一个常见问题,驴就此机会大致描述一下常见的做法.对于这个问题 ...
- app内区域截图利用html2Canvals保存到手机 截屏 (html2Canvas使用版本是:0.5.0-beta3。)
app内区域截图利用html2Canvals保存到手机 app内有时候需要区域内的截图保存dom为图像,我们可以使用html2Canvas将dom转换成base64图像字符串,然后再利用5+api保存 ...
- IDF-CTF-简单的Elf逆向Writeup
ElfCrackMe1 *:first-child { margin-top: 0 !important; } body>*:last-child { margin-bottom: 0 !imp ...
- linux 配置php环境变量
vim /etc/profile //加上 export PATH=$PATH:/usr/local/php/bin 保存退出 source /etc/profile php -v 注:该配置对所有用 ...
- php 接口获取公网ip并获取天气接口信息
<?php function get_ip(){ //判断服务器是否允许$_SERVER if(isset($_SERVER)){ if(isset($_SERVER['HTTP_X_FORWA ...
- Spring中如何使用工厂模式实现程序解耦?
目录 1. 啥是耦合.解耦? 2. jdbc程序进行解耦 3.传统dao.service.controller的程序耦合性 4.使用工厂模式实现解耦 5.工厂模式改进 6.结语 @ 1. 啥是耦合.解 ...
- RS422接口与RS485接口
RS422具体接线参考网站 RS485接口 RS485设备为半双工设备,RS485收发器信号相关引脚包括控制引脚.485A.485B,其中控制引脚的高低电平决定当前处于接收模式还是发送模式. RS48 ...