思路：数据结构
提交：Inf次
题解：
树状数组套主席树
考虑静态区间第k大是一个前缀主席树，但是如果修改是 (O(nlogn)) 的，查询时 (O(logn)) ，考虑去均衡两部分的复杂度，如何均衡的维护前缀和？于是上了树状数组。于是乎主席树 (i) 维护的是 ([i-lowbit(i)+1,i]) 的区间信息。这样修改时 (O(log^2n)) ,查询也是 (O(log^2n))。
// luogu-judger-enable-o2
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#define R register int
const int N=100010,Inf=0x3f3f3f3f;
using namespace std;
inline int g() {
R ret=0,fix=1; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-1:fix;
do ret=ret*10+(ch^48); while(isdigit(ch=getchar())); return ret*fix;
}
int n,m,mxn,tot,cnt;
int a[N],add[N<<8],ll[N<<8],rr[N<<8],sum[N<<8],c[N],LL[30],RR[30],rt[N];
struct que{int k,x,y,z;que() {}
que(int kk,int xx,int yy,int zz) {k=kk,x=xx,y=yy,z=zz;}
} q[N];
inline int lbt(int x) {return x&-x;}
inline void change(int& tr,int l,int r,int pos,int v) {
if(!tr) tr=++tot; sum[tr]+=v; if(l==r) return ; R md=l+r>>1;
pos<=md?change(ll[tr],l,md,pos,v):change(rr[tr],md+1,r,pos,v);
}
inline void calc(int pos,int v,int d) {for(;pos<=n;pos+=lbt(pos)) change(rt[pos],1,mxn,v,d);}
inline int query(int l,int r,int k) {
if(l==r) return l; R md=l+r>>1,ans=0;
for(R i=1;i<=LL[0];++i) ans-=sum[ll[LL[i]]];
for(R i=1;i<=RR[0];++i) ans+=sum[ll[RR[i]]];
if(ans>=k) {
for(R i=1;i<=LL[0];++i) LL[i]=ll[LL[i]];
for(R i=1;i<=RR[0];++i) RR[i]=ll[RR[i]];
return query(l,md,k);
} else {
for(R i=1;i<=LL[0];++i) LL[i]=rr[LL[i]];
for(R i=1;i<=RR[0];++i) RR[i]=rr[RR[i]];
return query(md+1,r,k-ans);
}
}
inline int ask(int l,int r,int k) {
--l; LL[0]=RR[0]=0; memset(LL,0,sizeof(LL)),memset(RR,0,sizeof(RR));
for(;l;l-=lbt(l)) LL[++LL[0]]=rt[l];
for(;r;r-=lbt(r)) RR[++RR[0]]=rt[r];
return query(1,mxn,k);
}
signed main() {
n=g(),m=g(); for(R i=1;i<=n;++i) a[i]=g(),add[++cnt]=a[i];
for(R i=1;i<=m;++i) { register char ch;
while(!isalpha(ch=getchar())); R x=g(),y=g(),z;
if(ch=='Q') z=g(),q[i]=que(1,x,y,z);
else q[i]=que(0,x,y,0),add[++cnt]=y;
} sort(add+1,add+cnt+1); mxn=unique(add+1,add+cnt+1)-add-1;
for(R i=1;i<=n;++i)
a[i]=lower_bound(add+1,add+mxn+1,a[i])-add,calc(i,a[i],1);
for(R i=1;i<=m;++i) {
R x=q[i].x,y=q[i].y,z;
if(q[i].k) z=q[i].z,printf("%d
",add[ask(x,y,z)]);
else calc(x,a[x],-1),a[x]=lower_bound(add+1,add+mxn+1,y)-add,calc(x,a[x],1);
} //system("pause");
}
权值线段树套区间线段树
然而他T了（BZOJ上不会T）
外层权值自带二分属性，可以(O(log^2n))查询，(O(log^2n))修改。
#include<bits/stdc++.h>
#define R register int
using namespace std;
namespace Luitaryi {
static char B[1<<22],*S=B,*T=B;
#define getchar() (S==T&&(T=(S=B)+fread(B,1,1<<15,stdin),S==T)?EOF:*S++)
inline int g() { R x=0,f=1;
register char ch; while(!isdigit(ch=getchar())) f=ch=='-'?-1:f;
do x=x*10+(ch^48); while(isdigit(ch=getchar())); return x*f;
} const int N=100010;
int n,m,a[N],b[N<<1],cnt;
namespace Int{
const int M=N*17*16; int tot,top,stk[N*17];
struct node {int ls,rs,sum;} t[M];
#define ls t[tr].ls
#define rs t[tr].rs
#define sum(tr) t[tr].sum
inline void change(int& tr,int l,int r,int p,int d) {
if(!tr) {if(top) tr=stk[top--]; else tr=++tot;}
if(l==r) {sum(tr)+=d; if(!sum(tr)) stk[++top]=tr,tr=0; return ;}
R md=l+r>>1; if(p<=md) change(ls,l,md,p,d); else change(rs,md+1,r,p,d);
sum(tr)=sum(ls)+sum(rs);
}
inline int query(int tr,int l,int r,int LL,int RR) {
if(LL<=l&&r<=RR) return sum(tr); R md=l+r>>1,ret=0;
if(ls&&LL<=md) ret+=query(ls,l,md,LL,RR);
if(rs&&RR>md) ret+=query(rs,md+1,r,LL,RR); return ret;
}
#undef ls
#undef rs
}
#define ls (tr<<1)
#define rs (tr<<1|1)
int rt[N<<2];
struct node { int l,r,k; node() {}
node(int _l,int _r,int _k) {l=_l,r=_r,k=_k;}
}c[N];
int p,vl,d;
inline void change(int tr,int l,int r) {
Int::change(rt[tr],1,n,p,d); if(l==r) return ; R md=l+r>>1;
if(vl<=md) change(ls,l,md); if(vl>md) change(rs,md+1,r);
}
int LL,RR;
inline int query(int tr,int l,int r,int k) {
if(l==r) return l; R md=l+r>>1,tmp=0;
if(rt[ls]) tmp=Int::query(rt[ls],1,n,LL,RR);
if(tmp<k) return query(rs,md+1,r,k-tmp);
return query(ls,l,md,k);
}
inline void main() { //freopen("in.in","r",stdin); freopen("out.out","w",stdout);
n=g(),m=g(); for(R i=1;i<=n;++i) b[++cnt]=a[i]=g();
for(R i=1,l,r,k;i<=m;++i) {
register char ch; while(!isalpha(ch=getchar())); l=g(),r=g();
if(ch=='Q') k=g(),c[i]=node(l,r,k);
else c[i]=node(l,r,0),b[++cnt]=r;
} sort(b+1,b+cnt+1),cnt=unique(b+1,b+cnt+1)-b-1;
for(R i=1;i<=n;++i) a[i]=lower_bound(b+1,b+cnt+1,a[i])-b;
for(R i=1;i<=m;++i) if(!c[i].k) c[i].r=lower_bound(b+1,b+cnt+1,c[i].r)-b;
d=1; for(R i=1;i<=n;++i) p=i,vl=a[i],change(1,1,cnt);
for(R i=1;i<=m;++i) {
if(c[i].k) LL=c[i].l,RR=c[i].r,printf("%d
",b[query(1,1,cnt,c[i].k)]);
else { p=c[i].l,vl=a[p],d=-1;
change(1,1,cnt),vl=a[p]=c[i].r,d=1,change(1,1,cnt);
}
}
}
} signed main() {Luitaryi::main(); return 0;}
整体二分
值域上整体二分，落实操作区间 ([LL,RR]) 中小于 (md) 的修改操作，扔到左边，大于 (md) 的修改忽略并扔到右边；并落实 ([LL,RR]) 中所有查询操作，若 (k) 比当前自己的排名大，那就把 (k) 剪掉当前的贡献，扔到右边，否则扔到左边。递归左边和右边。复杂度(O(nlog^2n))
#include<bits/stdc++.h>
#define R register int
using namespace std;
namespace Luitaryi {
inline int g() { R x=0,f=1;
register char ch; while(!isdigit(ch=getchar())) f=ch=='-'?-1:f;
do x=x*10+(ch^48); while(isdigit(ch=getchar())); return x*f;
} const int N=100010,Inf=1e9;
int n,m,cnt,tot,ans[N],b[N],c[N];
struct node {int l,r,k,w,op; node() {}
node(int _l,int _r,int _k,int _w,int _op) {l=_l,r=_r,k=_k,w=_w,op=_op;}
} a[N*3],mem1[3*N],mem2[3*N];
inline void add(int p,int d) {for(;p<=n;p+=p&-p) c[p]+=d;}
inline int query(int p) { R ret=0;
for(;p;p-=p&-p) ret+=c[p]; return ret;
}
inline void solve(int l,int r,int LL,int RR) {
if(LL>RR) return ;
if(l==r) {for(R i=LL;i<=RR;++i) if(!a[i].op) ans[a[i].w]=l; return ;}
R md=l+r>>1,p=0,q=0;
for(R i=LL,tmp;i<=RR;++i) {
if(a[i].op) {
if(a[i].r<=md) add(a[i].l,a[i].w),mem1[++p]=a[i];
else mem2[++q]=a[i];
} else {
tmp=query(a[i].r)-query(a[i].l-1);
if(tmp<a[i].k) a[i].k-=tmp,mem2[++q]=a[i];
else mem1[++p]=a[i];
}
} for(R i=1;i<=p;++i) if(mem1[i].op) add(mem1[i].l,-mem1[i].w);
for(R i=1;i<=p;++i) a[LL-1+i]=mem1[i];
for(R i=1;i<=q;++i) a[LL-1+p+i]=mem2[i];
solve(l,md,LL,LL+p-1),solve(md+1,r,LL+p,RR);
}
inline void main() {
n=g(),m=g(); for(R i=1;i<=n;++i) b[i]=g(),a[++cnt]=node(i,b[i],0,1,1);
for(R i=1,l,r,k;i<=m;++i) { register char ch; while(!isalpha(ch=getchar()));
if(ch=='Q') l=g(),r=g(),k=g(),a[++cnt]=node(l,r,k,++tot,0);
else l=g(),r=g(),a[++cnt]=node(l,b[l],0,-1,1),a[++cnt]=node(l,r,0,1,1),b[l]=r;
} solve(1,Inf,1,cnt); for(R i=1;i<=tot;++i) printf("%d
",ans[i]);
}
} signed main() {Luitaryi::main(); return 0;}
2019.09.14
62