整体二分

带修改区间第k小，使用整体二分。

#include <iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define lowbit(x) (x&-(x))
using namespace std;
#define MAXN 100005
struct node
{
    int tp,pos,val,ed,cur,addr;
}que[MAXN],q1[MAXN],q2[MAXN];
struct optnode
{
    int opt,pos1,pos2,val;
}opts[MAXN];
int head,tail,cnt;
int ans[MAXN];
int n,nn,num[MAXN],q,jump[MAXN],realval[MAXN];
int maxz;
char opt[3];
int tmp[MAXN];
struct lsh
{
    int val,id;
    bool operator <(lsh &t)const
    {
        return val<t.val;
    }
}num2[MAXN];
int tree[MAXN];
void init(int x)
{
    maxz=0;
    for(int i=1;i<=x;i++)
        num2[i].val=num[i],num2[i].id=i;
    sort(num2+1,num2+x+1);
    for(int i=1;i<=x;i++)
    {
        if(num2[i].val!=num2[i-1].val)
           ++maxz,realval[maxz]=num2[i].val;
        num[num2[i].id]=maxz;

    }
}
void add(int pos,int val)
{
    while(pos<=n)
    {tree[pos]+=val;
    pos+=lowbit(pos);
    }
}
int getsum(int pos)
{
    int tmp=0;
    while(pos)
    {
        tmp+=tree[pos];
        pos-=lowbit(pos);
    }
    return tmp;
}
void divide(int l,int r,int head,int tail)
{
    if(head>tail)return;
    if(l==r)
    {
        for(int i=head;i<=tail;i++)
        {
            if(que[i].tp==3)
                ans[que[i].addr]=l;
        }
        return;
    }
    int mid=(l+r)/2;
    for(int i=head;i<=tail;i++)
    {
        if(que[i].tp==1)if(num[que[i].val]<=mid)add(que[i].pos,1);
        if(que[i].tp==2)if(num[que[i].val]<=mid)add(que[i].pos,-1);
        if(que[i].tp==3)
            tmp[i]=getsum(que[i].ed)-getsum(que[i].pos-1);
    }
    for(int i=head;i<=tail;i++)
        {if(que[i].tp==1)if(num[que[i].val]<=mid)add(que[i].pos,-1);
         if(que[i].tp==2)if(num[que[i].val]<=mid)add(que[i].pos,1);
        }

        int l1=0,r1=0;
    for(int i=head;i<=tail;i++)
    {
        if(que[i].tp==3)
        {
            if(que[i].cur+tmp[i]>=que[i].val)
            q1[l1++]=que[i];
            else
            {
                que[i].cur+=tmp[i];
                q2[r1++]=que[i];
            }
        }
        else
        {
            if(num[que[i].val]<=mid)
                q1[l1++]=que[i];
            else q2[r1++]=que[i];
        }
    }
    for(int i=0;i<l1;i++)que[head+i]=q1[i];
    for(int i=0;i<r1;i++)que[head+l1+i]=q2[i];
    divide(l,mid,head,head+l1-1);
    divide(mid+1,r,head+l1,tail);
}

int main()
{
    int t1,t2,t3;
    scanf("%d%d",&n,&q);
    for(int i=1;i<=n;i++)scanf("%d",&num[i]),jump[i]=i;
    int nn=n;
    for(int i=1;i<=n;i++)
    que[++tail].pos=i,que[tail].tp=1,que[tail].val=i;
    for(int i=1;i<=q;i++)
    {
        scanf("%s",opt);
        if(opt[0]=='Q')
        {
            scanf("%d %d %d",&t1,&t2,&t3);
            que[++tail].tp=3,que[tail].pos=t1,que[tail].ed=t2,que[tail].val=t3;
            que[tail].addr=cnt++;
        }
        else{
            scanf("%d%d",&t1,&t2);
            num[++nn]=t2;
            que[++tail].tp=2,que[tail].pos=t1,que[tail].val=jump[t1];
            que[++tail].tp=1,que[tail].pos=t1,que[tail].val=nn;
            jump[t1]=nn;
        }
    }
    init(nn);
    divide(1,maxz,1,tail);
    for(int i=0;i<cnt;i++)
        printf("%d
",realval[ans[i]]);
    return 0;
}