BZOJ_3585_mex && BZOJ_3339_Rmq Problem_莫隊+分塊

Description

　　有一個長度為n的數組{a1,a2,...,an}。m次詢問，每次詢問一個區間內最小沒有出現過的自然數。

Input

　　第一行n,m。
　　第二行為n個數。
　　從第三行開始，每行一個詢問l,r。

Output

　　一行一個數，表示每個詢問的答案。

Sample Input

Sample Output

HINT

數據規模和約定
　　對於100%的數據：
　　1<=n,m<=200000
　　0<=ai<=109
　　1<=l<=r<=n

　　對於30%的數據：
　　1<=n,m<=1000

我的做法比較$sb$ 也比較裸，只能處理離線不修改的。

詢問莫隊，把權值分塊，找到第一個不滿的塊，暴力查即可。

好像主席樹也能做，挖坑待填。

代碼(3585&&3339)：

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
using namespace std;
char nc() {
	static char buf[100000],*p1,*p2;
	return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int rd() {
	register int x=0;
	register char s=nc();
	while(s<‘0‘||s>‘9‘)s=nc();
	while(s>=‘0‘&&s<=‘9‘)x=(x<<3)+(x<<1)+s-‘0‘,s=nc();
	return x;
}
#define N 200050
int pos[N],L[N],R[N],size,block,n,c[N],ans[N],h[N],ansblo[N],m;
struct A {
	int l,r,id;
}q[N];
bool cmp(const A &x,const A &y) {
	if(pos[x.l]!=pos[y.l]) return pos[x.l]<pos[y.l];
	return x.r<y.r;
}
void del(int x) {
	h[x]--;
	if(h[x]==0) ansblo[pos[x]]--;
}
void add(int x) {
	h[x]++;
	if(h[x]==1) ansblo[pos[x]]++;
}
int query() {
	int i,j;
	for(i=1;i<=block;i++) {
		if(ansblo[i]<R[i]-L[i]+1) {
			for(j=L[i];j<=R[i];j++) {
				if(!h[j]) return j;
			}
		}
	}
	return n;
}
void solve() {
	int i,l=1,r=0;
	for(i=1;i<=m;i++) {
		while(l<q[i].l) del(c[l]),l++;
		while(r>q[i].r) del(c[r]),r--;
		while(l>q[i].l) l--,add(c[l]);
		while(r<q[i].r) r++,add(c[r]);
		ans[q[i].id]=query();
	}
}
int main() {
	n=rd(); m=rd();
	int i,j;
	size=sqrt(n); block=n/size;
	for(i=1;i<=block;i++) {
		L[i]=R[i-1]+1; R[i]=size*i;
		for(j=L[i];j<=R[i];j++) {
			c[j]=rd(); c[j]=min(c[j],n); pos[j]=i;
		}
	}
	if(R[block]!=n) {
		block++; L[block]=R[block-1]+1; R[block]=n;
		for(i=L[block];i<=n;i++) {
			c[i]=rd(); c[i]=min(c[i],n); pos[i]=block;
		}
	}
	for(i=1;i<=m;i++) {
		q[i].l=rd(); q[i].r=rd(); q[i].id=i;
	}
	L[1]=0; pos[0]=1;
	sort(q+1,q+m+1,cmp);
	solve();
	for(i=1;i<=m;i++) printf("%d\n",ans[i]);
}

BZOJ_3585_mex && BZOJ_3339_Rmq Problem_莫隊+分塊