平方分割

本文介绍了一种高效的区间第K小数查询算法,利用预处理和分桶技术提高查询效率。通过将数据分为多个桶,并对每个桶进行部分排序,可以在O(logN)的时间复杂度内完成查询。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

#include<cstdio>
#include<algorithm>
#include<vector>
#include<iostream>
using namespace std;
#define maxn 100005
int n;
int m;
const int B=1000;
int a[maxn];
int I[maxn],J[maxn],K[maxn];
int num[maxn];
vector<int>bucket[maxn/B];
void solve(){
for(int i=0;i<n;i++){  
        bucket[i/B].push_back(a[i]);  
        num[i]=a[i];  
    }  
	sort(num,num+n);
	  
    for(int i=0;i<n/B;i++){  
        sort(bucket[i].begin(),bucket[i].end());  
    }  
	for(int i=0;i<m;i++){
	int l=I[i],r=J[i],k=K[i];  
        int lb=-1,ub=n-1;  
		while(ub-lb>1){
			int md=(lb+ub)/2;
			int x=num[md];
			int  tl=l-1,tr=r;
			int c=0;
			while(tl<tr&&tl%B!=0)if(a[tl++]<=x)c++;
			while(tl<tr&&tr%B!=0)if(a[--tr]<=x)c++;
			while(tl<tr){
				int b=tl/B;
				c+=upper_bound(bucket[b].begin(),bucket[b].end(),x)
				-bucket[b].begin();
				tl+=B;
			}
		if(c>=k)ub=md;
		else lb=md;
		}
		cout<<num[ub]<<endl;
	}
}

int main(){
		cin>>n>>m;
for(int  i=0;i<n;i++){
	cin>>a[i];
}	

for(int i=0;i<m;i++){
	cin>>I[i]>>J[i]>>K[i];
	
}
solve();
return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值