ZOO

题目大意

M个询问询问区间K小值。

整体二分

我们用K大数查询这道题的思路整体二分。

这道题直接套用K大数查询需要带两个log。

我们把询问l~r拆为两个询问1~r和1~l-1，那么每个询问都变成一个前缀，于是我们想着搞一波线性作法。

观察到对于答案区间为l~r，我们只需要扫描在这个区间内的数。

也就是说对于l~r中的数，<=mid的数继续归到左区间其余归到右区间。

一个区间内的询问按照右端点排序，然后就可以线性做。

总复杂度n log n。

注意本题答案范围比较大故先进行离散化。

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
using namespace std;
const int maxn=100000+10;
struct dong{
int a,cnt,x;
bool p;
};
dong ask[maxn*2];
int h1[maxn*5],h2[maxn*5],h3[maxn*5],la1[maxn*5],la2[maxn*5],la3[maxn*5],n1[maxn*20],n2[maxn*20],n3[maxn*20],g1[maxn*20],g2[maxn*20],g3[maxn*20];
int a[maxn],h[maxn],ans[maxn],rank[maxn],id[maxn*2];
int i,j,k,l,t,n,m,t1,t2,t3;
//h1=dl h2=ld h3=que
void add1(int x,int y){
g1[++t1]=y;
n1[la1[x]]=t1;
la1[x]=t1;
if (!h1[x]) h1[x]=t1;
}
void add2(int x,int y){
g2[++t2]=y;
n2[la2[x]]=t2;
la2[x]=t2;
if (!h2[x]) h2[x]=t2;
}
void add3(int x,int y){
g3[++t3]=y;
n3[la3[x]]=t3;
la3[x]=t3;
if (!h3[x]) h3[x]=t3;
}
void solve(int p,int l,int r){
int j,i,k,t,num=0;
if (l==r){
t=h3[p];
while (t){
j=g3[t];
ans[j]=h[l];
t=n3[t];
}
return;
}
int mid=(l+r+1)/2;
i=h2[p];k=h1[p];
while (i){
t=g2[i];
while (k){
j=g1[k];
if (ask[j].a>=t) break;
ask[j].x=num;
k=n1[k];
}
if (a[t]<=mid-1) num++;
i=n2[i];
}
while (k){
j=g1[k];
ask[j].x=num;
k=n1[k];
}
i=h3[p];
while (i){
t=g3[i];
if (ask[t*2].cnt+ask[t*2].x-ask[t*2-1].cnt-ask[t*2-1].x<=rank[t]-1){
ask[t*2].cnt+=ask[t*2].x;
ask[t*2-1].cnt+=ask[t*2-1].x;
ask[t*2].x=ask[t*2-1].x=0;
ask[t*2-1].p=ask[t*2].p=1;
add3(p*2+1,t);
}
else{
ask[t*2].x=ask[t*2-1].x=0;
ask[t*2-1].p=ask[t*2].p=0;
add3(p*2,t);
}
i=n3[i];
}
i=h1[p];
while (i){
t=g1[i];
if (ask[t].p) add1(p*2+1,t);else add1(p*2,t);
i=n1[i];
}
i=h2[p];
while (i){
t=g2[i];
if (a[t]<=mid-1) add2(p*2,t);else add2(p*2+1,t);
i=n2[i];
}
solve(p*2,l,mid-1);
solve(p*2+1,mid,r);
}
bool cmp(int a,int b){
return ask[a].a<ask[b].a;
}
int main(){
scanf("%d%d",&n,&m);
fo(i,1,n){
scanf("%d",&a[i]);
h[i]=a[i];
}
sort(h+1,h+n+1);
fo(i,1,n) a[i]=lower_bound(h+1,h+n+1,a[i])-h,add2(1,i);
fo(i,1,m){
scanf("%d%d%d",&j,&k,&l);
rank[i]=l;
add3(1,i);
ask[i*2-1].a=j-1;
ask[i*2].a=k;
}
fo(i,1,m*2) id[i]=i;
sort(id+1,id+m*2+1,cmp);
fo(i,1,m*2) add1(1,id[i]);
solve(1,1,n);
fo(i,1,m) printf("%d\n",ans[i]);
}