P3834 （静态区间第k大，整体二分）

题意：

给定长度为n的序列a，m次询问，每次询问给出L,R,k，要求计算[L,R]的第k小数

数据范围：n,m<=2e5，-1e9<=a(i)<=1e9

解法：

整体二分的主要思想：
对需要在值域[l,r]内二分的所有询问，共用一个二分，相比每个询问于单独二分，减少了二分总次数

整体二分是离线算法。

需要注意的点：
二分的值域存在负数，二分需要写mid=(l+r)>>1，而不能是mid=(l+r)/2
原因：
正数二分，[1,6]分为[1,3]和[4,6]，因为(1+6)/2=3
负数二分，[-6,-1]会分为[-6,-3]和[-2,-1]，因为(-6±1)/2=-3
发现负数二分，区间大小并不是均分的，正确的应该是[-6,-4]和[-3,-1]

因为多了一个符号，"向下取整"变成"向上取整"了，所以会错
正确做法是用位运算>>1，这样就对了（负数是对补码右移，而正数的补码就是原码，不一样）

ps：
利用树状数组离线的算法，好像常常会将BIT存值变为BIT存下标。

code：

#include<bits/stdc++.h>
using namespace std;
const int maxm=2e5+10;
struct Q{
    int l,r,k,op,id;
}q[maxm<<1],q1[maxm<<1],q2[maxm<<1];//q1:[l,mid],q2:[mid+1,r]
int cnt;
int ans[maxm];
int n,m;
//BIT
int c[maxm];
int lowbit(int i){
    return i&-i;
}
void add(int i,int t){
    while(i<maxm){
        c[i]+=t,i+=lowbit(i);
    }
}
int ask(int i){
    int ans=0;
    while(i){
        ans+=c[i],i-=lowbit(i);
    }
    return ans;
}
//
void solve(int l,int r,int L,int R){
    if(L>R)return ;
    if(l==r){
        for(int i=L;i<=R;i++){
            if(q[i].op==2){
                ans[q[i].id]=l;
            }
        }
        return ;
    }
    int mid=(l+r)>>1,cnt1=0,cnt2=0;
    for(int i=L;i<=R;i++){
        if(q[i].op==1){//number
            if(q[i].l<=mid){
                add(q[i].id,q[i].r),q1[++cnt1]=q[i];
            }else{
                q2[++cnt2]=q[i];
            }
        }else{//query
            int x=ask(q[i].r)-ask(q[i].l-1);
            if(x>=q[i].k){
                q1[++cnt1]=q[i];
            }else{
                q[i].k-=x;
                q2[++cnt2]=q[i];
            }
        }
    }
    for(int i=1;i<=cnt1;i++){//clear BIT
        if(q1[i].op==1)add(q1[i].id,-q1[i].r);
    }
    for(int i=1;i<=cnt1;i++){
        q[L-1+i]=q1[i];
    }
    for(int i=1;i<=cnt2;i++){
        q[L+cnt1-1+i]=q2[i];
    }
    solve(l,mid,L,L+cnt1-1);
    solve(mid+1,r,L+cnt1,R);
}
signed main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        int x;scanf("%d",&x);
        q[++cnt]={x,1,0,1,i};
    }
    for(int i=1;i<=m;i++){
        int l,r,k;scanf("%d%d%d",&l,&r,&k);
        q[++cnt]={l,r,k,2,i};
    }
    solve(-1e9,1e9,1,cnt);
    for(int i=1;i<=m;i++){
        printf("%d\n",ans[i]);
    }
    return 0;
}

---