BZOJ3065:带插入区间K小值
浅谈树状数组与主席树:https://www.cnblogs.com/AKMer/p/9946944.html
题目传送门:https://www.lydsy.com/JudgeOnline/problem.php?id=3065
去%了一发hzwer的博客学会的。
替罪羊树套线段树,外层维护位置,内层维护权值。
时间复杂度:\(O(nlog^2n)\)
空间复杂度:\(O(nlog^2n)\)
代码如下:
#include <cstdio>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
const double alpha=0.75;
const int maxn=7e4+5,maxsz=2e7;
char s[5];
int n,m,lstans;
int read() {
int x=0,f=1;char ch=getchar();
for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;
for(;ch>='0'&&ch<='9';ch=getchar())x=x*10+ch-'0';
return x*f;
}
struct Segment_tree {
int tot;
vector<int>rub;
int sum[maxsz],ls[maxsz],rs[maxsz];
int newnode() {
if(rub.empty())return ++tot;
else {
int res=rub.back();rub.pop_back();
return res;
}
}
void change(int p,int l,int r,int pos,int v) {
while(1) {
sum[p]+=v;
if(l==r)return;
int mid=(l+r)>>1;
if(pos<=mid) {
if(!ls[p])ls[p]=newnode();
p=ls[p],r=mid;
}
else {
if(!rs[p])rs[p]=newnode();
p=rs[p],l=mid+1;
}
}
}
void recycle(int u) {
if(!u)return;
rub.push_back(u);
recycle(ls[u]),recycle(rs[u]);
sum[u]=ls[u]=rs[u]=0;
}
}T2;
struct TiZuiYang_tree {
vector<int>num;
int tot,root,cnt;
int son[maxn][2];
vector<int>::iterator it;
int rt[maxn],siz[maxn],tmp[maxn],a[maxn],id[maxn];
int build(int l,int r) {
if(r<l)return 0;
if(l==r) {
int u=id[l];son[u][0]=son[u][1]=0;
rt[u]=T2.newnode(),siz[u]=1;
T2.change(rt[u],0,70000,a[u],1);
return u;
}
int mid=(l+r)>>1,u=id[mid];
rt[u]=T2.newnode();siz[u]=1;
son[u][0]=build(l,mid-1);
son[u][1]=build(mid+1,r);
siz[u]+=siz[son[u][0]]+siz[son[u][1]];
for(int i=l;i<=r;i++)T2.change(rt[u],0,70000,a[id[i]],1);
return u;
}
void find(int u,int l,int r) {
if(l==1&&r==siz[u]) {tmp[++cnt]=rt[u];return;}
int sz=siz[son[u][0]];
if(l<=sz+1&&r>=sz+1)num.push_back(a[u]);
if(r<=sz)find(son[u][0],l,r);
else if(l>sz+1)find(son[u][1],l-sz-1,r-sz-1);
else {
if(l<=sz)find(son[u][0],l,sz);
if(r>sz+1)find(son[u][1],1,r-sz-1);
}
}
int query(int l,int r,int rk) {
num.clear(),cnt=0,find(root,l,r);
int L=0,R=70000;
while(L!=R) {
int res=0,mid=(L+R)>>1;
for(int i=1;i<=cnt;i++)
res+=T2.sum[T2.ls[tmp[i]]];
for(it=num.begin();it!=num.end();it++)
if(*it>=L&&(*it)<=mid)res++;
if(res>=rk) {
for(int i=1;i<=cnt;i++)
tmp[i]=T2.ls[tmp[i]];
R=mid;
}
else {
for(int i=1;i<=cnt;i++)
tmp[i]=T2.rs[tmp[i]];
L=mid+1;rk-=res;
}
}
return L;
}
int change(int u,int pos,int x) {
T2.change(rt[u],0,70000,x,1);
int sz=siz[son[u][0]],res;
if(sz+1==pos) {res=a[u];a[u]=x;}
else if(sz>=pos)res=change(son[u][0],pos,x);
else res=change(son[u][1],pos-sz-1,x);
T2.change(rt[u],0,70000,res,-1);
return res;
}
void insert(int &u,int rk,int v) {
if(!u) {
u=++tot;rt[u]=T2.newnode();
T2.change(rt[u],0,70000,v,1);
siz[u]=1;a[u]=v;return;
}
T2.change(rt[u],0,70000,v,1);
int sz=siz[son[u][0]];siz[u]++;
if(sz>=rk)insert(son[u][0],rk,v);
else insert(son[u][1],rk-sz-1,v);
}
void dfs_rebuild(int u) {
if(!u)return;
dfs_rebuild(son[u][0]);
T2.recycle(rt[u]);
id[++cnt]=u;rt[u]=0;
dfs_rebuild(son[u][1]);
}
void find_rebuild(int &u,int rk) {
if(alpha*siz[u]<max(siz[son[u][0]],siz[son[u][1]])) {
cnt=0;dfs_rebuild(u);u=build(1,cnt);
return;
}
if(siz[son[u][0]]>=rk)find_rebuild(son[u][0],rk);
if(rk>siz[son[u][0]]+1)find_rebuild(son[u][1],rk-siz[son[u][0]]-1);
}
}T1;
int main() {
T1.tot=n=read();
for(int i=1;i<=n;i++)
T1.id[i]=i,T1.a[i]=read();
T1.root=T1.build(1,n);
m=read();
for(int i=1;i<=m;i++) {
scanf("%s",s+1);
if(s[1]=='Q') {
int l=read()^lstans,r=read()^lstans,k=read()^lstans;
lstans=T1.query(l,r,k);printf("%d\n",lstans);
}
if(s[1]=='M') {
int pos=read()^lstans,x=read()^lstans;
T1.change(T1.root,pos,x);
}
if(s[1]=='I') {
int pos=read()^lstans,x=read()^lstans;
T1.insert(T1.root,pos-1,x);
T1.find_rebuild(T1.root,pos);
}
}
return 0;
}
BZOJ3065:带插入区间K小值的更多相关文章
- [BZOJ3065]带插入区间K小值 解题报告 替罪羊树+值域线段树
刚了一天的题终于切掉了,数据结构题的代码真**难调,这是我做过的第一道树套树题,做完后感觉对树套树都有阴影了......下面写一下做题记录. Portal Gun:[BZOJ3065]带插入区间k小值 ...
- BZOJ3065 带插入区间K小值 || 洛谷P4278
这是一道让我崩溃的题...... 然鹅洛谷上时限被改然后只有20分......好像所有人都被卡了(雾) 由于替罪羊树不是依靠旋转操作而是依靠暴力重构的方式维护树的平衡,所以我们可以考虑使用替罪羊树套区 ...
- 【函数式权值分块】【块状链表】bzoj3065 带插入区间K小值
显然是块状链表的经典题.但是经典做法的复杂度是O(n*sqrt(n)*log^2(n))的,出题人明确说了会卡掉. 于是我们考虑每个块内记录前n个块的权值分块. 查询的时候差分什么的,复杂度就是O(n ...
- bzoj3065: 带插入区间K小值
无聊来写了下 一开始发现树高是O(n)的,然后就MLE了,进去看了下发现没有重构! 看了半天发现调用错了函数 然后进去又发现不满足sz = ch[0]->sz + ch[1]->sz + ...
- 【BZOJ3065】带插入区间K小值 替罪羊树+权值线段树
[BZOJ3065]带插入区间K小值 Description 从前有n只跳蚤排成一行做早操,每只跳蚤都有自己的一个弹跳力a[i].跳蚤国王看着这些跳蚤国欣欣向荣的情景,感到非常高兴.这时跳蚤国王决定理 ...
- bzoj 3065: 带插入区间K小值 替罪羊树 && AC300
3065: 带插入区间K小值 Time Limit: 60 Sec Memory Limit: 512 MBSubmit: 1062 Solved: 253[Submit][Status] Des ...
- 【题解】BZOJ 3065: 带插入区间K小值——替罪羊树套线段树
题目传送门 题解 orz vfk的题解 3065: 带插入区间K小值 系列题解 一 二 三 四 惨 一开始用了一种空间常数很大的方法,每次重构的时候merge两颗线段树,然后无限RE(其实是MLE). ...
- 3065: 带插入区间K小值_树套树_替罪羊树_权值线段树
经过周六一天,周一3个小时的晚自习,周二2个小时的疯狂debug,终于凭借自己切掉了这道树套树题. Code: #include <cstdio> #include <algorit ...
- BZOJ 3065 带插入区间K小值(sag套线段树)
3065: 带插入区间K小值 Time Limit: 60 Sec Memory Limit: 512 MBSubmit: 4696 Solved: 1527[Submit][Status][Di ...
- 【学习笔记】浅析平衡树套线段树 & 带插入区间K小值
常见的树套树 一般来说,在嵌套数据结构中,线段树多被作为外层结构使用. 但线段树毕竟是 静态 的结构,导致了一些不便. 下面是一个难以维护的例子: 带插入区间 \(k\) 小值问题 来源:Luogu ...
随机推荐
- 【特征匹配】SIFT原理与C源代码剖析
相关: KD树+BBF算法解析 SURF原理与源代码解析 SIFT的原理已经有非常多大牛的博客上做了解析,本文重点将以Rob Hess等人用C实现的代码做解析,结合代码SIFT原理会更easy理解.一 ...
- web翻译——插件
很多时候,可能我们web项目中需要的只是机械式的翻译,并不需要什么利用xml或者js json等等实现逼真翻译,那样工作量太大.这时候可能你就需要这几款小工具来帮助你.当然,如果 对翻译或者你的项目外 ...
- JQuery如何获取按键的unicode编码?
$("selector").keyup(function(xxx){ var myEvent = xxx; var code = myEvent.keyCode; alert(co ...
- Tsung 初步介绍安装
tsung是erlang的一个开源的一个压力测试工具,可以测试包括HTTP, WebDAV, Mysql, PostgreSQL, LDAP, and XMPP/Jabber等服务器.针对 HTTP ...
- mysql-mongdb-redis
千万级别:mysql 千万以及亿级别:mongdb
- sublime3 支持 jsx 语法
添加几个插件即可在js中快速写html babel 可以识别React,并高亮显示ES6 command+shift+p -> install package -> babel 使用 在打 ...
- 12 redis之aof日志持久化
Aof 的配置 appendonly no # 是否打开 aof日志功能 appendfsync always # 每1个命令,都立即同步到aof. 安全,速度慢 appendfsync everys ...
- 九度OJ 1077:最大序列和 (DP)
时间限制:1 秒 内存限制:32 兆 特殊判题:否 提交:5600 解决:1637 题目描述: 给出一个整数序列S,其中有N个数,定义其中一个非空连续子序列T中所有数的和为T的"序列和&qu ...
- 【题解】[CF718C Sasha and Array]
[题解]CF718C Sasha and Array 对于我这种喜欢写结构体封装起来的选手这道题真是太对胃了\(hhh\) 一句话题解:直接开一颗线段树的矩阵然后暴力维护还要卡卡常数 我们来把\(2 ...
- 【题解】P2602[JZOI2010]数字计数
[题解][P2602ZJOI2010]数字计数 乍看此题,感觉直接从数字的位上面动手,感觉应该很容易. 但是仔细看数据范围,发现如果不利用计数原理,肯定会超时,考虑数码出现的特征: \(A000\)到 ...