题目背景
亲爱的哥哥:
你在那个城市里面过得好吗?
我在家里面最近很开心呢。昨天晚上奶奶给我讲了那个叫「绝望」的大坏蛋的故事的说!它把人们的房子和田地搞坏,还有好多小朋友也被它杀掉了。我觉得把那么可怕的怪物召唤出来的那个坏蛋也很坏呢。不过奶奶说他是很难受的时候才做出这样的事的……
最近村子里长出了一大片一大片的蒲公英。一刮风,这些蒲公英就能飘到好远的地方了呢。我觉得要是它们能飘到那个城市里面,让哥哥看看就好了呢!
哥哥你要快点回来哦!
爱你的妹妹 V i o l e t Violet Violet
A z u r e Azure Azure 读完这封信之后微笑了一下。
“蒲公英吗……”
题目描述
在乡下的小路旁种着许多蒲公英,而我们的问题正是与这些蒲公英有关。
为了简化起见,我们把所有的蒲公英看成一个长度为n的序列 a 1 , a 2 . . a n a_1,a_2..a_n a1,a2..an,其中 a i a_i ai为一个正整数,表示第 i i i棵蒲公英的种类编号。
而每次询问一个区间 [ l , r ] [l,r] [l,r],你需要回答区间里出现次数最多的是哪种蒲公英,如果有若干种蒲公英出现次数相同,则输出种类编号最小的那个。
注意,你的算法必须是在线的
输入格式
第一行两个整数 n , m n,m n,m ,表示有 n n n株蒲公英, m m m 次询问。
接下来一行 n n n个空格分隔的整数 a i a_i ai表示,蒲公英的种类
再接下来 m m m行每行两个整数 l 0 , r 0 l_0,r_0 l0,r0我们令上次询问的结果为 x x x(如果这是第一次询问, 则 x = 0 x=0 x=0)。
令 l = ( l 0 + x − 1 ) m o d n + 1 , r = ( r 0 + x − 1 ) m o d n + 1 l=(l_0+x-1)\bmod n + 1,r=(r_0+x-1) \bmod n + 1 l=(l0+x−1)modn+1,r=(r0+x−1)modn+1如果 l > r l>r l>r,则交换 l , r l,r l,r 。
最终的询问区间为 [ l , r ] [l,r] [l,r]。
输出格式
输出 m m m行。每行一个整数,表示每次询问的结果。
输入输出样例
输入 #1
6 3
1 2 3 2 1 2
1 5
3 6
1 5
输出 #1
1
2
1
说明/提示
对于 20% 的数据,保证 1 ≤ n , m ≤ 3000 1\le n,m \le 3000 1≤n,m≤3000
对于 100% 的数据,保证 1 ≤ n ≤ 40000 , 1 ≤ m ≤ 50000 , 1 ≤ a i ≤ 1 0 9 1\le n \le 40000,1\le m \le 50000,1\le a_i \le 10^9 1≤n≤40000,1≤m≤50000,1≤ai≤109
解法一
- 为什么快读要比 s c a n f scanf scanf慢啊啊啊啊
- s c a n f scanf scanf
- r e a d read read
- 定义 c [ i ] [ j ] [ k ] c[i][j][k] c[i][j][k]表示从第 i i i个块到第 j j j个块数字 k k k出现的次数, f [ i ] [ j ] [ 0 ] f[i][j][0] f[i][j][0]表示区间数字最多的出现次数, f [ i ] [ j ] [ 1 ] f[i][j][1] f[i][j][1]表示区间众数。
- 对于每次询问直接按照分块的思想,对于 [ l , L ) [l,L) [l,L)和 ( r , R ] (r,R] (r,R]朴素扫描累加次数更新答案,询问后减少次数复原数组即可。
c o d e 1 code1 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;
}