[CSP-S模拟测试]:大新闻(主席树)

题目传送门(内部题20)


输入格式

第一行为两个数$n,m$,意义如题所述。
接下来一行$n$个数,代表一开始$n$条大新闻的$naive$值。
接下来$m$行,每行一个操作,输入格式如下:
读入$1$,代表第一种事件。
读入$2,x$,代表第二种事件。
读入$3,l,r,k$,代表第三种事件。


输出格式

对于每一个第三种事件输出一行整数,代表答案


样例

样例输入:

6 8
2 7 4 3 5 9
3 2 5 3
1
2 4
3 1 4 2
2 6
3 1 7 5
1
3 3 6 4

样例输出:

5
4
6
9


数据范围与提示

样例解释:

初始序列为$2,7,4,3,5,9$。
询问$7,4,3,5$的第$3$小值,答案为$5$。
操作后序列变成$7,4,3,5,9$。
操作后序列变成$4,7,4,3,5,9$。
询问$4,7,4,3$的第$2$小值,答案为$4$。
操作后序列变成$6,4,7,4,3,5,9$。
询问序列$6,4,7,4,3,5,9$的第$5$小值,答案为$6$。
操作后序列变成$4,7,4,3,5,9$。
询问序列$4,3,5,9$的第$4$小值,答案为$9$。

数据范围:

对于$30\%$的数据,满足$1leqslant n,mleqslant 2,000$。
对于$50\%$的数据,满足$1leqslant n,mleqslant 30,000$。
对于$70\%$的数据,满足$1leqslant n,mleqslant 80,000$。
对于$100\%$的数据,满足$1leqslant n,mleqslant 2 imes 10^5$。
保证数据合法,当没有新闻时不存在删除操作。
保证$1leqslant naive$值$leqslant 10^9$,保证$1leqslant lleqslant rleqslant $当前新闻数,$kleqslant r−l+1$。


题解

观察数据范围和操作,发现是主席树板子题(然而并不会打……)。

不会的先学一下主席树吧。

正解直接用的树上差分,少个$log$。

时间复杂度:$(n+m)log^2n$。

期望得分:$100$分。

实际得分:$100$分。


代码时刻

#include<bits/stdc++.h>
using namespace std;
int n,m;
int a[200001],root[500000],tot;
int trsum[10000001],lson[10000001],rson[10000001];
void insert(int &x,int pre,int l,int r,int w)
{
	x=++tot;
	trsum[x]=trsum[pre]+1;
	lson[x]=lson[pre];
	rson[x]=rson[pre];
	if(l==r)return;
	int mid=(l+r)>>1;
	if(w<=mid)insert(lson[x],lson[pre],l,mid,w);
	else insert(rson[x],rson[pre],mid+1,r,w);
}
int ask(int x,int pre,int l,int r,int k)
{
	if(l==r)return l;
	int mid=(l+r)>>1;
	if(trsum[lson[x]]-trsum[lson[pre]]>=k)return ask(lson[x],lson[pre],l,mid,k);
	else return ask(rson[x],rson[pre],mid+1,r,k-trsum[lson[x]]+trsum[lson[pre]]);
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=n;i;i--)
		scanf("%d",&a[i]);
	for(int i=1;i<=n;i++)
		insert(root[i],root[i-1],1,1<<30,a[i]);
	while(m--)
	{
		int opt;
		scanf("%d",&opt);
		switch(opt)
		{
			case 1:n--;break;
			case 2:
				int x;
				scanf("%d",&x);
				n++;
				insert(root[n],root[n-1],1,1<<30,x);
				break;
			case 3:
				int l,r,k;
				scanf("%d%d%d",&l,&r,&k);
				printf("%d
",ask(root[n-l+1],root[n-r],1,1<<30,k));
				break;
		}
	}
	return 0;
}

rp++

原文地址:https://www.cnblogs.com/wzc521/p/11447424.html