带修主席树（Dynamic Rankings）

原题： https://www.luogu.org/problemnew/show/P2617

题意：

区间第k大，支持单点修改

解析：

还是这张图，下面是普通的主席树：

下面是带修主席树：

请注意，原来的主席树每棵树只有一个版本，即不可修改。
现在每棵树都可以修改，那么每棵树会有多个版本，但是一般来说需要用到的只是最后一个版本。

当需要在第一棵树的位置p加上v时，我不仅需要将第一棵树维护至新的版本（在原来的第一棵树的基础上），还需要将树状数组的对应位置维护（1，2，4，8…）。注意，此时的树不需要从前面一棵树转移，只需要一直这样log维护。

通俗的理解： 原来的主席树是sum[i]=sum[i-1]+a[i]，现在变成while(rt){ans+=tr[rt];rt-=lowbit(rt)}了。那么此时，我每棵树就变成了所管辖范围（lowbit）的区间和了。

更新过程：

int update(int last, int now, int l, int r, int p, int v) {
    siz[now]=siz[last]+v;
    ls[now]=ls[last],rs[now]=rs[last];
    if(l==r)
        return now;
    if(p<=mid)
        ls[now]=update(ls[now],++cnt,l,mid,p,v);
    else
        rs[now]=update(rs[now],++cnt,mid+1,r,p,v);
    siz[now]=siz[ls[now]]+siz[rs[now]];
    return now;
}
void add(int _rt,int idx,int v) {   // 第_rt棵树，区间1-num上的点idx，加上值v
    int tt=_rt;
    while(_rt<=n) {
        rt[_rt]=update(rt[_rt],++cnt,1,num,idx,v);
        _rt+=_rt&-_rt;
    }
}

查询过程： （结合此题版本，其实就是一般的查询，现在变成log棵树的查询，每次变到左儿子则log棵树一起变成左儿子。L和R记录的是哪log棵树，由给出的区间算出。右节点y，左节点x：ans[r]-ans[l-1]就是区间ans）

int query(int l,int r,int k) {
    if(l==r)    // 返回时，l和r停在第k小的位置
        return l;
    int v=0;
    rep(i,1,ctl)v-=siz[ls[L[i]]];
    rep(i,1,ctr)v+=siz[ls[R[i]]];
    if(v>=k) {  // 左边已经有至少k个，多了
        rep(i,1,ctl)L[i]=ls[L[i]];
        rep(i,1,ctr)R[i]=ls[R[i]];
        return query(l,mid,k);
    } else {    // 不够，再往右找
        rep(i,1,ctl)L[i]=rs[L[i]];
        rep(i,1,ctr)R[i]=rs[R[i]];
        return query(mid+1,r,k-v);
    }
}

完整代码：

#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define mid (l+r>>1)

const int maxn=200009;
const int maxe=4e7+9;
int n, q;
bool f[maxn];
int x[maxn],y[maxn],z[maxn];
int a[maxn];
int tmp[maxn],num;
inline int getidx(int v) {
    return lower_bound(tmp+1,tmp+1+num,v)-tmp;
}

int rt[maxn];
int ls[maxe], rs[maxe], siz[maxe], cnt;

int update(int last, int now, int l, int r, int p, int v) {
    siz[now]=siz[last]+v;
    ls[now]=ls[last],rs[now]=rs[last];
    if(l==r)
        return now;
    if(p<=mid)
        ls[now]=update(ls[now],++cnt,l,mid,p,v);
    else
        rs[now]=update(rs[now],++cnt,mid+1,r,p,v);
    siz[now]=siz[ls[now]]+siz[rs[now]];
    return now;
}
void add(int _rt,int idx,int v) {   // 第_rt棵树，区间1-num上的点idx，加上值v
    int tt=_rt;
    while(_rt<=n) {
        rt[_rt]=update(rt[_rt],++cnt,1,num,idx,v);
        _rt+=_rt&-_rt;
    }
}
int L[maxn],R[maxn],ctl,ctr;
int query(int l,int r,int k) {
    if(l==r)    // 返回时，l和r停在第k小的位置
        return l;
    int v=0;
    rep(i,1,ctl)v-=siz[ls[L[i]]];
    rep(i,1,ctr)v+=siz[ls[R[i]]];
    if(v>=k) {  // 左边已经有至少k个，多了
        rep(i,1,ctl)L[i]=ls[L[i]];
        rep(i,1,ctr)R[i]=ls[R[i]];
        return query(l,mid,k);
    } else {    // 不够，再往右找
        rep(i,1,ctl)L[i]=rs[L[i]];
        rep(i,1,ctr)R[i]=rs[R[i]];
        return query(mid+1,r,k-v);
    }
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cin>>n>>q;
    num=0;
    rep(i,1,n)cin>>a[i],tmp[++num]=a[i];
    rep(i,1,q) {
        char ch[2];
        cin>>ch;
        if(ch[0]=='C') {
            f[i]=1;
            cin>>x[i]>>y[i];
            tmp[++num]=y[i];
        } else {
            cin>>x[i]>>y[i]>>z[i];
        }
    }
    sort(tmp+1,tmp+1+num);
    num=unique(tmp+1,tmp+1+num)-tmp-1;
    rep(i,1,n)add(i,getidx(a[i]),1);
    rep(i,1,q) {
        if(f[i]) {
            add(x[i],getidx(a[x[i]]),-1);
            a[x[i]]=y[i];
            add(x[i],getidx(y[i]),1);
        } else {
            ctl=ctr=0;
            for(int j=x[i]-1; j; j-=j&-j)
                L[++ctl]=rt[j];
            for(int j=y[i]; j; j-=j&-j)
                R[++ctr]=rt[j];
            int idx=query(1,num,z[i]);
            cout<<tmp[idx]<<'\n';
        }
    }
}