整体二分之静态区间第K大数

对于整体二分,其他的帖子讲的更详细,而我比较煞笔不看代码根本不知道是什么意思,所以留下代码以供自己后续理解。
我的理解是,整体二分是对于所有询问区间的二分。
我们先假设一个值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]);
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值