JZOJ5933. 【NOIP2018模拟10.27】百鸽笼

题目

在JOU管理员群里一共有n（n≤200000）个管理员，为了容纳这些管理员，vfk准备了n个鸽笼。每个 鸽笼中都装有一个咕咕能力值为vi的鸽子，每只鸽子能力值不一定相同。每当RU开始或结束或 咕咕咕时，管理员们就会对这些鸽笼进行操作。操作包括三种：
1 vfk和管理员们取出最左端的鸽笼。
2 v vfk和管理员们在最左侧新放入一个装有咕咕能力值为v的鸽子的鸽笼。
3 l r k vfk想要知道从左到右第l个到第r个鸽笼中鸽子的第k小的咕咕能力值。 小∞ ⇔想要知道每个3操作（即询问）的答案，请你来帮帮他。

题解

求区间第k小，很自然的想到主席树，
可是这里面有修改操作，
带修主席树？
显然这个并不是联赛考点，
注意观察可以知道，它只是修改第一个位置，
如果将a反过来，那么修改就是最后一个位置，
这样就可以解决了修改的问题。
于是就变成了一个简单的主席树问题。

code

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#define N 200003
#define P putchar
#define G getchar
using namespace std;
char ch;
void read(int &n)
{
	n=0;
	ch=G();
	while((ch<'0' || ch>'9') && ch!='-')ch=G();
	int w=1;
	if(ch=='-')w=-1,ch=G();
	while('0'<=ch && ch<='9')n=(n<<3)+(n<<1)+ch-'0',ch=G();
	n*=w;
}

void write(int x){if(x>9) write(x/10);P(x%10+'0');}

struct node
{
	int l,r,s;
}t[N*100];

int n,m,a[N],all,tot,root[N*2],v[N*2],k,V;
int cz[4][N];

void build(int pre,int x,int l,int r)
{
	if(l==r)
	{
		t[x].s=t[pre].s+1;
		return;
	}
	int m=(l+r)>>1;
	if(V<=m)
	{
		t[x].l=++tot;t[x].r=t[pre].r;
		build(t[pre].l,t[x].l,l,m);
	}
	else
	{
		t[x].l=t[pre].l;t[x].r=++tot;
		build(t[pre].r,t[x].r,m+1,r);
	}
	t[x].s=t[t[x].l].s+t[t[x].r].s;
}

void find(int x,int y,int l,int r)
{
	if(l==r)
	{
		V=l;
		return;
	}
	int v=t[t[y].l].s-t[t[x].l].s,m=(l+r)>>1;
	if(k<=v)find(t[x].l,t[y].l,l,m);
		else k=k-v,find(t[x].r,t[y].r,m+1,r);
}

int main()
{
	freopen("pigeon.in","r",stdin);
	freopen("pigeon.out","w",stdout);
	
	read(n);read(m);all=n;
	for(int i=1;i<=n;i++)read(a[i]),v[i]=a[i];
	for(int i=1;i<=m;i++)
	{
		read(cz[0][i]);
		if(cz[0][i]==2)read(cz[1][i]),v[++all]=cz[1][i];else
		if(cz[0][i]==3)read(cz[1][i]),read(cz[2][i]),read(cz[3][i]);
	}
	sort(v+1,v+1+all);
	all=unique(v+1,v+1+all)-v-1;
	for(int i=1;i<=n;i++)
		root[i]=++tot,V=lower_bound(v+1,v+1+all,a[n-i+1])-v,
			build(root[i-1],root[i],1,all);
	
	for(int i=1;i<=m;i++)
	{
		if(cz[0][i]==1)n--;else
		if(cz[0][i]==2)root[++n]=++tot,V=lower_bound(v+1,v+1+all,cz[1][i])-v,
			build(root[n-1],root[n],1,all);
		else
		{
			k=cz[3][i];
			find(root[n-cz[2][i]],root[n-cz[1][i]+1],1,all);
			write(v[V]);P('\n');
		}
	}
	
	return 0;
}