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);elseadd(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;}