虽然AC了但是时间惨不忍睹...不科学....怎么会那么慢呢...

无修改的区间众数..分块, 预处理出Mode[i][j]表示第i块到第j块的众数, sum[i][j]表示前i块j出现次数(前缀和,事实上我是写后缀和..因为下标从0开始..), cnt[i][j][k]表示第i块中的前j个数中,k出现次数。预处理O(N1.5), 询问每次O(N0.5), 总O((N+M)N0.5)

-------------------------------------------------------------------

#include<cstdio>

#include<cstring>
#include<algorithm>
#include<cmath>
#include<cctype>
 
using namespace std;
 
typedef pair<int, int> pii;
 
const int maxb = 209;
const int maxn = 40009;
 
struct HASH {
int h[maxn], n;
HASH() : n(0) {
}
void Add(int v) {
h[n++] = v;
}
void Work() {
sort(h, h + n);
n = unique(h, h + n) - h;
}
inline int Hash(int v) {
return lower_bound(h, h + n, v) - h;
}
inline int _Hash(int v) {
return h[v];
}
} H;
 
int read() {
char c = getchar();
int ret = 0;
for(; !isdigit(c); c = getchar()) c = getchar();
for(; isdigit(c); c = getchar()) ret = ret * 10 + c - '0';
return ret;
}
 
int seq[maxn];
int sum[maxb][maxn], cnt[maxb][maxb][maxb], c[maxn], Id[maxb][maxn];
int N, Q, B, n;
pii Mode[maxb][maxb]; // <numId, cnt>
 
void Init() {
N = read(); Q = read();
B = (int) sqrt(N);
n = N / B;
if(N % B) n++;
for(int i = 0; i < N; i++)
H.Add(seq[i] = read());
H.Work();
for(int i = 0; i < N; i++)
seq[i] = H.Hash(seq[i]);
memset(sum, 0, sizeof sum);
for(int i = 0; i < N; i++)
sum[i / B][seq[i]]++;
for(int i = n; i--; )
for(int j = 0; j < H.n; j++)
sum[i][j] += sum[i + 1][j];
for(int i = 0; i < n; i++) {
memset(c, 0, sizeof c);
int Max = 0, numId;
for(int j = i; j < n; j++) {
int p = j * B;
for(int k = 0; k < B; k++, p++) if(++c[seq[p]] > Max) {
numId = seq[p];
Max = c[seq[p]];
} else if(c[seq[p]] == Max) {
numId = min(numId, seq[p]);
}
Mode[i][j] = make_pair(numId, Max);
}
}
memset(cnt, 0, sizeof cnt);
memset(Id, -1, sizeof Id);
for(int i = 0; i < n; i++) {
Id[i][H.n] = 0;
int p = i * B;
for(int j = 0; j < B; j++, p++) {
int v = seq[p];
if(!~Id[i][v])
Id[i][v] = Id[i][H.n]++;
cnt[i][j][Id[i][v]]++;
}
for(int j = B; j--; )
for(int k = 0; k < Id[i][H.n]; k++)
cnt[i][j][k] += cnt[i][j + 1][k];
}
}
 
inline int getSum(int l, int r, int v) {
return sum[l][v] - sum[r + 1][v];
}
 
inline int getCnt(int Block, int l, int r, int numId) {
return cnt[Block][l][Id[Block][numId]] - cnt[Block][r + 1][Id[Block][numId]];
}
 
int solve(int l, int r) {
if(l > r) swap(l, r);
int lb = l / B, rb = r / B, ans, CNT = 0;
memset(c, 0, sizeof c);
if(rb - lb > 1) {
pii &t = Mode[lb + 1][rb - 1];
ans = t.first;
CNT = t.second;
for(int i = l; i / B == lb; i++) {
int v = seq[i], Cnt = getCnt(lb, l % B, B - 1, v) + getCnt(rb, 0, r % B, v);
Cnt += getSum(lb + 1, rb - 1, v);
if(Cnt > CNT)
ans = v, CNT = Cnt;
else if(Cnt == CNT)
ans = min(ans, v);
}
for(int i = rb * B; i <= r; i++) {
int v = seq[i], Cnt = getCnt(lb, l % B, B - 1, v) + getCnt(rb, 0, r % B, v);
Cnt += getSum(lb + 1, rb - 1, v);
if(Cnt > CNT)
ans = v, CNT = Cnt;
else if(Cnt == CNT)
ans = min(ans, v);
}
} else if(lb + 1 == rb) {
for(int i = l; i / B == lb; i++) {
int v = seq[i], Cnt = getCnt(lb, l % B, B - 1, v) + getCnt(rb, 0, r % B, v);
if(Cnt > CNT)
ans = v, CNT = Cnt;
else if(Cnt == CNT)
ans = min(ans, v);
}
for(int i = rb * B; i <= r; i++) {
int v = seq[i], Cnt = getCnt(lb, l % B, B - 1, v) + getCnt(rb, 0, r % B, v);
if(Cnt > CNT)
ans = v, CNT = Cnt;
else if(Cnt == CNT)
ans = min(ans, v);
}
} else if(lb == rb) {
for(; l <= r; l++) {
int v = seq[l], Cnt = getCnt(lb, l % B, r % B, v);
if(Cnt > CNT)
ans = v, CNT = Cnt;
else if(Cnt == CNT)
ans = min(ans, v);
}
}
return H._Hash(ans);
}
 
void Work() {
int ans = 0;
while(Q--) {
int l = read(), r = read();
printf("%d\n", ans = solve((l + ans - 1) % N, (r + ans - 1) % N));
}
}
 
int main() {
Init();
Work();
return 0;
}

-------------------------------------------------------------------

2724: [Violet 6]蒲公英

Time Limit: 40 Sec  Memory Limit: 512 MB
Submit: 1140  Solved: 373
[Submit][Status][Discuss]

Description

Input

修正一下

l = (l_0 + x - 1) mod n + 1, r = (r_0 + x - 1) mod n + 1

Output

Sample Input

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

Sample Output

1
2
1

HINT

修正下:

n <= 40000, m <= 50000

Source

BZOJ 2724: [Violet 6]蒲公英( 分块 )的更多相关文章

  1. BZOJ 2724: [Violet 6]蒲公英 [分块 区间众数]

    传送门 题面太美不忍不放 分块分块 这种题的一个特点是只有查询,通常需要预处理:加入修改的话需要暴力重构预处理 预处理$f[i][j]$为第i块到第j块的众数,显然$f[i][j]=max{f[i][ ...

  2. BZOJ 2724 [Violet 6]蒲公英(分块)

    题意 在线区间众数 思路 预处理出 f[i][j] 即从第 i 块到第 j 块的答案.对于每个询问,中间的整块直接用预处理出的,两端的 sqrtn 级别的数暴力做,用二分查找它们出现的次数.每次询问的 ...

  3. BZOJ 2724: [Violet 6]蒲公英

    2724: [Violet 6]蒲公英 Time Limit: 40 Sec  Memory Limit: 512 MBSubmit: 1633  Solved: 563[Submit][Status ...

  4. [BZOJ 2724] [Violet 6] 蒲公英 【分块】

    题目链接:BZOJ - 2724 题目分析 这道题和 BZOJ-2821 作诗 那道题几乎是一样的,就是直接分块,每块大小 sqrt(n) ,然后将数字按照数值为第一关键字,位置为第二关键字排序,方便 ...

  5. BZOJ.2724.[Violet 6]蒲公英(静态分块)

    题目链接 区间众数 强制在线 考虑什么样的数会成为众数 如果一个区间S1的众数为x,那么S1与新区间S2的并的众数只会是x或S2中的数 所以我们可以分块先预处理f[i][j]表示第i到第j块的众数 对 ...

  6. 【刷题】BZOJ 2724 [Violet 6]蒲公英

    Description Input 修正一下 l = (l_0 + x - 1) mod n + 1, r = (r_0 + x - 1) mod n + 1 Output Sample Input ...

  7. 【BZOJ 2724】 2724: [Violet 6]蒲公英 (区间众数不带修改版本)

    2724: [Violet 6]蒲公英 Time Limit: 40 Sec  Memory Limit: 512 MBSubmit: 1908  Solved: 678 Description In ...

  8. 【BZOJ】2724: [Violet 6]蒲公英

    2724: [Violet 6]蒲公英 Time Limit: 40 Sec  Memory Limit: 512 MBSubmit: 2900  Solved: 1031[Submit][Statu ...

  9. 【BZOJ2724】[Violet 6]蒲公英 分块+二分

    [BZOJ2724][Violet 6]蒲公英 Description Input 修正一下 l = (l_0 + x - 1) mod n + 1, r = (r_0 + x - 1) mod n ...

随机推荐

  1. 九度OnlineJudge之1021:统计字符

    题目描述:     统计一个给定字符串中指定的字符出现的次数. 输入:     测试输入包含若干测试用例,每个测试用例包含2行,第1行为一个长度不超过5的字符串,第2行为一个长度不超过80的字符串.注 ...

  2. Unix/Linux环境C编程入门教程(8) FreeBSD CCPP开发环境搭建

    1. FreeBSD是一种自由类Unix操作系统,是由经过BSD.386BSD和4.4BSD发展而来的类Unix的一个重要分支.FreeBSD拥有超过200名活跃开发者和上千名贡献者.FreeBSD被 ...

  3. 《Linux 设备驱动程序》读后感。 并发,竞态,死锁。

    1. 概念 并发:在操作系统中,是指一个时间段中有几个程序都处于已启动运行到运行完毕之间,且这几个程序都是在同一个处理机上运行,但任一个时刻点上只有一个程序在处理机上运行. 来源: 1. Linux ...

  4. oracle字符集

    oracle server端字符集

  5. javascript 中 keyup、keypress和keydown事件

    keyup.keypress和keydown事件都是有关于键盘的事件 1. keydown事件在键盘的键被按下的时候触发,keyup 事件在按键被释放的时候触发    keydown.keypress ...

  6. poj2186 Popular Cows --- 强连通

    给一个有向图,问有多少结点是其它全部结点都能够到达的. 等价于,在一个有向无环图上,找出度为0 的结点.假设出度为0的结点仅仅有一个,那么这个就是答案.假设大于1个.则答案是0. 这题有环.所以先缩点 ...

  7. 1.一步一步学c#(一):.NET体系结构(知识点)

    一.C#和.NET关系     1. 首先C#语言有两个很重要的方面:第一它是为了Microsoft的.NET Framework的结合而设计的,其次它是一种为问题而设计解决问题的方法的语言,它有很多 ...

  8. 内容高度小于窗口高度时版权div固定在底部

    <!doctype html><html><head><meta charset="utf-8"><title>文档内容 ...

  9. C#实现多态之一抽象

    1. 抽象类.抽象方法.抽象属性的特点 (1)      关键字:abstract (2)      抽象类只能是其他类的基类 (3)      抽象成员必须存在于抽象类中,但抽象类可以没有抽象成员, ...

  10. Unity5UGUI 官方教程学习笔记(四)UI Image

    Image Source image:源图片  需要显示的图片 Color:颜色  会与图片进行颜色的混合 Material:材质 Image Type:  Simple   精灵只会延伸到适合Rec ...