P4137 Rmq Problem / mex【权值线段树】【莫队】

传送门
给定长度为$N$的数组，求$l\sim r$ 的$mex$(从0开始第一个没有在区间出现的数)

分析

这种位置会覆盖前一个位置影响的，都考虑左边最近一个的位置

静态区间在线查询 $MEX$

主席树，维护最近值域出现的位置

在 $r$ 树上二分找值域的出现位置小于 $l$ 的最小值

证明：

假如一个数$x$ 在区间 $l\sim r$ 出现过，那么它在$r$数上对应值域的值，一定是大于等于$l$的

推得，一个数在$l\sim r$区间出现过的条件是，在$r$树上值域的值一定大于$l$

反之，则不出现。$mex$则是不出现的最小数，对应在树上就是值域最小且出现时间小于$l$

主席树

主席树第$r$个版本维护，区间$1\sim r$中，所有数出现的最后一个位置
考虑，建立一颗权值线段树，维护的是某个数最后一次的下标
分析查询，这样的话我们要找到$mex$，意味着要找，出现下标小于$l$且值最小的数，权值线段树找值最小的，肯定是在左边，什么时候找右边呢？当左边出现的所有数的下标都是大于$l$的，也就是最小值都大于$l$，说明这些数都在区间出现过，分析到这里就能写了

莫队

简单莫队，不带修改
同样分块$\sqrt{N}$，不同块$l$从小到大，相同块$r$从小到大奇偶排序
增加某个数，如果此数没出现过，查看当前数是否就是不增加时的$mex$，如果是，$while$一直找到$mex$，如果不是增加出现次数即可
删除某个数，如果次数出现过，但是删除后没了，跟$mex$比较大小，如果$mex$更小就跳过，如果$mex$更大，跟新$mex$为当前值，如果删除后还有，减少出现次数即可
这样就能维护$mex$了

• 注意，这里增加数的时候，会出现一个$while$循环，这里可能会导致复杂度变高，可以使用回滚莫队（不增莫队）解决

代码

主席树

//P4137 
/*
  @Author: YooQ
*/
#include <bits/stdc++.h>
using namespace std;
#define sc scanf
#define pr printf
#define ll long long
#define FILE_OUT freopen("out", "w", stdout);
#define FILE_IN freopen("P4137_1.in", "r", stdin);
#define debug(x) cout << #x << ": " << x << "\n";
#define AC 0
#define WA 1
#define INF 0x3f3f3f3f
const ll MAX_N = 1e6+5;
const ll MOD = 1e9+7;
int N, M, K;

int arr[MAX_N];
struct Tr {
	int k, l, r;
}tr[MAX_N<<2];
int indx = 0;
int root[MAX_N];

int clone(int rt) {
	tr[++indx] = tr[rt];
	return indx;
}

#define mid l + ((r-l)>>1)
void update(int& rt, int l, int r, int x, int k) {
	rt = clone(rt);
	if (l == r) {
		tr[rt].k = k;
		return;
	}
	if (x <= mid) update(tr[rt].l, l, mid, x, k);
	if (x  > mid) update(tr[rt].r, mid+1, r, x, k);
	tr[rt].k = min(tr[tr[rt].l].k, tr[tr[rt].r].k);
}

int query(int rt, int l, int r, int x) {
	if (l == r || !rt) return l;
	if (tr[tr[rt].l].k >= x) return query(tr[rt].r, mid+1, r, x);
	return query(tr[rt].l, l, mid, x);
}

void solve(){
	sc("%d%d", &N, &M);
	for (int i = 1; i <= N; ++i) {
		sc("%d", &arr[i]);
		++arr[i];
		root[i] = root[i-1];
		update(root[i], 1, 200001, arr[i], i);
	}
	
	int l, r;
	for (int i = 1; i <= M; ++i) {
		sc("%d%d", &l, &r);
		pr("%d\n", query(root[r], 1, 200001, l)-1);
	}
}

signed main()
{
	#ifndef ONLINE_JUDGE
//	FILE_IN
	FILE_OUT
	#endif
	int T = 1;//cin >> T;
	while (T--) solve();

	return AC;
}

莫队

pending...