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[maxn],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;
}