BZOJ3524 [Poi2014]Couriers

标签：主席树

题目

题目传送门

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

分析

写完树剖和LCT后觉得主席树代码量还是很短很美滋滋的

区间中出现次数超过(r−l+1)/2的数也只可能有一个

建n棵权值线段树，然后查询的时候类似前缀和相减就可以了

第一次写主席树，对着hzw的代码调了好长时间

抱怨网上也没有一个正常点入门的主席树资料

code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define dep(i,a,b) for(int i=a;i>=b;i--)
#define ll long long
#define mem(x,num) memset(x,num,sizeof x)
using namespace std;
inline int read()
{
int f=1,x=0;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
const int maxn=1e7+6;
int root[maxn],ls[maxn],rs[m
d784
axn],sum[maxn];
int n,m,sz;
void update(int l,int r,int x,int &y,int v){
y=++sz;
sum[y]=sum[x]+1;
if(l==r)return;
ls[y]=ls[x];rs[y]=rs[x];int mid=(l+r)>>1;
if(v<=mid)update(l,mid,ls[x],ls[y],v);else update(mid+1,r,rs[x],rs[y],v);
}

int query(int L,int R){
int l=1,r=n,mid,x,y,temp=(R-L+1)/2;
x=root[L-1];y=root[R];
while(l<r){
if(sum[y]-sum[x]<=temp)return 0;
mid=(l+r)>>1;
if(sum[ls[y]]-sum[ls[x]]>temp){r=mid;x=ls[x];y=ls[y];}
else if(sum[rs[y]]-sum[rs[x]]>temp){l=mid+1;x=rs[x];y=rs[y];}
else return 0;
}
return l;
}
int main()
{
n=read(),m=read();
rep(i,1,n){
int x=read();
update(1,n,root[i-1],root[i],x);
}
rep(i,1,m){
int l=read(),r=read();
printf("%d\n",query(l,r));
}
return 0;
}