[HEOI2016/TJOI2016]排序 二分+线段树 || 80分的桶排序

原题:https://www.luogu.org/problemnew/show/P2824

题解:两种操作，区间升序，降序排序。先写一种暴力的桶排序，每次排序的时间复杂度为O(n)。原理和基数排序类似，开一个桶(数组)，扫一遍区间在对应的下标+1，然后在从前往后(后往前查）一遍就行了，只用扫两遍就行了，比快排还快。

80分的桶排序。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
const int N=1e5+10;
int a[N],cnt[N],n,q,m,p;

inline int rd(){
	int x=0;int f=1;char s=getchar();
	while(!isdigit(s)){if(s=='-')f=-1;s=getchar();}
	while(isdigit(s)) x=(x<<1)+(x<<3)+s-'0',s=getchar();
	return x*f;
}
inline void work(int opt,int l,int r){
	int lm=1e18,rm=0;
	for(int i=l;i<=r;i++){
		cnt[a[i]]++;lm=min(lm,a[i]);rm=max(rm,a[i]);
	}
	int hd=l;
	if(opt==0){
		for(int i=lm;i<=rm;i++)
			while(cnt[i]>0)cnt[i]--,a[hd++]=i;
	}else{
		for(int i=rm;i>=lm;i--)
			while(cnt[i]>0)cnt[i]--,a[hd++]=i; 
	}
}
int main(){
//	freopen("t.in","r",stdin);
	n=rd();q=rd();
	for(int i=1;i<=n;i++) a[i]=rd();
	for(int i=1,opt,x,y;i<=q;i++){
		opt=rd();x=rd();y=rd();
		work(opt,x,y);
	}
	p=rd();
	printf("%d\n",a[p]);
	return 0;
}

这道题的正解是要把问题转换成01排序，若这个序列只有0,1那么在处理区间时可以用线段树，只要查区间1的个数，区间修改就行了。考虑如何转化，二分第p个数的排名，将大于等于mid标记成1,小于mid的标记成0，01排序用线段树修改，查第p个位置是否为1就行了。单调性的简单证明，若第p个位置为1意味着这个数的排名应该是>=mid的，我们要找到最大的mid满足，改点为1。具有单调性。

细节:

1.区间修改是注意全部是0或1，避免查询时又l>r的情况。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#define ls now<<1
#define rs now<<1|1
using namespace std;
const int N=2e5+10;
int a[N],n,q,m,p,x[N],y[N],opt[N];
int sum[N<<2],tt[N<<2];
inline int rd(){
	int x=0;int f=1;char s=getchar();
	while(!isdigit(s)){if(s=='-')f=-1;s=getchar();}
	while(isdigit(s)) x=(x<<1)+(x<<3)+s-'0',s=getchar();
	return x*f;
}
inline void upd(int now){sum[now]=sum[ls]+sum[rs];}
inline void pushdown(int now,int l,int r){
	if(!tt[now]) return;
	tt[ls]=tt[rs]=tt[now];
	int mid=(l+r)>>1;
	if(tt[now]==1){//1为区间修改为1 ,-1为区间修改为0 
		int mid=(l+r)>>1;
		sum[ls]=(mid-l+1);sum[rs]=(r-mid);
	}else sum[ls]=sum[rs]=0;
	tt[now]=0;
}
void build(int now,int l,int r,int x){
	if(l==r){
		sum[now]=(a[l]>=x);tt[now]=0;
		return ;
	}
	int mid=(l+r)>>1;
	build(ls,l,mid,x);
	build(rs,mid+1,r,x);
	upd(now);tt[now]=0;
}
int query(int now,int l,int r,int x,int y){ 
	if(x<=l && r<=y){
		return sum[now];
	}
	int mid=(l+r)>>1;
	pushdown(now,l,r);
	int ans=0;
	if(x<=mid) ans+=query(ls,l,mid,x,y);
	if(mid<y)  ans+=query(rs,mid+1,r,x,y);
	return ans;
}
void change(int now,int l,int r,int x,int y,int c){
	if(x>y) return ;//一定要加 
	if(x<=l && r<=y){
		sum[now]=c*(r-l+1);
		tt[now]=(c?1:-1); 
		return ;
	}
	int mid=(l+r)>>1;
	pushdown(now,l,r);
	if(x<=mid) change(ls,l,mid,x,y,c);
	if(mid<y)  change(rs,mid+1,r,x,y,c);
	upd(now);
}
bool check(int val){
	build(1,1,n,val);
	for(int i=1;i<=q;i++){
		int num=query(1,1,n,x[i],y[i]);
		if(opt[i]==0){
			change(1,1,n,x[i],y[i]-num,0);
			change(1,1,n,y[i]-num+1,y[i],1); 
		}
		else {
			change(1,1,n,x[i],x[i]+num-1,1);
			change(1,1,n,x[i]+num,y[i],0);
		}
	}
	return query(1,1,n,p,p);;
}
int main(){
//	freopen("sort1.in","r",stdin);
	n=rd();q=rd();
	for(int i=1;i<=n;i++) a[i]=rd();
	for(int i=1;i<=q;i++){
		opt[i]=rd();x[i]=rd();y[i]=rd();
	}
	p=rd();
	int l=1;int r=n;
	while(l<=r){
		int mid=(l+r)>>1;
		if(check(mid))l=mid+1;
		else r=mid-1;
	}
	printf("%d\n",l-1);
	return 0;
}