【题解】BZOJ 3065: 带插入区间K小值——替罪羊树套线段树
题解
orz vfk的题解
3065: 带插入区间K小值 系列题解
惨
一开始用了一种空间常数很大的方法,每次重构的时候merge两颗线段树,然后无限RE(其实是MLE)。
后来改成枚举子树元素插入,空间缩小为约 \(\frac 1 4\) ,然而TLE。
然后把替罪羊树的 \(\alpha\) 从 0.6改成0.75,就卡过了。
代码
#include <bits/stdc++.h>
using namespace std;
const int MAXN=140005, MAXM=3e7, MAXB=2e7, MX=70000;
const double al=0.75;
char BUF[MAXB], *cp=BUF;
void rd(int &x){
x=0;
while(*cp<'0'||'9'<*cp)cp++;
while('0'<=*cp&&*cp<='9')x=x*10+*cp++-'0';
}
char rc(){while(*cp<'A'||'Z'<*cp)cp++; return *cp++;}
int N, M, L, R, TOP, top, ov;
int nt, tot, ok, A[MAXN];
struct Seg{
Seg *lc, *rc;
int s;
void *operator new(size_t);
void operator delete(void *);
Seg();
void up(){s=lc->s+rc->s;}
}tr[MAXM], *ST[MAXM], *tmp[MAXN];
void *Seg::operator new(size_t size){return ST[--TOP];}
void Seg::operator delete(void *p){ST[TOP++]=(Seg*)p;}
Seg::Seg(){lc=rc=tr;s=0;}
void dec(Seg *&x){
if(x==tr) return;
dec(x->lc); dec(x->rc); delete(x); x=tr;
}
void upd(Seg *&x, int l, int r, int k, int v){
if(x==tr) x=new Seg;
if(l==r){x->s+=v; return;}
int mid=(l+r)>>1;
if(k<=mid) upd(x->lc,l,mid,k,v);
else upd(x->rc,mid+1,r,k,v);
x->up();
if(!x->s) dec(x);
}
struct Node{
Node *lc, *rc;
Seg *c, *s;
int sz, v;
void up(){sz=1+lc->sz+rc->sz;}
}nd[MAXN], *st[MAXN], *root;
void dfs(Node *x){
if(x==nd) return; dec(x->s);
dfs(x->lc); st[top++]=x; dfs(x->rc);
}
void bu(Node *&x, int l, int r){
if(l>r){x=nd; return;}
int mid=(l+r)>>1; x=st[mid];
bu(x->lc,l,mid-1); bu(x->rc,mid+1,r);
x->up();
for(int i=l; i<=r; ++i) upd(x->s,0,MX,st[i]->v,1);
}
void rebu(Node *&x){top=0; dfs(x); bu(x,0,top-1);}
void ins(Node *&x, int k, int v, int d=0){
if(x==nd){
x=&nd[++tot]; x->v=v;
upd(x->c,0,MX,v,1);
upd(x->s,0,MX,v,1);
return;
}
upd(x->s,0,MX,v,1);
int t=x->lc->sz+1;
if(k<=t){
ins(x->lc,k,v,d+1); x->up();
if(x->lc->sz>=al*x->sz) rebu(x);
}else{
ins(x->rc,k-t,v,d+1); x->up();
if(x->rc->sz>=al*x->sz) rebu(x);
}
}
void md(Node *&x, int k, int v){
if(x==nd) return;
int t=x->lc->sz;
if(k==t+1){
ov=x->v; x->v=v;
dec(x->c); upd(x->c,0,MX,v,1);
}else if(k<=t) md(x->lc,k,v);
else md(x->rc,k-t-1,v);
upd(x->s,0,MX,ov,-1);
upd(x->s,0,MX,v,1);
}
void qry(Node *x, int l, int r){
if(x==nd) return;
if(L<=l&&r<=R){
tmp[nt++]=x->s;
return;
}
int t=x->lc->sz;
if(L<=l+t&&l+t<=R) tmp[nt++]=x->c;
if(L<l+t) qry(x->lc,l,l+t-1);
if(l+t<R) qry(x->rc,l+t+1,r);
}
int kth(int k){
int l=0, r=MX;
while(l<r){
int t=0, mid=(l+r)>>1;
for(int i=0; i<nt; ++i) t+=tmp[i]->lc->s;
if(k<=t){
r=mid;
for(int i=0; i<nt; ++i) tmp[i]=tmp[i]->lc;
}else{
l=mid+1; k-=t;
for(int i=0; i<nt; ++i) tmp[i]=tmp[i]->rc;
}
}
return l;
}
void init(){
tr[0].lc=tr[0].rc=tr;
for(int i=MAXM-1; i>0; --i) ST[TOP++]=tr+i;
root=nd[0].lc=nd[0].rc=nd;
nd[0].c=nd[0].s=tr;
for(int i=1; i<MAXN; ++i){
nd[i].c=nd[i].s=tr;
nd[i].lc=nd[i].rc=nd;
nd[i].sz=1;
}
for(int i=1; i<=N; ++i){
upd(nd[i].c,0,MX,A[i],1);
st[top++]=nd+i; nd[i].v=A[i];
}
bu(root,0,top-1); tot=N;
}
int main(){
fread(BUF, 1, MAXB, stdin);rd(N);
for(int i=1; i<=N; ++i) rd(A[i]);
init(); rd(M);int last=0;
while(M--){
char ch=rc();
int x,y,k; rd(x),rd(y);
x^=last; y^=last;
if(ch=='Q'){
L=x, R=y;
rd(k); k^=last;
nt=0; qry(root,1,N);
printf("%d\n", last=kth(k));
}else if(ch=='M')md(root,x,y);
else if(ch=='I')ins(root,x,y),N++;
}
return 0;
}
【题解】BZOJ 3065: 带插入区间K小值——替罪羊树套线段树的更多相关文章
- bzoj 3065: 带插入区间K小值 替罪羊树 && AC300
3065: 带插入区间K小值 Time Limit: 60 Sec Memory Limit: 512 MBSubmit: 1062 Solved: 253[Submit][Status] Des ...
- BZOJ 3065 带插入区间K小值 (替罪羊树套线段树)
毒瘤题.参考抄自博客:hzwer 第一次写替罪羊树,完全是照着题解写的,发现这玩意儿好强啊,不用旋转每次都重构还能nlognnlognnlogn. 还有外面二分和里面线段树的值域一样,那么r = mi ...
- BZOJ 3065 带插入区间K小值(sag套线段树)
3065: 带插入区间K小值 Time Limit: 60 Sec Memory Limit: 512 MBSubmit: 4696 Solved: 1527[Submit][Status][Di ...
- bzoj 3065: 带插入区间K小值(分块)
Description 从前有n只跳蚤排成一行做早操,每只跳蚤都有自己的一个弹跳力a[i].跳蚤国王看着这些跳蚤国欣欣向荣的情景,感到非常高兴.这时跳蚤国王决定理性愉悦一下,查询区间k小值.他每次向它 ...
- BZOJ 3065 带插入区间K小值
http://www.lydsy.com/JudgeOnline/problem.php?id=3065 思路:替罪羊树套权值线段树. 当替罪羊树某个子树大于某个比利(比例)时就暴力重构,本题时间复杂 ...
- 3065: 带插入区间K小值_树套树_替罪羊树_权值线段树
经过周六一天,周一3个小时的晚自习,周二2个小时的疯狂debug,终于凭借自己切掉了这道树套树题. Code: #include <cstdio> #include <algorit ...
- 【BZOJ3065】带插入区间K小值 替罪羊树+权值线段树
[BZOJ3065]带插入区间K小值 Description 从前有n只跳蚤排成一行做早操,每只跳蚤都有自己的一个弹跳力a[i].跳蚤国王看着这些跳蚤国欣欣向荣的情景,感到非常高兴.这时跳蚤国王决定理 ...
- 【BZOJ】3065: 带插入区间K小值
http://www.lydsy.com/JudgeOnline/problem.php?id=3065 题意:带插入.修改的区间k小值在线查询.(原序列n<=35000, 询问<=175 ...
- 【bzoj3065】带插入区间K小值 替罪羊树套权值线段树
题目描述 从前有n只跳蚤排成一行做早操,每只跳蚤都有自己的一个弹跳力a[i].跳蚤国王看着这些跳蚤国欣欣向荣的情景,感到非常高兴.这时跳蚤国王决定理性愉悦一下,查询区间k小值.他每次向它的随从伏特提出 ...
随机推荐
- AtCoder Grand Contest 019 A: Ice Tea Store
tourist出的题诶!想想就很高明,老年选手可能做不太动.不过A题还是按照惯例放水的. AtCoder Grand Contest 019 A: Ice Tea Store 题意:买0.25L,0. ...
- 【bzoj2829】信用卡凸包 凸包
题目描述 输入 输出 样例输入 26.0 2.0 0.00.0 0.0 0.02.0 -2.0 1.5707963268 样例输出 21.66 题解 凸包 傻逼题,答案显然为:所有圆心构成的凸包周长+ ...
- 【数据库_Postgresql】sql查询结果添加序号列
ROW_NUMBER () OVER (ORDER BY A .ordernumber ASC) AS 序号
- LCA的倍增算法
LCA,即树上两点之间的公共祖先,求这样一个公共祖先有很多种方法: 暴力向上:O(n) 每次将深度大的点往上移动,直至二者相遇 树剖:O(logn) 在O(2n)预处理重链之后,每次就将深度大的沿重链 ...
- 【原创】【目录】实现rich editor(富文本编辑器)教程,深入理解selectiona/range操作与浏览器差异
日常工作中,接触富文本编辑的次数还是很多,特发此教程,希望可以改变富文本编辑的技术门槛较高的现状,解决这部分的坑. 前提准备,兼容获取range,统一回车行为,前期准备工作 了解document.ex ...
- 【loj2461】【2018集训队互测Day 1】完美的队列
#2461. 「2018 集训队互测 Day 1」完美的队列 传送门: https://loj.ac/problem/2461 题解: 直接做可能一次操作加入队列同时会弹出很多数字,无法维护:一个操作 ...
- CSS中background样式的repeat和no-repeat
自http://blog.sina.com.cn/s/blog_98eef7830101cv8t.html获取的知识,特此感谢. <!DOCTYPE html> <html lang ...
- selenium - switch_to_alert() - 警告框处理
在WebDriver中处理JavaScript所生成的alert.confirm以及prompt十分简单,具体做法是使用 switch_to.alert 方法定位到 alert/confirm/pro ...
- C/C++ Volatile关键词深度剖析
文章来源:http://hedengcheng.com/?p=725 背景 此微博,引发了朋友们的大量讨论:赞同者有之:批评者有之:当然,更多的朋友,是希望我能更详细的解读C/C++ Volatile ...
- centos7通过yum安装MySQL
一:去官网查看最新安装包 https://dev.mysql.com/downloads/repo/yum/ 二:下载MySQL源安装包 wget http://dev.mysql.com/get/m ...