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

题解：主席树模板。
区间内大于$\frac{r-l+1}{2}\,\,\,\,\,\,\,$ 的数只会有一个。

事实证明pool比newnode()快。。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
using namespace std;
const int Maxn=5e5+50;
inline int read()
{
    char ch=getchar();
    int i=0,f=1;
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){i=(i<<3)+(i<<1)+ch-'0';ch=getchar();}
    return i*f;
}
int n,m;

struct node
{
    node *lc,*rc;
    int sum;
    node():lc(NULL),rc(NULL),sum(0){}
}*tr[Maxn],POOL[10000000],*pool=POOL;
inline node* newnode()
{
    return ++pool;
}
inline void build(int l,int r,node *&x)
{
    if(l==r)return;
    x->lc=new node();
    x->rc=new node();
    int mid=(l+r)>>1;
    build(l,mid,x->lc);
    build(mid+1,r,x->rc);
}
inline void insert(int l,int r,node *x,node *&y,int v)
{
    y=newnode();
    y->lc=x->lc;
    y->rc=x->rc;
    y->sum=x->sum+1;
    if(l==r)return;
    int mid=(l+r)>>1;
    if(v<=mid)insert(l,mid,x->lc,y->lc,v);
    else insert(mid+1,r,x->rc,y->rc,v);
}
inline int query(int l,int r,node *x,node *y,int sum)
{
    if(l==r)return l;
    int mid=(l+r)>>1;
    if(y->lc->sum-x->lc->sum>sum)return query(l,mid,x->lc,y->lc,sum);
    else if(y->rc->sum-x->rc->sum>sum) query(mid+1,r,x->rc,y->rc,sum);
    else return 0;
}
int main()
{
    n=read(),m=read();
    tr[0]=newnode();
    build(1,n,tr[0]);
    for(int i=1;i<=n;i++)
    {
        int x=read();
        insert(1,n,tr[i-1],tr[i],x);
    }
    while(m--)
    {
        int l=read(),r=read();
        printf("%d\n",query(1,n,tr[l-1],tr[r],(r-l+1)>>1));
    }
}