对于整体二分,其他的帖子讲的更详细,而我比较煞笔不看代码根本不知道是什么意思,所以留下代码以供自己后续理解。
我的理解是,整体二分是对于所有询问区间的二分。
我们先假设一个值mid,如果给出的数列,它比mid要小那么就将它的id(或者说是位置)以1的形式加入到树状数组(表示这个位置内存在这个数),反之则从树状数组中删掉。这样一来,树状数组就可以查询当前询问区间内存在多少个数,如果当前询问区间存在的数的个数比其要询问的第K大数要多(即意味着第K大数必然再当前询问区间内),那么将该询问区间划分到二分的左边,二分的右边则表示个数比第K大数要小(意味着第K大数还不在这个区间内)。
然后不断的二分假设的值和二分左右两边询问区间,最终,假设的值终究会符合其中一些询问区间的要求,也就是说这个假设的值就是一些询问期间的第K大数。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<algorithm>
#include<bitset>
#define scand(a) scanf("%d",&a)
#define scandd(a,b) scanf("%d%d",&a,&b)
#define scanddd(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define mst(a,b) memset(a,b,sizeof(a))
#define lowbit(x) x&-x
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int maxn=1e5+5;
struct nod
{
int x,id;
}a[maxn];
struct query
{
int x,y,k,id;
}que[maxn],tmp[maxn];
int n,m,ans[maxn],cnt,tre[maxn];
int sum(int x)
{
int sum=0;
while(x>0)
{
sum+=tre[x];
x-=lowbit(x);
}
return sum;
}
void add(int x,int k)
{
while(x<=n)
{
tre[x]+=k;
x+=lowbit(x);
}
}
bool cmp(nod x,nod y)
{
return x.x<y.x;
}
void ef(int l,int r,int ql,int qr)
{
if(ql>qr)return;
int p=ql,q=qr,mid=l+r>>1;
if(l==r)//从-inf到inf不断二分,当l==r时,说明此时这个值符合区间询问的答案
{
for(int i=ql;i<=qr;i++)
ans[que[i].id]=l;
return;
}
while(a[cnt+1].x<=mid&&cnt<n)
{
//对于当前区间[l,r],如果询问区间左端的数小于当前区间的中间值,那么就将其加进来
add(a[++cnt].id,1);
}
while(cnt&&a[cnt].x>mid)
{
add(a[cnt--].id,-1);
}
for(int i=ql;i<=qr;i++)
{
//然后通过树状数组来查询当前询问区间右端到左端这一部分,已经有多少个数了
int all=sum(que[i].y)-sum(que[i].x-1);
if(all>=que[i].k)tmp[p++]=que[i];//若这一部分的数大过询问区间的要求的值那么就将p右移
else tmp[q--]=que[i];//反之,q左移,如此一来,当这个for循环跑完时,p和q的位置就是询问区间的中间值
}//也就是说中间值的左端,全部都是询问区间的k值小于询问区间内所拥有的数的数量,中间值右端同理
for(int i=ql;i<=qr;i++)
que[i]=tmp[i];//相当于对其进行划分,大于询问区间中间值的在右边,反之在左边
ef(l,mid,ql,p-1);
//然后就从当前区间左边到中间值去搜索符合询问区间的值(此时的询问区间的范围变为原询问区间的左端和其中间值)
ef(mid+1,r,q+1,qr);
}
int main()
{
#ifdef local
freopen("1.txt","r",stdin);
#endif
scandd(n,m);
for(int i=1;i<=n;i++)
{
scand(a[i].x);
a[i].id=i;
}
sort(a+1,a+1+n,cmp);
for(int i=1;i<=m;i++)
{
scanddd(que[i].x,que[i].y,que[i].k);
que[i].id=i;
}
ef(-inf,inf,1,m);
for(int i=1;i<=m;i++)
printf("%d\n",ans[i]);
}