BZOJ3524 POI2014 Couriers 题解&代码

题意:给出一个长度为n的序列a满足1≤a[i]≤n。
又有m组询问,每次对于一个区间[l,r]问是否存在一个数在[l,r]中出现的次数大于(r-l+1)/2。如果存在,输出这个数,否则输出0。
思路:嘛,很纯粹的主席树…只不过空间限定太严格了点,竟然卡空间,我卡了好几次才AC…
不过我现在特别纠结为啥我BZOJ2223过不了

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn = 500005;
int n,m,l,r,x,cnt,rt[maxn],son[2][maxn*20],sum[maxn*20];
void Insert(int rl,int rn,int l,int r,int v)
{
    sum[rn]=sum[rl]+1;
    if(l==r)return;
    int mid = (l+r)/2;
    if(mid>=v) son[1][rn]=son[1][rl],son[0][rn]=++cnt,Insert(son[0][rl],son[0][rn],l,mid,v);
    else son[0][rn]=son[0][rl],son[1][rn]=++cnt,Insert(son[1][rl],son[1][rn],mid+1,r,v);

}
int query(int rl,int rn,int l,int r,int k)
{
    if(l==r)return l;
    int mid = (l+r)/2;
    if(sum[son[0][rn]]-sum[son[0][rl]]>k) return query(son[0][rl],son[0][rn],l,mid,k);
    if(sum[son[1][rn]]-sum[son[1][rl]]>k) return query(son[1][rl],son[1][rn],mid+1,r,k);
    return 0;
}
int main(void)
{
    scanf("%d%d",&n,&m);
    for(int i = 1; i <= n; i++)rt[i]=++cnt;
    for(int i = 1; i <= n; i++)
    {
        scanf("%d",&x);
        Insert(rt[i-1],rt[i],1,n,x);
    }
    for(int i = 1; i <= m; i++)
    {
        scanf("%d%d",&l,&r);
        printf("%d\n",query(rt[l-1],rt[r],1,n,(r-l+1)/2));
    }
    return 0;
}

你可能感兴趣的:(poi)