整体二分--HDU5412 动态区间第k小

传送门

整体二分！带修改的和不带修改的其实差不多啦

详细解法在https://blog.csdn.net/sizeof_you/article/details/81738217

 

#include<iostream>
#include<cstdio>
#include<cstdio>
#include<cstring>
#include<cmath>
#define maxn 100005
using namespace std;
int n,m,ans[maxn],f[maxn],tot,a[maxn],now1,now2;

inline int rd(){
	int x=0,f=1;char c=' ';
	while(c<'0' || c>'9') {if(c=='-')f=-1;c=getchar();}
	while(c<='9' && c>='0') x=x*10+c-'0',c=getchar();
	return x*f;
}

struct Query{
	int l,r,k,cnt,num;
	Query(){}
	Query(int L,int R,int K,int CNT,int NUM){
		l=L,r=R,k=K,cnt=CNT,num=NUM;
	} 
}q[maxn*3],q1[maxn*3],q2[maxn*3];

inline void add(int x,int y){
	for(;x<=n;x+=x&-x)
		f[x]+=y;
}

inline int query(int x){
	int sum=0;
	for(;x;x-=x&-x)
		sum+=f[x];
	return sum;
}

inline void calc(int vl,int vr,int l,int r){
	for(int i=l;i<=r;i++){
		if(q[i].k)
			q[i].cnt=query(q[i].r)-query(q[i].l-1);//找到这个询问区间内有多少比mid小的数 
		else if(q[i].r<=vr)//如果修改的数小于mid说明对后面有影响要修改 
			add(q[i].l,q[i].cnt);
	}
	for(int i=l;i<=r;i++)
		if(!q[i].k && q[i].r<=vr) add(q[i].l,-q[i].cnt);//把改了的改回去 
	now1=now2=0;
	for(int i=l;i<=r;i++)
		if(q[i].k){//把询问和操作分成两部分 
			if(q[i].cnt>=q[i].k) q1[++now1]=q[i];
			else q[i].k-=q[i].cnt,q2[++now2]=q[i];
		}
		else{
			if(q[i].r<=vr) q1[++now1]=q[i];
			else q2[++now2]=q[i];
		}
	int now=l;
	for(int i=1;i<=now1;i++) q[now++]=q1[i];//将两部分重新合并 
	for(int i=1;i<=now2;i++) q[now++]=q2[i];
	return;
}

inline void solve(int l,int r,int vl,int vr){
	if(vl==vr){
		for(int i=l;i<=r;i++)
			if(q[i].k) ans[q[i].num]=vl;
		return;
	}
	int mid=(vl+vr)>>1;
	calc(vl,mid,l,r);//计算vl~mid的，将询问分开 
	int tmp=now1;
	if(tmp) solve(l,l+tmp-1,vl,mid);//不能直接写now1 
	if(l+tmp<=r) solve(l+tmp,r,mid+1,vr);//递归求左右两部分 
	return;
}

int main(){
	while(~scanf("%d",&n)){
		tot=now1=now2=0; memset(f,0,sizeof f); memset(ans,0,sizeof ans);
		for(int i=1;i<=n;i++){
			a[i]=rd(); q[++tot]=Query(i,a[i],0,1,0);
		}
		m=rd();
		for(int i=1;i<=m;i++){
			int ty=rd(),x,y,z;
			if(ty==1){
				x=rd(),y=rd();
				q[++tot]=Query(x,a[x],0,-1,0);//修改相当于删掉一个数再添上另一个 
				q[++tot]=Query(x,y,0,1,0);
				a[x]=y; ans[i]=0;
			}
			else {
				x=rd(),y=rd(),z=rd();
				q[++tot]=Query(x,y,z,0,i);
			}
		}
		solve(1,tot,1,1e9);
		for(int i=1;i<=m;i++)
			if(ans[i]>0) printf("%d\n",ans[i]);
	}
	return 0;
}