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;
}