bzoj 3524: [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

和2223是双倍经验。

直接上主席树用区间k大做法做就行了。

可是我的主席树直接被卡内存了。连左右端点都不能记录。。= =b

#include<cstdio>#include<algorithm>using namespace std;struct tree{     int ll,rr;     int s;}tr[10000001];int tot;struct number{     int x,t;}a[500001];inline bool cmp(number x,number y){     if(x.x<y.x)          return true;     return false;}int f[500001],fx[500001];int pl[500001];inline void build(int l,int r){     int p=tot;     if(l!=r)     {          int mid=(l+r)/2;          tot++;          tr[p].ll=tot;          build(l,mid);          tot++;          tr[p].rr=tot;          build(mid+1,r);     }}inline int add(int p,int l,int r,int loc,int x){ if(p==0)      return 0;     if(l==r)     {          tot++;          tr[tot].s=tr[p].s+x;          return tot;     }     else     {          int mid=(l+r)/2;          tot++;          int pp=tot;          if(loc<=mid)               tr[pp].ll=add(tr[p].ll,l,mid,loc,x);          else               tr[pp].ll=tr[p].ll;          if(loc>mid)               tr[pp].rr=add(tr[p].rr,mid+1,r,loc,x);          else               tr[pp].rr=tr[p].rr;          tr[pp].s=tr[tr[pp].ll].s+tr[tr[pp].rr].s;          return pp;     }}inline int find(int p1,int p2,int k,int l,int r){     if(l==r)     {          if(tr[p2].s-tr[p1].s>k)               return l;          return 0;     }     else     {       int mid=(l+r)/2;          int tmp=tr[tr[p2].ll].s-tr[tr[p1].ll].s;          if(tmp>k)               return find(tr[p1].ll,tr[p2].ll,k,l,mid);          return find(tr[p1].rr,tr[p2].rr,k,mid+1,r);     }}int main(){     int n,m;     scanf("%d%d",&n,&m);     int i;     for(i=1;i<=n;i++)     {          scanf("%d",&a[i].x);          a[i].t=i;     }     sort(a+1,a+1+n,cmp);     int px=1;     f[px]=a[1].x;     fx[a[1].t]=px;     for(i=2;i<=n;i++)     {          if(a[i].x!=a[i-1].x)               px++;          f[px]=a[i].x;          fx[a[i].t]=px;     }     tot++;     pl[0]=tot;     build(1,px);     for(i=1;i<=n;i++)     {       pl[i]=tot+1;          add(pl[i-1],1,px,fx[i],1);     }     int s,t;     for(i=1;i<=m;i++)     {          scanf("%d%d",&s,&t);          printf("%d\n",f[find(pl[s-1],pl[t],(t-s+1)/2,1,px)]);     }     return 0;}