BZOJ 3489 A simple rmq problem —— KD - tree

题目链接：点我啊╭(╯^╰)╮

题目大意：

     多次询问 $（l,r）$ 内只出现一次的最大值是多少

解题思路：

    由于强制在线，但给了 $40sec$
    记 $L[i]$ 为数字 $i$ 上次出现的位置，$R[i]$ 为下次出现的位置
    对于数字 $a[i]$ 是否满足条件，即满足 $(i,L[a[i]],R[a[i]])$ 与 $(l,r)$ 的限制
    由于限制关系只与最大最小值有关
    因此可以用三维的 $KD-tree$ 来维护
    改一下 $query$ 里的估价函数即可

#include<bits/stdc++.h>
#define rint register int
#define deb(x) cerr<<#x<<" = "<<(x)<<'\n';
using namespace std;
typedef long long ll;
const int inf = 0x3f3f3f3f;
const int maxn = 1e5 + 5;
int n, m, rt, tot, now, ans;
int v[maxn], L[maxn], R[maxn], pos[maxn];
struct node{
	int pla[3], val;
} a[maxn];
struct tree{
	int mx[3], mn[3], sz, ls, rs;
	int val, ans;
	node place;
} t[maxn];

inline int New(){
	return ++tot;
}
inline bool cmp(const node &A, const node &B){
	return A.pla[now] < B.pla[now];
}

inline void update(int p){
	for(int i=0; i<=2; i++){
		t[p].mx[i] = t[p].mn[i] = t[p].place.pla[i];
		if(t[p].ls) t[p].mx[i] = max(t[p].mx[i], t[t[p].ls].mx[i]), \
					t[p].mn[i] = min(t[p].mn[i], t[t[p].ls].mn[i]);
		if(t[p].rs) t[p].mx[i] = max(t[p].mx[i], t[t[p].rs].mx[i]), \
					t[p].mn[i] = min(t[p].mn[i], t[t[p].rs].mn[i]);
	}
	t[p].sz = t[t[p].ls].sz + t[t[p].rs].sz + 1;
	t[p].ans = max(t[p].ans, max(t[t[p].ls].ans, t[t[p].rs].ans));
}

inline int build(int l, int r, int opt){
	if(l > r) return 0;
	int x = New(), mid = l + r >> 1; now = opt;
	nth_element(a+l, a+mid, a+r+1, cmp); t[x].place = a[mid];
	t[x].val = t[x].ans = a[mid].val;
	t[x].ls = build(l, mid-1, (opt+1)%3);
	t[x].rs = build(mid+1, r, (opt+1)%3);
	update(x); return x;
}

inline bool ck(int l, int r, int p){
	return t[p].ans>ans && t[p].mx[0]>=l && t[p].mn[0]<=r \
			&& t[p].mn[1]<l && t[p].mx[2]>r;
}

inline void query(int l, int r, int p){
	if(!p || !ck(l, r, p)) return;
	if(t[p].place.pla[0]>=l && t[p].place.pla[0]<=r \
		&& t[p].place.pla[1]<l && t[p].place.pla[2]>r)
			ans = max(ans, t[p].val);
	query(l, r, t[p].ls), query(l, r, t[p].rs);
}

int main() {
	scanf("%d%d", &n, &m);;
	for(int i=1; i<=n; i++) 
		scanf("%d", v+i), L[i] = pos[v[i]], R[pos[v[i]]] = i, pos[v[i]] = i;
	for(int i=1; i<=n; i++){
		a[i].val = v[i];
		a[i].pla[0] = i;
		a[i].pla[1] = L[i];
		a[i].pla[2] = R[i] ? R[i] : n + 1;
	}
	rt = build(1, n, 0);
	for(int i=1, l, r; i<=m; i++){
		scanf("%d%d", &l, &r);
		l = (l + ans) % n + 1, r = (r + ans) % n + 1;
		if(l > r) swap(l, r);
		ans = 0; 
		query(l, r, rt);
		printf("%d\n", ans);
	}
}