BZOJ3065(替罪羊树套线段树)
以前看到这题都瑟瑟发抖,终于过了心情舒畅。
按下标为关键字建替罪羊树,每个结点开一个权值线段树,维护的这个结点代表的子树的信息。
这题还得垃圾回收,自己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(替罪羊树套线段树)的更多相关文章
- [BZOJ 1901] Dynamic Rankings 【树状数组套线段树 || 线段树套线段树】
题目链接:BZOJ - 1901 题目分析 树状数组套线段树或线段树套线段树都可以解决这道题. 第一层是区间,第二层是权值. 空间复杂度和时间复杂度均为 O(n log^2 n). 线段树比树状数组麻 ...
- ZJOI 2017 树状数组(线段树套线段树)
题意 http://uoj.ac/problem/291 思路 不难发现,九条カレン醬所写的树状数组,在查询区间 \([1,r]\) 的时候,其实在查询后缀 \([r,n]\) :在查询 \([l,r ...
- BZOJ4317Atm的树&BZOJ2051A Problem For Fun&BZOJ2117[2010国家集训队]Crash的旅游计划——二分答案+动态点分治(点分树套线段树/点分树+vector)
题目描述 Atm有一段时间在虐qtree的题目,于是,他满脑子都是tree,tree,tree…… 于是,一天晚上他梦到自己被关在了一个有根树中,每条路径都有边权,一个神秘的声音告诉他,每个点到其他的 ...
- hdu-4819-线段树套线段树
http://acm.hdu.edu.cn/showproblem.php?pid=4819 给出一个N*N的矩阵,每次询问一个m*m的子矩阵里的floor((maxv+minv)/2)并把中间的元素 ...
- dfs序+主席树 或者 树链剖分+主席树(没写) 或者 线段树套线段树 或者 线段树套splay 或者 线段树套树状数组 bzoj 4448
4448: [Scoi2015]情报传递 Time Limit: 20 Sec Memory Limit: 256 MBSubmit: 588 Solved: 308[Submit][Status ...
- 【vijos】1750 建房子(线段树套线段树+前缀和)
https://vijos.org/p/1750 是不是我想复杂了.... 自己yy了个二维线段树,然后愉快的敲打. 但是wa了两法.......sad 原因是在处理第二维的更新出现了个小问题,sad ...
- [TJOI2017][bzoj4889] 不勤劳的图书管理员 [线段树套线段树]
题面 传送门 思路 考虑两本书的位置交换对答案的贡献: (为了方便描述,用"左边那本"和"右边那本"称呼两本我们要交换的书,"中间那本"是我 ...
- 【bzoj4785】[Zjoi2017]树状数组 线段树套线段树
题目描述 漆黑的晚上,九条可怜躺在床上辗转反侧.难以入眠的她想起了若干年前她的一次悲惨的OI 比赛经历.那是一道基础的树状数组题.给出一个长度为 n 的数组 A,初始值都为 0,接下来进行 m 次操作 ...
- BZOJ1901 ZOJ2112 线段树+treap (线段树套线段树)
BZOJ1901: 线段树套线段树做法: (外层线段树 里层动态开节点的权值线段树) 有一个小小的trick 可以省掉二分变成nlog^2n的 就是把查询的区间都取出来- logn个一起走- 2016 ...
- [BZOJ 3110] [luogu 3332] [ZJOI 2013]k大数查询(权值线段树套线段树)
[BZOJ 3110] [luogu 3332] [ZJOI 2013]k大数查询(权值线段树套线段树) 题面 原题面有点歧义,不过从样例可以看出来真正的意思 有n个位置,每个位置可以看做一个集合. ...
随机推荐
- C简单实现双向链表
<pre name="code" class="cpp">//链表结构 typedef struct DulNode { DataType data ...
- 【iOS】swift-ObjectC 在iOS 8中使用UIAlertController
iOS 8的新特性之一就是让接口更有适应性.更灵活,因此许多视图控制器的实现方式发生了巨大的变化.全新的UIPresentationController在实现视图控制器间的过渡动画效果和自适应设备尺寸 ...
- github感悟
刚学GitHub进入网页全英文的,感觉很惊讶,自己竟然要在全英文的网站上学习,虽然是英文的但并不感觉有压力,可能之前用eclipse就是全英文的现在除了惊讶,没太多的想法.然后就是我的GitHub地址 ...
- day-3 python多线程编程知识点汇总
python语言以容易入门,适合应用开发,编程简洁,第三方库多等等诸多优点,并吸引广大编程爱好者.但是也存在一个被熟知的性能瓶颈:python解释器引入GIL锁以后,多CPU场景下,也不再是并行方式运 ...
- Django-rest-framework源码分析----权限
添加权限 (1)API/utils文件夹下新建premission.py文件,代码如下: message是当没有权限时,提示的信息 # utils/permission.py class SVIPPr ...
- SpringBoot14 SpringBoot整合mybatis
1 版本说明 springboot:2.0 jdk:1.8 2 创建springBoot项目 创建项目时勾选必要web,MySQL,mybatis相关依赖 创建完成后再pom文件中添加自动部署.lom ...
- Mego开发文档 - 索引
Mego 开发文档 Mego 快速概述 主要特性 获取Mego 使用流程 模型 查询 保存数据 入门 Mego 快速开始 创建项目 安装Nuget包 创建连接字符串 创建模型及数据上下文(添加引用) ...
- MySQL中使用sql语句获得表结构
最近在研究PHP,那么就必须涉及到mysql.其中一个功能通过表数据自动生成页面,紧接着发现在一张空表中无法读取数据(因为人家刚刚新建,就是空的没有数据) 延伸出来便是直接查表结构获得字段名,再进行处 ...
- RSA的公钥、私钥
一.举个例子 1.发消息 用对方的公钥给对方发消息 2.发公告 发公告的时候,用自己的私钥形成签名! 二.加密和签名 RSA的公钥.私钥是互相对应的,RSA会生成两个密钥,你可以把任何一个用于公钥,然 ...
- 老帖收藏,留供参考:SpringMvc2.5+Mybatis3.2.7
一.项目背景 SpringMvc+Mybatis 数据库连接池是阿里巴巴的druid.日志框架式logback 二.配置文件 1.SpringMvc-servlet.xml <?xml vers ...