ZOJ - 2112 主席树套树状数组
题意:动态第k大,可单点更新,操作+原数组范围6e4
年轻人的第一道纯手工树套树
静态第k大可以很轻易的用权值主席树作差而得
而动态第k大由于修改第i个数会影响[i...n]棵树,因此我们不能在原主席树T上扩展,
而是另开一个表示影响的主席树dT并由树状数组来控制更新和查询前缀线段树和
注意更新时dT的下标要遵循树状数组的子集表示方法
查询时留意整个dT下移到儿子的操作
最后祝ZOJ全家Segmentation Fault
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<map>
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define rrep(i,j,k) for(int i=j;i>=k;i--)
#define erep(i,u) for(int i=head[u];~i;i=nxt[i])
#define print(a) printf("%lld",(ll)(a))
#define printbk(a) printf("%lld ",(ll)(a))
#define println(a) printf("%lld\n",(ll)(a))
using namespace std;
const int MAXN = 6e4+11;
const int MAXM = MAXN<<5;
const int MAXT = MAXN; //!!!!
typedef long long ll;
ll read(){
ll x=0,f=1;register char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
int a[MAXN],b[MAXN],c[MAXN],Q[10006][3],n,m,nn,N;
map<int,int> mp;
char str[1<<5];
const int FLAG = 1e9+7;
struct FST{
int num[MAXM];
int lc[MAXM],rc[MAXM];
int T[MAXN],dT[MAXN],tot;
void init(){
memset(T,0,sizeof T);
memset(dT,0,sizeof dT);
memset(lc,0,sizeof lc);
memset(rc,0,sizeof rc);
memset(num,0,sizeof num);
tot=0;
}
int build(int l,int r){
int cur=++tot;
lc[cur]=rc[cur]=num[cur]=0;
if(l==r) return cur;
int mid=l+r>>1;
lc[cur]=build(l,mid);
rc[cur]=build(mid+1,r);
num[cur]=num[lc[cur]]+num[rc[cur]];
return cur;
}
inline void copy(int cur,int old){
num[cur]=num[old];
lc[cur]=lc[old];
rc[cur]=rc[old];
}
int update(int old,int l,int r,int k,int v){
int cur=++tot;
copy(cur,old);
num[cur]+=v;
if(l==r) return cur;
int mid=l+r>>1;
if(k<=mid) lc[cur]=update(lc[old],l,mid,k,v);
else rc[cur]=update(rc[old],mid+1,r,k,v);
return cur;
}
}fst;
struct FT{
int fstid[MAXT][2];
void init(){
memset(fstid,0,sizeof fstid);
}
inline int lowbit(int x){
return x&-x;
}
void update(int k,int v){
int t=lower_bound(b+1,b+1+N,a[k])-b;
for(int i=k;i<MAXT;i+=lowbit(i)){
fst.dT[i]=fst.update(fst.dT[i],1,N,t,v);
}
}
int query(int k,bool who){
int res=0;
for(int i=k;i>0;i-=lowbit(i)){
res+=fst.num[fst.lc[fstid[i][who]]];
}
return res;
}
}ft;
int query(int st,int ed,int L,int R,int l,int r,int k){
for(int i=L;i>0;i-=ft.lowbit(i)) ft.fstid[i][0]=fst.dT[i];
for(int i=R;i>0;i-=ft.lowbit(i)) ft.fstid[i][1]=fst.dT[i];
while(1){
if(l==r) return l;
int t=fst.num[fst.lc[ed]]-fst.num[fst.lc[st]];
t+=ft.query(R,1);
t-=ft.query(L,0);
if(k<=t){
ed=fst.lc[ed];
st=fst.lc[st];
for(int i=L;i>0;i-=ft.lowbit(i)) ft.fstid[i][0]=fst.lc[ft.fstid[i][0]];
for(int i=R;i>0;i-=ft.lowbit(i)) ft.fstid[i][1]=fst.lc[ft.fstid[i][1]];
r=l+r>>1;
}else{
k-=t;
ed=fst.rc[ed];
st=fst.rc[st];
for(int i=L;i>0;i-=ft.lowbit(i)) ft.fstid[i][0]=fst.rc[ft.fstid[i][0]];
for(int i=R;i>0;i-=ft.lowbit(i)) ft.fstid[i][1]=fst.rc[ft.fstid[i][1]];
l=l+r>>1;l++;
}
}
}
int main(){
//freopen("t.in","r",stdin);
//freopen("t.out","w",stdout);
int T=read();
while(T--){
n=read(); m=read();
mp.clear(); nn=n; ft.init();
rep(i,1,n) b[i]=a[i]=read();
rep(i,1,m){
scanf("%s",str);
Q[i][2]=FLAG;
Q[i][0]=read();
Q[i][1]=read();
if(str[0]=='Q') Q[i][2]=read();
else{
++nn;
b[nn]=a[nn]=Q[i][1];
}
}
sort(b+1,b+1+nn);
N=unique(b+1,b+1+nn)-b-1;
rep(i,1,nn) c[i]=lower_bound(b+1,b+1+N,a[i])-b;
rep(i,1,nn) mp[c[i]]=a[i];
fst.init(); fst.T[0]=fst.build(1,N);
rep(i,1,n) fst.T[i]=fst.update(fst.T[i-1],1,N,c[i],1);
rep(i,1,n) fst.copy(fst.dT[i]=++fst.tot,fst.T[0]);
rep(i,1,m){
if(Q[i][2]!=FLAG){
int L=Q[i][0];
int R=Q[i][1];
int k=Q[i][2];
println(mp[query(fst.T[L-1],fst.T[R],L-1,R,1,N,k)]);
}else{
int k=Q[i][0];
int v=Q[i][1];
ft.update(k,-1);
a[k]=v;
ft.update(k,1);
}
}
}
return 0;
}
ZOJ - 2112 主席树套树状数组的更多相关文章
- [luogu2617][bzoj1901][Zju2112]Dynamic Rankings【树套树+树状数组+主席树】
题目网址 [传送门] 题目大意 请你设计一个数据结构,支持单点修改,区间查询排名k. 感想(以下省略脏话inf个字) 真的强力吹爆洛谷数据,一般的树套树还给我T了一般的点,加强的待修主席树还给我卡了几 ...
- [bzoj3196][Tyvj1730]二逼平衡树_树套树_位置线段树套非旋转Treap/树状数组套主席树/权值线段树套位置线段树
二逼平衡树 bzoj-3196 Tyvj-1730 题目大意:请写出一个维护序列的数据结构支持:查询给定权值排名:查询区间k小值:单点修改:查询区间内定值前驱:查询区间内定值后继. 注释:$1\le ...
- UOJ#291. 【ZJOI2017】树状数组 树套树
原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ291.html 题解 结论:这个写错的树状数组支持的是后缀加和后缀求和.这里的后缀求和在 x = 0 的时 ...
- Codechef EDGEST 树套树 树状数组 线段树 LCA 卡常
原文链接http://www.cnblogs.com/zhouzhendong/p/9016579.html 题目传送门 - Codechef EDGEST 题意 给定相同点集上的两棵生成树$T_1$ ...
- BZOJ3110 [Zjoi2013]K大数查询 树套树 线段树 整体二分 树状数组
欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ3110 题意概括 有N个位置,M个操作.操作有两种,每次操作如果是1 a b c的形式表示在第a个位 ...
- BZOJ1452 [JSOI2009]Count 【树套树 (树状数组)】
1452: [JSOI2009]Count Time Limit: 10 Sec Memory Limit: 64 MB Submit: 2693 Solved: 1574 [Submit][St ...
- [Bzoj3262]陌上花开(CDQ分治&&树状数组||树套树)
题目链接 题目就是赤裸裸的三维偏序,所以用CDQ+树状数组可以比较轻松的解决,但是还是树套树好想QAQ CDQ+树状数组 #include<bits/stdc++.h> using nam ...
- 2019南昌网络赛-I. Yukino With Subinterval 线段树套树状数组,CDQ分治
TMD...这题卡内存卡的真优秀... 所以以后还是别用主席树的写法...不然怎么死的都不知道... 树套树中,主席树方法开权值线段树...会造成空间的浪费...这道题内存卡的很紧... 由于树套树已 ...
- 洛谷P3380 【模板】二逼平衡树(树套树)(线段树+树状数组)
P3380 [模板]二逼平衡树(树套树) 题目描述 您需要写一种数据结构(可参考题目标题),来维护一个有序数列,其中需要提供以下操作: 查询k在区间内的排名 查询区间内排名为k的值 修改某一位值上的数 ...
随机推荐
- Oracle——约束
NOT NULLUNIQUE PRIMARY KEYFOREIGN KEYCHECK 如果不指定约束名 ,Oracle server 自动按照 SYS_Cn 的格式指定约束名 --指定约束名 CREA ...
- Photo2
Story: 驯鹿:“其实我只是想要一个肩膀而已.” 小男孩:“当你需要我的时候,我会在你身边.” Profession: 页面的主色调是淡黄色,这种柔和的色调表达出了柔和的气氛,整个画面颜色的运用都 ...
- js中的函数参数问题
js函数没有Java中的重载现象.js函数的参数是放在arguments的容器里面的. <script type="text/javascript"> functio ...
- 【转】ConcurrentMap 分析和思考
预备知识:Java HashMap and HashSet 的实现机制 由预备知识可以知道hashmap 的存储结构为: (图像来自http://www.ibm.com/developerworks/ ...
- 敏捷软件开发:原则、模式与实践——第11章 DIP:依赖倒置原则
第11章 DIP:依赖倒置原则 DIP:依赖倒置原则: a.高层模块不应该依赖于低层模块.二者都应该依赖于抽象. b.抽象不应该依赖于细节.细节应该依赖于抽象. 11.1 层次化 下图展示了一个简单的 ...
- 编写高质量代码改善C#程序的157个建议——建议108:将类型标识为sealed
建议108:将类型标识为sealed sealed能够阻止类型被其他类型继承.代码如下: sealed class SampleClass { } class OtherClass : SampleC ...
- python文件操作os模块
Python 统计某一文件夹下文件数量 使用python pathlib模块 from pathlib import Path dir_path = ' ' print(len(list(Path( ...
- C#操作系统计划任务
首先需要引用system32下的taskschd.dll 测试环境:win8+vs2012+.net4.0 /// <summary> /// 删除任务 /// </summary& ...
- .Net Core使用OpenXML导出,导入Excel
导出Excel是程序很常用到的功能,.Net Core可以借助Open-XML-SDK来导出Excel. Open-XML-SDK open-xml-sdk是是微软开源的项目.Open XML SDK ...
- solr7.4 centos7安装
环境:centos7.JDK1.8.solr 自带Jetty启动 一.安装JDK1.8环境 1.下载JDK jdk-8u172-linux-x64.rpm 下载地址:http://www.oracle ...