BZOJ3524 [Poi2014]Couriers 【主席树】

题目

给一个长度为n的序列a。1≤a[i]≤n。

m组询问，每次询问一个区间[l,r]，是否存在一个数在[l,r]中出现的次数大于(r-l+1)/2。如果存在，输出这个数，否则输出0。

输入格式

第一行两个数n，m。

第二行n个数，a[i]。

接下来m行，每行两个数l,r，表示询问[l,r]这个区间。

输出格式

m行，每行对应一个答案。

输入样例

7 5

1 1 3 2 3 4 3

1 3

1 4

3 7

1 7

6 6

输出样例

1

0

3

0

4

提示

【数据范围】

n,m≤500000

2016.7.9重设空间，但未重测！

题解

主席树板题

我都不好意思拿出来写博客

#include<iostream>
#include<cstdio>
#include<algorithm>
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define lbt(x) (x & -x)
using namespace std;
const int maxn = 500005,maxm = 10000005,INF = 1000000000;
inline int read(){
	int out = 0,flag = 1; char c = getchar();
	while (c < 48 || c > 57) {if (c == '-') flag = -1; c = getchar();}
	while (c >= 48 && c <= 57) {out = (out << 3) + (out << 1) + c - '0'; c = getchar();}
	return out * flag;
}
int rt[maxn],siz = 0,n,m,sum[maxm],ls[maxm],rs[maxm];
void update(int& u,int pre,int l,int r,int pos){
	sum[u = ++siz] = sum[pre] + 1;
	if (l == r) return;
	ls[u] = ls[pre]; rs[u] = rs[pre];
	int mid = l + r >> 1;
	if (mid >= pos) update(ls[u],ls[pre],l,mid,pos);
	else update(rs[u],rs[pre],mid + 1,r,pos);
}
int query(int u,int v,int l,int r,int len){
	if (l == r) return sum[u] - sum[v] > len ? l : 0;
	int mid = l + r >> 1,t = sum[ls[u]] - sum[ls[v]];
	if (t > len) return query(ls[u],ls[v],l,mid,len);
	else return query(rs[u],rs[v],mid + 1,r,len);
}
int main(){
	n = read(); m = read();
	REP(i,n) update(rt[i],rt[i - 1],1,n,read());
	while (m--){
		int l = read(),r = read();
		printf("%d\n",query(rt[r],rt[l - 1],1,n,(r - l + 1) >> 1));
	}
	return 0;
}