[Luogu P4168] [BZOJ 2724] [Violet]蒲公英

洛谷传送门

BZOJ传送门

题目背景

亲爱的哥哥：

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

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

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

哥哥你要快点回来哦！

                                           爱你的妹妹 Violet

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

“蒲公英吗……”

题目描述

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

为了简化起见，我们把所有的蒲公英看成一个长度为nnn的序列 (a1,a2..an)" role="presentation" style="position: relative;">(a1,a2..an)(a1,a2..an)(a_1,a_2..a_n) ，其中 aiaia_i 为一个正整数，表示第iii棵蒲公英的种类编号。

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

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

输入输出格式

输入格式：

第一行两个整数 n,mn,mn,m ，表示有nnn株蒲公英，m" role="presentation" style="position: relative;">mmm 次询问。

接下来一行nnn个空格分隔的整数 ai" role="presentation" style="position: relative;">aiaia_i ，表示蒲公英的种类

再接下来 mmm 行每行两个整数 l0,r0" role="presentation" style="position: relative;">l0,r0l0,r0l_0,r_0 ，我们令上次询问的结果为 xxx（如果这是第一次询问，则 x=0" role="presentation" style="position: relative;">x=0x=0x=0）。

令 l=(l0+x−1)modn+1,r=(r0+x−1)modn+1l=(l0+x−1)modn+1,r=(r0+x−1)modn+1l=(l_0+x-1) \mod n + 1,r=(r_0+x-1) \mod n + 1 ，如果 l>rl>rl>r，则交换 l,rl,rl,r 。

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

输出格式：

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

输入输出样例

输入样例#1：

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

输出样例#1：

1
2
1

说明

对于 20% 的数据，保证 1≤n,m≤3000" role="presentation" style="position: relative;">1≤n,m≤30001≤n,m≤30001\le n,m \le 3000 。

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

解题分析

这道题和作诗这道题很像，处理方式也大致相同，就不废话了…

不过aiaia_i的范围很大，需要先离散化一下。

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cctype>
#include <algorithm>
#include <cmath>
#define R register
#define IN inline
#define gc getchar()
#define W while
#define MX 50500
template <class T>
IN void in(T &x)
{x = 0; R char c = gc;for (; !isdigit(c); c = gc);for (;  isdigit(c); c = gc)x = (x << 1) + (x << 3) + c - 48;
}
int bel[MX], dat[MX], cpy[MX], sum[230][MX], ans[230][230], cnt[MX];
int dot, q, ctyp, lastans, dif, siz, bd;
IN int query(R int lef, R int rig)
{if(lef > rig) std::swap(lef, rig);int l = bel[lef], r = bel[rig];if(r - l <= 1){for (R int i = lef; i <= rig; ++i){++cnt[dat[i]];if(cnt[dat[i]] > cnt[lastans] || (cnt[dat[i]] == cnt[lastans] && dat[i] < lastans)) lastans = dat[i];}for (R int i = lef; i <= rig; ++i) --cnt[dat[i]];return lastans = cpy[lastans];}lastans = ans[l + 1][r - 1];int lb = l * siz, buf, tar;for (R int i = lef; i <= lb; ++i){++cnt[dat[i]];tar = cnt[lastans] + sum[r - 1][lastans] - sum[l][lastans];buf = cnt[dat[i]] + sum[r - 1][dat[i]] - sum[l][dat[i]];if(buf > tar || (buf == tar && dat[i] < lastans)) lastans = dat[i];}for (R int i = (r - 1) * siz + 1; i <= rig; ++i){++cnt[dat[i]];tar = cnt[lastans] + sum[r - 1][lastans] - sum[l][lastans];buf = cnt[dat[i]] + sum[r - 1][dat[i]] - sum[l][dat[i]];if(buf > tar || (buf == tar && dat[i] < lastans)) lastans = dat[i];}for (R int i = lef; i <= lb; ++i) --cnt[dat[i]];for (R int i = (r - 1) * siz + 1; i <= rig; ++i) --cnt[dat[i]];return lastans = cpy[lastans];
}
int main(void)
{int a, b, res;in(dot), in(q); siz = std::sqrt(dot) + 1;for (R int i = 1; i <= dot; ++i) in(dat[i]), cpy[i] = dat[i];std::sort(cpy + 1, cpy + 1 + dot);dif = std::unique(cpy + 1, cpy + 1 + dot) - cpy - 1;for (R int i = 1; i <= dot; ++i)//离散化{dat[i] = std::lower_bound(cpy + 1, cpy + 1 + dif, dat[i]) - cpy;bel[i] = (i - 1) / siz + 1;++sum[bel[i]][dat[i]];} bd = bel[dot];for (R int i = 1; i <= dif; ++i)for (R int j = 1; j <= bd; ++j) sum[j][i] += sum[j - 1][i];for (R int i = 1; i <= bd; ++i){res = 0;for (R int j = (i - 1) * siz + 1; j <= dot; ++j){++cnt[dat[j]];if(cnt[dat[j]] > cnt[res] || (cnt[dat[j]] == cnt[res] && dat[j] < res)) res = dat[j];ans[i][bel[j]] = res;}for (R int j = (i - 1) * siz + 1; j <= dot; ++j) --cnt[dat[j]];}W (q--){in(a), in(b);a = (a + lastans - 1) % dot + 1, b = (b + lastans - 1) % dot + 1;lastans = 0; printf("%d\n", query(a, b));}
}