题意:给定一个序列n
,有m
次查询,每次查询一个区间[l,r]
,求区间中每一种数在区间中第一次出现的位置的中位数,强制在线。
如果考虑离线做法,从大到小枚举 r
,用线段树对每个 l
维护[l,r]
中数的种类数,查询时在线段树上二分,那么支持在线只需要简单可持久化一下,复杂度是 O((n+q)logn)
。
空间需要开到 2 * nlogn ,因为每次修改需要-1, 再 + 1,所以一共修改了两次,所以还有一个常数2。
代码:
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <iostream>
#include <string>
#include <cmath>
#include <vector>
#include <set>
#include <map>
#include <bitset>
#include <stack>
using namespace std;
#define REP(i,n) for ( int i=1; i<=int(n); i++ )
#define MP make_pair
#define PB push_back
#define SZ(x) (int((x).size()))
#define ALL(x) (x).begin(), (x).end()
#define X first
#define Y second
template<typename T> inline bool chkmin(T &a, const T &b) { return a > b ? a = b, 1 : 0; }
template<typename T> inline bool chkmax(T &a, const T &b) { return a < b ? a = b, 1 : 0; }
typedef long long LL;
typedef long double LD;
const int INF = 0x3f3f3f3f;
template <class T>
inline bool RD(T &ret) {
char c; int sgn;
if (c = getchar(), c == EOF) return 0;
while (c != '-' && (c<'0' || c>'9')) c = getchar();
sgn = (c == '-') ? -1 : 1 , ret = (c == '-') ? 0 : (c - '0');
while (c = getchar(), c >= '0'&&c <= '9') ret = ret * 10 + (c - '0');
ret *= sgn;
return 1;
}
template <class T>
inline void PT(T x) {
if (x < 0) putchar('-') ,x = -x;
if (x > 9) PT(x / 10);
putchar(x % 10 + '0');
}
typedef pair<LL, LL> pii;
const int N = 2e5 + 10;
const int M = 30 * 2 * N;
int ls[M], rs[M], root[N], tot, data[M];
int new_node(int lst = 0) {
ls[++ tot] = ls[lst];
rs[tot] = rs[lst];
data[tot] = data[lst];
return tot;
}
void build(int l, int r, int &rt) {
rt = new_node();
if(l == r) return ;
int m = (l + r) >> 1;
build(l, m, ls[rt]);
build(m + 1, r, rs[rt]);
}
void update(int pos, int val, int lst, int l, int r, int &rt) {
rt = new_node(lst);
data[rt] += val;
if(l == r) return ;
int m = (l + r) >> 1;
if(pos <= m) update(pos, val, ls[lst], l, m, ls[rt]);
else update(pos, val, rs[lst], m + 1, r, rs[rt]);
}
int query(int L, int R, int l, int r, int rt) {
if(L <= l && r <= R) return data[rt];
int m = (l + r) >> 1, ans = 0;
if(L <= m) ans += query(L, R, l, m, ls[rt]);
if(R > m) ans += query(L, R, m + 1, r, rs[rt]);
return ans;
}
int ask_kth(int k, int l, int r, int rt) {
if(l == r) return l;
int m = (l + r) >> 1;
return data[ls[rt]] >= k ? ask_kth(k, l, m, ls[rt]) : ask_kth(k - data[ls[rt]], m + 1, r, rs[rt]);
}
int a[N], pre[N];
int main() {
int T;
cin >> T;
for(int cas = 1; cas <= T; cas ++) {
int n, m;
scanf("%d %d", &n, &m);
memset(pre, -1, sizeof(pre[0]) * (n + 5));
for(int i = 1; i <= n; i ++) scanf("%d", &a[i]);
tot = 0;
build(1, n, root[n + 1]);
for(int i = n; i >= 1; i --) {
if(~ pre[a[i]]) update(pre[a[i]], -1, root[i + 1], 1, n, root[i]);
update(i, 1, root[~pre[a[i]] ? i : i + 1], 1, n, root[i]);
pre[a[i]] = i;
}
printf("Case #%d:", cas);
int ans = 0;
while(m --) {
int l, r;
scanf("%d %d", &l, &r);
l = (l + ans) % n + 1;
r = (r + ans) % n + 1;
if(l > r) swap(l, r);
int cnt = query(l, r, 1, n, root[l]);
ans = ask_kth((cnt + 1) >> 1, 1, n, root[l]);
printf(" %d", ans);
}
puts("");
}
}