[Violet]蒲公英（分块）

题目背景

亲爱的哥哥：

你在那个城市里面过得好吗？

我在家里面最近很开心呢。昨天晚上奶奶给我讲了那个叫「绝望」的大坏蛋的故事的说！它把人们的房子和田地搞坏，还有好多小朋友也被它杀掉了。我觉得把那么可怕的怪物召唤出来的那个坏蛋也很坏呢。不过奶奶说他是很难受的时候才做出这样的事的……

最近村子里长出了一大片一大片的蒲公英。一刮风，这些蒲公英就能飘到好远的地方了呢。我觉得要是它们能飘到那个城市里面，让哥哥看看就好了呢！

哥哥你要快点回来哦！

爱你的妹妹 $Violet$

$Azure$ 读完这封信之后微笑了一下。

“蒲公英吗……”

题目描述

在乡下的小路旁种着许多蒲公英，而我们的问题正是与这些蒲公英有关。

为了简化起见，我们把所有的蒲公英看成一个长度为n的序列 $a_1,a_2..a_n$，其中$a_i$为一个正整数，表示第$i$棵蒲公英的种类编号。

而每次询问一个区间 $[l,r]$，你需要回答区间里出现次数最多的是哪种蒲公英，如果有若干种蒲公英出现次数相同，则输出种类编号最小的那个。

注意，你的算法必须是在线的

输入格式

第一行两个整数 $n,m$ ，表示有$n$株蒲公英，$m$ 次询问。

接下来一行$n$个空格分隔的整数$a_i$表示，蒲公英的种类

再接下来$m$行每行两个整数$l_0,r_0$我们令上次询问的结果为 $x$（如果这是第一次询问， 则 $x=0$）。

令 $l=(l_0+x-1)\mathrm{m}\mathrm{o}\mathrm{d}n+1,r=(r_0+x-1)\mathrm{m}\mathrm{o}\mathrm{d}n+1$如果$l>r$，则交换$l,r$ 。

最终的询问区间为$[l,r]$。

输出格式

输出$m$行。每行一个整数，表示每次询问的结果。

输入输出样例

输入 #1

6 3 
1 2 3 2 1 2 
1 5 
3 6 
1 5

输出 #1

1 
2 
1

说明/提示

对于 20% 的数据，保证 $1\le n,m\le 3000$

对于 100% 的数据，保证 $1\le n\le 40000,1\le m\le 50000,1\le a_i\le 10^9$

解法一

• 为什么快读要比$scanf$慢啊啊啊啊
• $scanf$
• $read$
• 定义$c[i][j][k]$表示从第$i$个块到第$j$个块数字$k$出现的次数，$f[i][j][0]$表示区间数字最多的出现次数，$f[i][j][1]$表示区间众数。
• 对于每次询问直接按照分块的思想，对于$[l,L)$和$(r,R]$朴素扫描累加次数更新答案，询问后减少次数复原数组即可。

$code1$

#include <bits/stdc++.h> 
using namespace std; 
const int maxn = 4e4 + 100; 
const int maxm = 37; 

template <class T> 
inline void read(T &s) {
	s = 0; T w = 1, ch = getchar(); 
	while (!isdigit(ch)) { if (ch == '-') w = -1; ch = getchar(); }
	while (isdigit(ch)) { s = (s << 1) + (s << 3) + (ch ^ 48); ch = getchar(); }
	s *= w; 
} 

int n, m, tot, t; 
int a[maxn], b[maxn], L[maxn], R[maxn], pos[maxn], now[2]; 
int c[maxm][maxm][maxn], f[maxm][maxm][2]; 

inline void discretize() {
	memcpy(b, a, sizeof(b)); 
	sort(b + 1, b + n + 1); 
	tot = unique(b + 1, b + n + 1) - b - 1; 
	for (int i = 1; i <= n; ++i) 
		a[i] = lower_bound(b + 1, b + tot + 1, a[i]) - b; 	
}

inline void work(int x, int y, int num) {
	++c[x][y][num]; 
	if (c[x][y][num] > now[0] || (c[x][y][num] == now[0] && num < now[1])) {
		now[0] = c[x][y][num]; 
		now[1] = num; 
	}
}

int query(int l, int r) {
	int p = pos[l], q = pos[r]; 
	int x = 0, y = 0; 
	if (p + 1 <= q - 1) 
		x = p + 1, y = q -1; 
	memcpy(now, f[x][y], sizeof(now)); 
	if (p == q) {
		for (int i = l; i <= r; ++i) work(x, y, a[i]); 
		for (int i = l; i <= r; ++i) --c[x][y][a[i]]; 
	}
	else {
		for (int i = l; i <= R[p]; ++i) work(x, y, a[i]); 
		for (int i = L[q]; i <= r; ++i) work(x, y, a[i]); 
		for (int i = l; i <= R[p]; ++i) --c[x][y][a[i]]; 
		for (int i = L[q]; i <= r; ++i) --c[x][y][a[i]]; 
	}
	return b[now[1]]; 
}

int main() {
	// freopen("1.in", "r", stdin); 

	read(n), read(m); 
	for (int i = 1; i <= n; ++i) scanf("%d", &a[i]); 

	discretize(); 
	
	t = pow(double(n), (double)1 / 3); 
	int len = t ? n / t : n; 
	for (int i = 1; i <= t; ++i) {
		L[i] = (i - 1) * len + 1; 
		R[i] = i * len; 
	}
	if (R[t] < n) {
		++t; L[t] = R[t - 1] + 1; R[t] = n; 
	}
	for (int i = 1; i <= t; ++i) 
		for (int j = L[i]; j <= R[i]; ++j) 
			pos[j] = i; 
	memset(f, 0, sizeof(f)); 
	memset(c, 0, sizeof(c)); 
	for (int i = 1; i <= t; ++i) {
		for (int j = i; j <= t; ++j) {
			for (int k = L[i]; k <= R[j]; ++k) 
				++c[i][j][a[k]]; 
			for (int k = 1; k <= tot; ++k) {
				if (c[i][j][k] > f[i][j][0]) {
					f[i][j][0] = c[i][j][k]; 
					f[i][j][1] = k; 
				}
			}
		}
	}

	int x = 0; 
	while (m--) {
		int l, r; 
		scanf("%d %d", &l, &r); 
		l = (l + x - 1) % n + 1; 
		r = (r + x - 1) % n + 1; 
		if (l > r) swap(l, r); 
		x = query(l, r); 
		printf("%d\n", x); 
	}
	return 0; 
}