bzoj 2223: [Coci 2009]PATULJCI

Description

Input

Output

10 3 1 2 1 2 1 2 3 2 3 3 8 1 2 1 3 1 4 1 5 2 5 2 6 6 9 7 10

Sample Input

no
yes 1
no
yes 1
no
yes 2
no
yes 3

Sample Output

HINT

Notice:输入第二个整数是序列中权值的范围Lim，即1<=ai(1<=i<=n)<=Lim。

1<=Lim<=10000

题目特意规定了询问的是出现次数超过一半以上的数，这意味着什么呢。

我们可以根据这点想到区间k大的查询。不同的是，我们这次直接找出现一半以上的数就可以了。

用主席树就可以搞定。

#include<cstdio>
#include<algorithm>
using namespace std;
struct tree
{
     int l,r;
     int ll,rr;
     int s;
}tr[8000001];
int tot;
struct number
{
     int x,t;
     bool operator< (number xt) const
     {
          return x<xt.x;
     }
}a[500001];
int f[500001],fx[500001];
int pl[500001];
inline void build(int l,int r)
{
     int p=tot;
     tr[p].l=l;
     tr[p].r=r;
     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 x)
{
	 if(p==0)
	      return 0;
     if(tr[p].l==l&&tr[p].r==r)
     {
          tot++;
          tr[tot].l=l;
          tr[tot].r=r;
          tr[tot].s=tr[p].s+x;
          return tot;
     }
     else
     {
          int mid=(tr[p].l+tr[p].r)/2;
          tot++;
          int pp=tot;
          tr[pp].l=tr[p].l;
          tr[pp].r=tr[p].r;
          if(l<=mid)
               tr[pp].ll=add(tr[p].ll,l,r,x);
          else
               tr[pp].ll=tr[p].ll;
          if(r>mid)
               tr[pp].rr=add(tr[p].rr,l,r,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)
{
     if(tr[p1].l==tr[p1].r)
     {
          if(tr[p2].s-tr[p1].s>k)
               return tr[p1].l;
          return 0;
     }
     else
     {
          int tmp=tr[tr[p2].ll].s-tr[tr[p1].ll].s;
          if(tmp>k)
               return find(tr[p1].ll,tr[p2].ll,k);
          return find(tr[p1].rr,tr[p2].rr,k);
     }
}
int main()
{
     int n,m,lim;
     scanf("%d%d",&n,&lim);
     int i;
     for(i=1;i<=n;i++)
     {
          scanf("%d",&a[i].x);
          a[i].t=i;
     }
     sort(a+1,a+1+n);
     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],fx[i],fx[i],1);
     }
     int s,t;
     scanf("%d",&m);
     for(i=1;i<=m;i++)
     {
          scanf("%d%d",&s,&t);
          int xt=f[find(pl[s-1],pl[t],(t-s+1)/2)];
          if(xt==0)
               printf("no\n");
          else
               printf("yes %d\n",xt);
     }
     return 0;
}