浅谈暴力数据结构——分块
分块是一种思想,对于 n n n 个元素,分成 n \sqrt n n 个长度为 n \sqrt n n 的块,若非完全平方数,会有剩余的角块,分块实际上是一个只有三层的 n \sqrt n n 叉树,通过对完整块的整体操作和对不完整块的暴力操作而使复杂度尽可能的低。
通过适当的划分,预处理一部分信息保存下来,用空间换取时间,达到时空平衡。
1.块的划分
一般地,对于第 i i i 个元素,属于第 ⌊ i − 1 n ⌋ + 1 \lfloor\frac{i-1}{\sqrt n}\rfloor+1 ⌊ni−1⌋+1 个块,若块需要维护的信息一部分是原有的(如:和,极值等),那么进行预处理。
代码实现:
int siz=(int)sqrt(n);
for(int i=1;i<=n;i++)
belong[i]=(i-1)/siz+1;
其中 s i z siz siz 表示块的大小 , b e l o n g [ i ] belong[i] belong[i] 表示第 i i i 个元素属于哪个块。
2.例题
(Loj 6277 数列分块入门 1 )[https://loj.ac/p/6277]
-
对于添加中整块的部分在块内做标记,块外的部分直接暴力加。
-
查询 i i i 的结果就是 a r r [ i ] + a d d [ b e l o n g [ i ] ] arr[i]+add[belong[i]] arr[i]+add[belong[i]] , a d d [ i ] add[i] add[i]表示块内维护的加法标记。
解释: L L L 与 R R R 指块的左右边界
#include <cmath>
#include <cstdio>
#include <cstring>
#include <stdio.h>
using namespace std;
const int A = 1e5 + 11;
inline int read() {
char c = getchar(); int x = 0, f = 1;
for ( ; !isdigit(c); c = getchar()) if(c == '-') f = -1;
for ( ; isdigit(c); c = getchar()) x = x * 10 + (c ^ 48);
return x * f;
}
int n, m, t;
int a[A], add[A], L[A], R[A], pos[A];
void change(int l, int r, int c) {
int p = pos[l], q = pos[r];
if(p == q) {
for(int i = l; i <= r; i++) a[i] += c;
return;
}
for(int i = p + 1; i < q; i++) add[i] += c;
for(int i = l; i <= R[p]; i++) a[i] += c;
for(int i = L[q]; i <= r; i++) a[i] += c;
}
int ask(int now) {
int p = pos[now];
return a[now] + add[p];
}
int main() {
n = read();
for(int i = 1; i <= n; i++) a[i] = read();
t = sqrt(n);
for(int i = 1; i <= t; i++) {
L[i] = (i - 1) * sqrt(n) + 1;
R[i] = i * sqrt(n);
}
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;
}
}
m = n;
while(m--) {
int opt, l, r, c;
opt = read(), l = read(), r = read(), c = read();
if(opt == 0) change(l, r, c);
else printf("$d\n",ask(r));
}
}