[bzoj3524/2223][Poi2014]Couriers_主席树

Couriers bzoj-3524 Poi-2014

题目大意：给定n个数的序列，询问区间内是否存在一个在区间内至少出现了（区间长度>>1）次的数。如果有，输出该数，反之输出0。

注释：$1\le n,m\le 5\cdot 10^5$。

---

想法：主席树裸题。

主流做法就是弄一个Existence数组询问有没有这样的数，然后查询区间中位数即可。

但是可以在query的时候强行查询，因为没有输出0，直接输出即可。

最后，附上丑陋的代码... ...

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#define N 500010
using namespace std;
int cnt[N<<5],ls[N<<5],rs[N<<5],tot;
int a[N],b[N],root[N];
int build(int l,int r)
{
	int pos=++tot,mid=(l+r)>>1;
	if(l==r) return pos;
	ls[pos]=build(l,mid); rs[pos]=build(mid+1,r);
	return pos;
}
int update(int pre,int l,int r,int k)
{
	int pos=++tot,mid=(l+r)>>1;
	cnt[pos]=cnt[pre]+1;
	ls[pos]=ls[pre]; rs[pos]=rs[pre];
	if(l==r) return pos;
	if(k<=mid) ls[pos]=update(ls[pre],l,mid,k);
	else rs[pos]=update(rs[pre],mid+1,r,k);
	return pos;
}
// bool Existence(int x,int y,int l,int r,int Want)
// {
// 	if(l==r) return (cnt[y]-cnt[x])>=Want;
// 	int mid=(l+r)>>1;
// 	int LR=1;
// 	if(cnt[ls[y]]-cnt[ls[x]]>cnt[rs[y]]-cnt[rs[x]]) LR=0;
// 	else if(cnt[ls[y]]-cnt[ls[x]]<cnt[rs[y]]-cnt[rs[x]]) LR=2;
// 	if(LR==0) return Existence(ls[x],ls[y],l,mid,Want);
// 	if(LR==2) return Existence(rs[x],rs[y],mid+1,r,Want);
// 	return Existence(ls[x],ls[y],l,mid,Want)|Existence(rs[x],rs[y],mid+1,r,Want);
// }
// int query(int x,int y,int l,int r,int k)
// {
// 	int mid=(l+r)>>1;
// 	if(l==r) return b[l];
// 	int dlt=cnt[ls[y]]-cnt[ls[x]];
// 	if(dlt>=k) return query(ls[x],ls[y],l,mid,k);
// 	else return query(rs[x],rs[y],mid+1,r,k-dlt);
// }
int query(int x,int y,int l,int r,int k)
{
	int mid=(l+r)>>1;
	if(l==r) return l;
	int L=cnt[ls[y]]-cnt[ls[x]],R=cnt[rs[y]]-cnt[rs[x]];
	if(L>k) return query(ls[x],ls[y],l,mid,k);
	if(R>k) return query(rs[x],rs[y],mid+1,r,k);
	else return 0;
}
int main()
{
	int n,m; cin >> n >> m ;
	for(int i=1;i<=n;i++) scanf("%d",&a[i])/* ,b[i]=a[i] */;
	sort(b+1,b+n+1); /* int len=unique(b+1,b+n+1)-b-1; */
	root[0]=build(1,n);
	for(int i=1;i<=n;i++)
	{
		// int t=lower_bound(b+1,b+n+1,a[i])-b;
		root[i]=update(root[i-1],1,n,a[i]);
	}
	for(int x,y,i=1;i<=m;i++)
	{
		scanf("%d%d",&x,&y);
		// if(Existence(root[x-1],root[y],1,n,(y-x+1)/2+1)) printf("%d\n",query(root[x-1],root[y],1,n,(y-x+1)/2+1));
		// else printf("0\n");
		printf("%d\n",query(root[x-1],root[y],1,n,(y-x+1)>>1));
	}
	return 0;
}

小结：主席树真强... ...