bzoj3524 [Poi2014]Couriers 主席树

题目链接：传送门

3524: [Poi2014]Couriers

Time Limit: 20 Sec   Memory Limit: 256 MB
Submit: 2482   Solved: 983
[ Submit][ Status][ Discuss]

Description

给一个长度为n的序列a。1≤a[i]≤n。
m组询问，每次询问一个区间[l,r]，是否存在一个数在[l,r]中出现的次数大于(r-l+1)/2。如果存在，输出这个数，否则输出0。

Input

第一行两个数n，m。
第二行n个数，a[i]。
接下来m行，每行两个数l,r，表示询问[l,r]这个区间。

Output

m行，每行对应一个答案。

Sample Input

7 5
1 1 3 2 3 4 3
1 3
1 4
3 7
1 7
6 6

Sample Output

1
0
3
0
4

HINT

【数据范围】

n,m≤500000

2016.7.9重设空间，但未重测！

Source

By Dzy

此题可以当做主席树模板题。

一个数如果在一个区间内出现次数超过一半，那么排序后一定会在中间的位置。所以此题等同于求区间内第k大值，k=（r-l+1）/2。

直接上主席树即可。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<set>
using namespace std;
void read(int &x)
{
	x=0;
	int f=1;char c=getchar();
	while(c<'0' || c>'9'){if(c=='-')	f=-1;c=getchar();}
	while(c>='0' && c<='9')	{x=x*10+c-'0';c=getchar();}
	x*=f;
}
int p[2000000+500000*20],inde;
int lson[2000000+500000*20],rson[2000000+500000*20];
int rot[500005];
int n,q;
int a[500005];
void b_t(int l,int r,int id)
{
	if(l==r)	return;
	int mid=(l+r)>>1;
	lson[id]=++inde;
	rson[id]=++inde;
	b_t(l,mid,lson[id]);
	b_t(mid+1,r,rson[id]);
}
void add(int& now,int lasid,int l,int r,int x)
{
	now=++inde;
	if(l==r)
	{
		p[now]=p[lasid]+1;
		return;
	}
	lson[now]=lson[lasid];
	rson[now]=rson[lasid];
	int mid=(l+r)>>1;
	if(x<=mid)	add(lson[now],lson[lasid],l,mid,x);
	else	add(rson[now],rson[lasid],mid+1,r,x);
	p[now]=p[lson[now]]+p[rson[now]];
}
int que(int lid,int rid,int l,int r,int t)
{
	if(l==r)
	{
		if(p[rid]-p[lid]>t)	return l;
		return 0;
	}
	int mid=(l+r)>>1;
	if(p[lson[rid]]-p[lson[lid]]>t)	return que(lson[lid],lson[rid],l,mid,t);
	if(p[rson[rid]]-p[rson[lid]]>t)	return que(rson[lid],rson[rid],mid+1,r,t);
	return 0;
}
int main()
{
	read(n);
	read(q);
	rot[0]=++inde;
	b_t(1,n,rot[0]);
	for(int i=1;i<=n;i++)
	{
		read(a[i]);
		add(rot[i],rot[i-1],1,n,a[i]);
	}
	int l,r;
	for(int i=0;i<q;i++)
	{
		read(l);
		read(r);
		printf("%d\n",que(rot[l-1],rot[r],1,n,(r-l+1)/2));
	}
	return 0;
}