非旋treap模板
bzoj3580
非旋转treap
在大神教导下发现split一段区间时先split右边再split左边比较好写
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cctype>
using namespace std;
typedef long long LL;
const int M=200007;
inline int rd(){
int x=0;bool f=1;char c=getchar();
for(;!isdigit(c);c=getchar())if(c=='-')f=0;
for(;isdigit(c);c=getchar())x=x*10-48+c;
return f?x:-x;
}
struct D{
int c[2];
D(){c[0]=c[1]=0;}
};
struct treap{
int sz,val,fix,mn,rev,tag,c[2];
void input(){
val=mn=rd();
fix=rand();
sz=1;
rev=tag=c[0]=c[1]=0;
}
}a[M];
int n,m,rt;
inline int min(int x,int y){
return x<y?x:y;
}
void update(int x){
a[x].mn=a[x].val,a[x].sz=1;
if(a[x].c[0]) a[x].mn=min(a[x].mn,a[a[x].c[0]].mn),a[x].sz+=a[a[x].c[0]].sz;
if(a[x].c[1]) a[x].mn=min(a[x].mn,a[a[x].c[1]].mn),a[x].sz+=a[a[x].c[1]].sz;
}
void totag(int x,int z){
a[x].tag+=z;
a[x].mn+=z;
a[x].val+=z;
}
void pushdown(int x){
if(a[x].rev){
if(a[x].c[0])a[a[x].c[0]].rev^=1;
if(a[x].c[1])a[a[x].c[1]].rev^=1;
a[x].rev^=1;
swap(a[x].c[0],a[x].c[1]);
}
if(a[x].tag){
if(a[x].c[0]) totag(a[x].c[0],a[x].tag);
if(a[x].c[1]) totag(a[x].c[1],a[x].tag);
a[x].tag=0;
}
}
int build(int l,int r){
if(l>r) return 0;
int mid=l+r>>1;
a[mid].c[0]=build(l,mid-1);
a[mid].c[1]=build(mid+1,r);
update(mid);
return mid;
}
D split(int x,int k){
D nw;
if(x==0)return nw;
pushdown(x);
if(a[a[x].c[0]].sz+1==k){
nw.c[1]=a[x].c[1];
a[x].c[1]=0;
nw.c[0]=x;
}
else if(k<=a[a[x].c[0]].sz){
nw=split(a[x].c[0],k);
a[x].c[0]=nw.c[1];
nw.c[1]=x;
}
else{
nw=split(a[x].c[1],k-a[a[x].c[0]].sz-1);
a[x].c[1]=nw.c[0];
nw.c[0]=x;
}
update(x);
return nw;
}
int merge(int x,int y){
if(!x) return y;
if(!y) return x;
if(a[x].fix<a[y].fix){
pushdown(x);
a[x].c[1]=merge(a[x].c[1],y);
update(x);
return x;
}
else{
pushdown(y);
a[y].c[0]=merge(x,a[y].c[0]);
update(y);
return y;
}
}
void add(int x,int y,int z){
D t2=split(rt,y);
D t1=split(t2.c[0],x-1);
totag(t1.c[1],z);
rt=merge(merge(t1.c[0],t1.c[1]),t2.c[1]);
}
void ins(int x,int z){
D tt=split(rt,x);
rt=merge(merge(tt.c[0],z),tt.c[1]);
}
void del(int x){
D t2=split(rt,x);
D t1=split(t2.c[0],x-1);
rt=merge(t1.c[0],t2.c[1]);
}
void get(int x,int y){
D t2=split(rt,y);
D t1=split(t2.c[0],x-1);
printf("%d\n",a[t1.c[1]].mn);
rt=merge(merge(t1.c[0],t1.c[1]),t2.c[1]);
}
void reverse(int x,int y){
D t2=split(rt,y);
D t1=split(t2.c[0],x-1);
a[t1.c[1]].rev^=1;
rt=merge(merge(t1.c[0],t1.c[1]),t2.c[1]);
}
void revolve(int x,int y,int z){
D t2=split(rt,y);
D t1=split(t2.c[0],x-1);
D tt=split(t1.c[1],z);
rt=merge(merge(merge(t1.c[0],tt.c[1]),tt.c[0]),t2.c[1]);
}
int main(){
int i,x,y,z;
char s[9];
n=rd();
for(i=1;i<=n;i++) a[i].input();
rt=build(1,n);
m=rd();
for(i=1;i<=m;i++){
scanf("%s",s);
if(s[0]=='A'){
x=rd(),y=rd(),z=rd();
add(x,y,z);
}
else if(s[0]=='I'){
x=rd();
a[++n].input();
ins(x,n);
}
else if(s[0]=='D'){
x=rd();
del(x);
}
else if(s[0]=='M'){
x=rd(),y=rd();
get(x,y);
}
else if(s[3]=='E'){
x=rd(),y=rd();
reverse(x,y);
}
else{
x=rd(),y=rd(),z=rd();
z=(z%(y-x+1)+(y-x+1))%(y-x+1);
z=(y-x+1)-z;
revolve(x,y,z);
}
}
return 0;
}
非旋treap模板的更多相关文章
- 平衡树简单教程及模板(splay, 替罪羊树, 非旋treap)
原文链接https://www.cnblogs.com/zhouzhendong/p/Balanced-Binary-Tree.html 注意是简单教程,不是入门教程. splay 1. 旋转: 假设 ...
- [模板] 平衡树: Splay, 非旋Treap, 替罪羊树
简介 二叉搜索树, 可以维护一个集合/序列, 同时维护节点的 \(size\), 因此可以支持 insert(v), delete(v), kth(p,k), rank(v)等操作. 另外, prev ...
- 非旋 treap 结构体数组版(无指针)详解,有图有真相
非旋 $treap$ (FHQ treap)的简单入门 前置技能 建议在掌握普通 treap 以及 左偏堆(也就是可并堆)食用本blog 原理 以随机数维护平衡,使树高期望为logn级别, FHQ ...
- 2018.08.27 rollcall(非旋treap)
描述 初始有一个空集,依次插入N个数Ai.有M次询问Bj,表示询问第Bj个数加入集合后的排名为j的数是多少 输入 第一行是两个整数N,M 接下来一行有N个整数,Ai 接下来一行有M个整数Bj,保证数据 ...
- 非旋Treap总结 : 快过Splay 好用过传统Treap
非旋$Treap$ 其高级名字叫$Fhq\ Treap$,既然叫$Treap$,它一定满足了$Treap$的性质(虽然可能来看这篇的人一定知道$Treap$,但我还是多说几句:$Fhp\ Treap$ ...
- 2827: 千山鸟飞绝 非旋treap
国际惯例的题面:看起来很不可做的样子,我们先来整理一下题意吧.就是,维护每个点曾经拥有过的最大的两个属性值,支持把点的位置移动.我们用map对每个位置进行离散化,对每个位置建立一个平衡树.为了方便分离 ...
- 2081.09.22 Kuma(非旋treap)
描述 有N张卡片,编号从0到n-1, 刚开始从0到n-1按顺序排好. 现有一个操作, 对于p. l,表示从第p张卡片之后的l张卡片拿到 最前面. 例如n=7的时候, 刚开始卡片序列为0 1 2 3 4 ...
- 2018.08.06 bzoj1500: [NOI2005]维修数列(非旋treap)
传送门 平衡树好题. 我仍然是用的fhqtreap,感觉速度还行. 维护也比线段树splay什么的写起来简单. %%%非旋treap大法好. 代码: #include<bits/stdc++.h ...
- 2018.08.05 bzoj3223: Tyvj 1729 文艺平衡树(非旋treap)
传送门 经典的平衡树问题,之前已经用splay写过一次了,今天我突发奇想,写了一发非旋treap的版本,发现挺好写的(虽然跑不过splay). 代码: #include<bits/stdc++. ...
随机推荐
- time step和采样频率的关系
当前的采样频率为11025HZ,overlap=0.5,取的是1024个采样点为1帧.则time step为256时的时间长度为11秒.即以11秒为单位分割原始的音频,生成一张语普图. >> ...
- 【BZOJ-4726】Sabota? 树形DP
4726: [POI2017]Sabota? Time Limit: 20 Sec Memory Limit: 128 MBSec Special JudgeSubmit: 128 Solved ...
- 架构师养成记--12.Concurrent工具类CyclicBarrier和CountDownLatch
java.util.concurrent.CyclicBarrier 一组线程共同等待,直到达到一个公共屏障点. 举个栗子,百米赛跑中,所有运动员都要等其他运动员都准备好后才能一起跑(假如没有发令员) ...
- Java防止SQL注入2(通过filter过滤器功能进行拦截)
首先说明一点,这个过滤器拦截其实是不靠谱的,比如说我的一篇文章是介绍sql注入的,或者评论的内容是有关sql的,那会过滤掉:且如果每个页面都经过这个过滤器,那么效率也是非常低的. 如果是要SQL注入拦 ...
- Python为8bit深度图像应用color map
图片中存在着色版的概念,二维矩阵的每个元素的值指定了一种颜色,因此可以显示出彩色. 迁移调色板 下述python代码将VOC数据集中的某个语义分割的图片的调色板直接应用在一个二维矩阵代表的图像上 #l ...
- Linux磁盘分区及配额
在现有磁盘的基础上进行分区格式化并为特定用户实施磁盘配额,使其对磁盘这一分区的写入有一定的限制 前期准备: 在我的虚拟机rhel7上有/dev/sda这一分区和fsy这一用户,我将对/dev/sda进 ...
- JAVA中保留小数的多种方法
// 方式一:double f = 3.1516;BigDecimal b = new BigDecimal(f);double f1 = b.setScale(2, BigDecimal.ROUND ...
- Average Precision of VOC
转一篇文章,主要是关于VOC中Average Precision指标的 原文出处:https://sanchom.wordpress.com/tag/average-precision/ 还有一篇文章 ...
- web前端基础知识-(六)web框架
一.web框架本质 众所周知,对于所有的Web应用,本质上其实就是一个socket服务端,用户的浏览器其实就是一个socket客户端. #!/usr/bin/env python #coding:ut ...
- Win10 磁盘占用 100% 有效解决办法
立即查看 任务管理器,看看是不是有一个 服务主机 unistack服务组或者找类似名称的,点开以后你会看到里面有同步主机 blablah请你毫不犹豫的结束它!结束它!结束它! 按下WIN+R调出运行, ...