Thụ trạng số tổ bộ quyền giá trị đoạn thẳng thụ

#include&LTalgorithm>
#include&LTcstdio>
#include&LTcstring>
#include&LTmap>
#include&LTset>
#define LC o
#define RC o
usingnamespacestd;
constintmaxn=1000010;
intn,m,a[maxn],u[maxn],x[maxn],l[maxn],r[maxn],k[maxn],cur,cur1,cur2,
q1[maxn],q2[maxn],v[maxn];
charop[maxn];
setint>ST;
mapint,int>mp;

structsegment_tree// phong trang động thái khai điểm quyền giá trị đoạn thẳng thụ
{
intcur,rt[maxn*4],sum[maxn*60],lc[maxn*60],rc[maxn*60];

voidbuild(int&o){o=++cur;}

voidprint(into,intl,intr){
if(!o)return;
if(l==r&&sum[o])printf("%d",l);
intmid=(l+r)>>1;
print(lc[o],l,mid);
print(rc[o],mid+1,r);
}

voidupdate(int&o,intl,intr,intx,intv){
if(!o)o=++cur;
sum[o]+=v;
if(l==r)return;
intmid=(l+r)>>1;
if(xmid)
update(lc[o],l,mid,x,v);
else
update(rc[o],mid+1,r,x,v);
}
}st;

// thụ trạng số tổ thực hiện
intlowbit(into){return(o&(-o));}

voidupd(into,intx,intv){
for(;on;o+=lowbit(o))st.update(st.rt[o],1,n,x,v);
}

voidgtv(into,int*A,int&p){
p=0;
for(;o;o-=lowbit(o))A[++p]=st.rt[o];
}

intqry(intl,intr,intk){
if(l==r)returnl;
intmid=(l+r)>>1,siz=0;
for(inti=1;icur1;i++)siz+=st.sum[st.lc[q1[i]]];
for(inti=1;icur2;i++)siz-=st.sum[st.lc[q2[i]]];
// printf( "j %d %d %d %d\n",cur1,cur2,siz,k);
if(siz>=k){
for(inti=1;icur1;i++)q1[i]=st.lc[q1[i]];
for(inti=1;icur2;i++)q2[i]=st.lc[q2[i]];
returnqry(l,mid,k);
}else{
for(inti=1;icur1;i++)q1[i]=st.rc[q1[i]];
for(inti=1;icur2;i++)q2[i]=st.rc[q2[i]];
returnqry(mid+1,r,k-siz);
}
}

/* đoạn thẳng thụ thực hiện
void build(int o,int l,int r)
{
st.build(st.rt[o]);
if(l==r)return;
int mid=(l+r)>>1;
build(LC,l,mid);
build(RC,mid+1,r);
}
void print(int o,int l,int r)
{
printf( "%d %d:",l,r);
st.print(st.rt[o],1,n);
printf( "\n" );
if(l==r)return;
int mid=(l+r)>>1;
print(LC,l,mid);
print(RC,mid+1,r);
}
void update(int o,int l,int r,int q,int x,int v)
{
st.update(st.rt[o],1,n,x,v);
if(l==r)return;
int mid=(l+r)>>1;
if(q
else update(RC,mid+1,r,q,x,v);
}
void getval(int o,int l,int r,int ql,int qr)
{
if(l>qr||r&LTql)return;
if(ql
int mid=(l+r)>>1;
getval(LC,l,mid,ql,qr);
getval(RC,mid+1,r,ql,qr);
}
int query(int l,int r,int k)
{
if(l==r)return l;
int mid=(l+r)>>1,siz=0;
for(int i=1;i
if(siz>=k)
{
for(int i=1;i
return query(l,mid,k);
}
else
{
for(int i=1;i
return query(mid+1,r,k-siz);
}
}
*/

intmain(){
scanf("%d%d",&n,&m);
for(inti=1;in;i++)scanf("%d",a+i),ST.insert(a[i]);
for(inti=1;im;i++){
scanf("%c",op+i);
if(op[i]=='C')
scanf("%d%d",u+i,x+i),ST.insert(x[i]);
else
scanf("%d%d%d",l+i,r+i,k+i);
}
for(setint>::iteratorit=ST.begin();it!=ST.end();it++)
mp[*it]=++cur,v[cur]=*it;
for(inti=1;in;i++)a[i]=mp[a[i]];
for(inti=1;im;i++)
if(op[i]=='C')x[i]=mp[x[i]];
n+=m;
// build(1,1,n);
for(inti=1;in;i++)upd(i,a[i],1);
// print(1,1,n);
for(inti=1;im;i++){
if(op[i]=='C'){
upd(u[i],a[u[i]],-1);
upd(u[i],x[i],1);
a[u[i]]=x[i];
}else{
gtv(r[i],q1,cur1);
gtv(l[i]-1,q2,cur2);
printf("%d\n",v[qry(1,n,k[i])]);
}
}
return0;
}