以前看到这题都瑟瑟发抖,终于过了心情舒畅。

按下标为关键字建替罪羊树,每个结点开一个权值线段树,维护的这个结点代表的子树的信息。

这题还得垃圾回收,自己yy的,不知对不对..

#include <cstdio>
#include <algorithm>
using namespace std;
#define M ((L+R)>>1)
#define l(x) t[x].s[0]
#define r(x) t[x].s[1]
#define L(x) t2[x].s[0]
#define R(x) t2[x].s[1] const int N=;
char op[];
int n,x,y,z,q,rt,tt,t1,la,tp,tp2,V,st[N],st2[N],rb[N*],rb2[N*];
struct nd {int x,v,sz,s[];}t[N*],t2[N*]; void dl(int x) {if(!x) return; rb2[++tt]=x,dl(L(x)),dl(R(x)),t2[x].x=L(x)=R(x)=;}
void upd(int &o,int L,int R,int x,int y) {
if(!o) o=rb2[tt--]; t2[o].x+=y;
if(L==R) return;
if(x<=M) upd(L(o),L,M,x,y); else upd(R(o),M+,R,x,y);
if(!t2[o].x) dl(o),o=;
}
void ge(int x) {if(!x) return; ge(l(x)),st[++tp]=x,ge(r(x));}
void mr(int l,int r,int &x) {
if(!l&&!r) return; t2[x=rb2[tt--]].x=t2[l].x+t2[r].x,mr(L(l),L(r),L(x)),mr(R(l),R(r),R(x));
}
void bd(int &x,int L,int R) {
if(L>R) return; t[x=rb[t1--]].v=t[st[M]].v;
if(L==R) t[x].sz=;
else bd(l(x),L,M-),bd(r(x),M+,R),mr(t[l(x)].x,t[r(x)].x,t[x].x),t[x].sz=t[l(x)].sz+t[r(x)].sz+;
upd(t[x].x,,N,t[x].v,);
}
void rm() {for(int i=;i<=tp;i++) dl(t[st[i]].x),l(st[i])=r(st[i])=t[st[i]].x=t[st[i]].v=t[st[i]].sz=;}
void ins(int &x,int p,int v) {
if(!x) x=rb[t1--]; upd(t[x].x,,N,v,),t[x].sz++;
if(!t[x].v) {t[x].v=v; return;}
if(t[l(x)].sz+>=p) ins(l(x),p,v); else ins(r(x),p-t[l(x)].sz-,v);
if(t[x].sz*0.75<max(t[l(x)].sz,t[r(x)].sz)) tp=,ge(x),bd(x,,tp),rm();
}
void gt(int x,int l,int r) {
if(r-l+==t[x].sz) {st[++tp]=t[x].x; return;}
if(l<=t[l(x)].sz) gt(l(x),l,min(r,t[l(x)].sz));
if(l<=t[l(x)].sz+&&r>=t[l(x)].sz+) st2[++tp2]=t[x].v;
if(r>t[l(x)].sz+) gt(r(x),max(l-t[l(x)].sz-,),r-t[l(x)].sz-);
}
int dfs(int L,int R,int k) {
if(L==R) return L; int s=;
for(int i=;i<=tp;i++) s+=t2[L(st[i])].x;
for(int i=;i<=tp2;i++) if(st2[i]<=M&&st2[i]>=L) s++;
if(s>=k) {
for(int i=;i<=tp;i++) st[i]=L(st[i]);
return dfs(L,M,k);
} else {
for(int i=;i<=tp;i++) st[i]=R(st[i]);
return dfs(M+,R,k-s);
}
}
int qr(int l,int r,int k) {tp=tp2=,gt(rt,l,r); return dfs(,N,k);}
void gai(int x,int p,int v) {
if(t[l(x)].sz+==p) V=t[x].v,t[x].v=v; else if(t[l(x)].sz>=p) gai(l(x),p,v); else gai(r(x),p-t[l(x)].sz-,v);
upd(t[x].x,,N,v,),upd(t[x].x,,N,V,-);
} int main() {
for(int i=;i<N*;i++) rb[i]=i;
for(int i=;i<N*;i++) rb2[i]=i;
scanf("%d",&n),tt=N*-,t1=N*-;
for(int i=;i<=n;i++) scanf("%d",&x),ins(rt,i,x+);
scanf("%d",&q);
while(q--) {
scanf("%s%d%d",op,&x,&y),x^=la,y^=la;
if(op[]=='Q') scanf("%d",&z),printf("%d\n",la=qr(x,y,z^la)-);
else if(op[]=='I') ins(rt,x,y+);
else gai(rt,x,y+);
}
return ;
}

BZOJ3065(替罪羊树套线段树)的更多相关文章

  1. [BZOJ 1901] Dynamic Rankings 【树状数组套线段树 || 线段树套线段树】

    题目链接:BZOJ - 1901 题目分析 树状数组套线段树或线段树套线段树都可以解决这道题. 第一层是区间,第二层是权值. 空间复杂度和时间复杂度均为 O(n log^2 n). 线段树比树状数组麻 ...

  2. ZJOI 2017 树状数组(线段树套线段树)

    题意 http://uoj.ac/problem/291 思路 不难发现,九条カレン醬所写的树状数组,在查询区间 \([1,r]\) 的时候,其实在查询后缀 \([r,n]\) :在查询 \([l,r ...

  3. BZOJ4317Atm的树&BZOJ2051A Problem For Fun&BZOJ2117[2010国家集训队]Crash的旅游计划——二分答案+动态点分治(点分树套线段树/点分树+vector)

    题目描述 Atm有一段时间在虐qtree的题目,于是,他满脑子都是tree,tree,tree…… 于是,一天晚上他梦到自己被关在了一个有根树中,每条路径都有边权,一个神秘的声音告诉他,每个点到其他的 ...

  4. hdu-4819-线段树套线段树

    http://acm.hdu.edu.cn/showproblem.php?pid=4819 给出一个N*N的矩阵,每次询问一个m*m的子矩阵里的floor((maxv+minv)/2)并把中间的元素 ...

  5. dfs序+主席树 或者 树链剖分+主席树(没写) 或者 线段树套线段树 或者 线段树套splay 或者 线段树套树状数组 bzoj 4448

    4448: [Scoi2015]情报传递 Time Limit: 20 Sec  Memory Limit: 256 MBSubmit: 588  Solved: 308[Submit][Status ...

  6. 【vijos】1750 建房子(线段树套线段树+前缀和)

    https://vijos.org/p/1750 是不是我想复杂了.... 自己yy了个二维线段树,然后愉快的敲打. 但是wa了两法.......sad 原因是在处理第二维的更新出现了个小问题,sad ...

  7. [TJOI2017][bzoj4889] 不勤劳的图书管理员 [线段树套线段树]

    题面 传送门 思路 考虑两本书的位置交换对答案的贡献: (为了方便描述,用"左边那本"和"右边那本"称呼两本我们要交换的书,"中间那本"是我 ...

  8. 【bzoj4785】[Zjoi2017]树状数组 线段树套线段树

    题目描述 漆黑的晚上,九条可怜躺在床上辗转反侧.难以入眠的她想起了若干年前她的一次悲惨的OI 比赛经历.那是一道基础的树状数组题.给出一个长度为 n 的数组 A,初始值都为 0,接下来进行 m 次操作 ...

  9. BZOJ1901 ZOJ2112 线段树+treap (线段树套线段树)

    BZOJ1901: 线段树套线段树做法: (外层线段树 里层动态开节点的权值线段树) 有一个小小的trick 可以省掉二分变成nlog^2n的 就是把查询的区间都取出来- logn个一起走- 2016 ...

  10. [BZOJ 3110] [luogu 3332] [ZJOI 2013]k大数查询(权值线段树套线段树)

    [BZOJ 3110] [luogu 3332] [ZJOI 2013]k大数查询(权值线段树套线段树) 题面 原题面有点歧义,不过从样例可以看出来真正的意思 有n个位置,每个位置可以看做一个集合. ...

随机推荐

  1. ThreadLocal就是这么简单

    前言 今天要研究的是ThreadLocal,这个我在一年前学习JavaWeb基础的时候接触过一次,当时在baidu搜出来的第一篇博文ThreadLocal,在评论下很多开发者认为那博主理解错误,给出了 ...

  2. MMA8451重力加速度计通过写内部校准寄存器进行校准

    |版权声明:本文为博主原创文章,未经博主允许不得转载. AN4069应用笔记中提到MMA8451的三个轴重力校准有两种方法, 第一种方法是简易校准,将贴有MMA8451的设备整体,Z轴正面朝上放在校准 ...

  3. PHP环境手动搭建wamp-----Apache+MySQL+PHP

    首先下载分别下载Apache+MySQL+PHP. 然后分别解压到文件夹中. 1.安装Apache 1)检查80端口是否占用 说明:apache软件占用80软件,在计算机中一个端口只能被一个软件占用 ...

  4. ASP.NET Web API编程——路由

    路由过程大致分为三个阶段: 1)请求URI匹配已存在路由模板 2)选择控制器 3)选择操作 1匹配已存在的路由模板 路由模板 在WebApiConfig.Register方法中定义路由,例如模板默认生 ...

  5. Web Api 使用模型验证

    public class Person { public int Id { get; set; } [Required(ErrorMessage = "姓名不能为空啊啊啊!")] ...

  6. jquery 实时监听输入框值变化方法

    $('.offers-number').bind('input propertychange', function (a, b) { var value = $(this).val() if (!va ...

  7. eclipse版本对应的jdk版本

    Installing Eclipse is relatively easy, but does involve a few steps and software from at least two d ...

  8. NHibernate与IbatisNet的简单比较

    NHibernate是当前最流行的Java O/R mapping框架Hibernate的移植版本,当前版本是1.0 rc-1.它出身于sf.net..IbatisNet是另外一种优秀的Java O/ ...

  9. C#制作ActiveX插件

    首先新建项目--->类库,取名:ActiveXDemo 右键项目属性:应用属性==>程序集信息=>使程序集Com可见, 生成==>输出==>为com互操作注册 新建接口类 ...

  10. python、java实现二叉树,细说二叉树添加节点、深度优先(先序、中序、后续)遍历 、广度优先 遍历算法

    数据结构可以说是编程的内功心法,掌握好数据结构真的非常重要.目前基本上流行的数据结构都是c和c++版本的,我最近在学习python,尝试着用python实现了二叉树的基本操作.写下一篇博文,总结一下, ...