区间第k大数

思路是分块，不然得树套树（我是蒟蒻不会）

用分块T2的思路+二分就能求出区间第k大数

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#include<cmath>
using namespace std;
#define ll long long
vector<ll>vec[1004];
ll n,blo,v[30005],m,bl[30005],lazy[1004];
void reset(ll x)
{
    vec[x].clear();
    
    for(ll i=(x-1)*blo+1;i<=min(x*blo,n);i++)
    vec[x].push_back(v[i]);
    sort(vec[x].begin(),vec[x].end());
}
void change(ll l,ll r,ll val)
{
    for(ll i=l;i<=min(blo*bl[l],r);i++)
    {
        v[i]+=val;
    }
    reset(bl[l]);
    if(bl[l]!=bl[r])
    {
        for(ll i=(bl[r]-1)*blo+1;i<=r;i++)
        v[i]+=val;
        reset(bl[r]);
    }
    for(ll i=bl[l]+1;i<=bl[r]-1;i++)
    lazy[i]+=val;
}
ll query(ll l,ll r,ll val)
{
    ll ans=0;
    for(ll i=l;i<=min(r,bl[l]*blo);i++)
    if(v[i]+lazy[bl[l]]<=val)ans++;
    if(bl[l]!=bl[r])
    {
        for(ll i=blo*(bl[r]-1)+1;i<=r;i++)
        if(v[i]+lazy[bl[r]]<=val) ans++;
    }
    for(ll i=bl[l]+1;i<=bl[r]-1;i++)
    ans+=upper_bound(vec[i].begin(),vec[i].end(),val-lazy[i])-vec[i].begin();
    return ans;
    
}
int main()
{
    scanf("%lld %lld",&n,&m);
    blo=sqrt(n);
    for(ll i=1;i<=n;i++)
    {
        scanf("%lld",&v[i]);
        bl[i]=(i-1)/blo+1;
        vec[bl[i]].push_back(v[i]);
    }
    for(ll i=1;i<=bl[n];i++)
    sort(vec[i].begin(),vec[i].end());
    char char_[2];
    //cout<<query(2,6,8)<<endl;
    for(ll i=1;i<=m;i++)
    {
        ll x,y,z;
        scanf("%s %lld %lld %lld",char_,&x,&y,&z);
        if(char_[0]=='C')
        {
            change(x,y,z);
        }else
        {
            ll l=0,r=1e15,mid;
            while(l<r)
            {
                
                
                mid=(l+r)>>1;
                //cout<<l<<' '<<r<<' '<<mid<<endl;
                ll tt=query(x,y,mid);
                //cout<<l<<' '<<r<<' '<<mid<<' '<<tt<<endl;
                printf("l:%lld , r:%lld , mid:%lld , query:%lld\n",l,r,mid,tt);
                if(tt<z)
                {
                    l=mid+1;
                }else r=mid;
            }
            printf("%lld\n",r);
        }
        
    }
    
} 

有点慢，可以先ST表预处理出l,r区间里的最大最小值，进行二分