分块

分块 · 浏览次数 : 7

小编点评

## 分块算法 **思路:** * 将序列分成几块,每个块大小为 \\(\\sqrt n\\)。 * 每块内,对每个元素进行暴力操作,并记录块内的和。 * 查询时,对每个块内所有元素的和进行计算。 * 分块的效率约为 \\(\\sqrt n\\),而线段树的效率约为 \\(\\log n\\)。 **代码:** ```c++ typedef long long ll;const int N = 100010, M = 350;int n, m, len;ll add[M], sum[M];int w[N];int get(int i) { //第i的位置映射到了哪一块 return i / len;}void change(int l, int r, int d) { if (get(l) == get(r)) { //段内,直接暴力 for (int i = l; i <= r; ++i) w[i] += d, sum[get(i)] += d; } else { int i = l, j = r; while (get(i) == get(l)) w[i] += d, sum[get(i)] += d, ++i; while (get(j) == get(r)) w[j] += d, sum[get(j)] += d, --j; for (int k = get(i); k <= get(j); ++k) { sum[k] += len * d; /* 这里有童鞋可能会问:如果最后一段有可能不足len长呢? 我们可以先将它补全,不影响最后的答案,因为题目是不会询问到数列以外的部分的 */ add[k] += d; } }}ll query(int l, int r) { ll res = 0; if (get(l) == get(r)) { //段内,直接暴力 for (int i = l; i <= r; ++i) res += w[i] + add[get(i)]; } else { int i = l, j = r; while (get(i) == get(l)) res += w[i] + add[get(i)], ++i; while (get(j) == get(r)) res += w[j] + add[get(j)], --j; for (int k = get(i); k <= get(j); ++k) res += sum[k]; } return res;}int main() { scanf(\"%d%d\", &n, &m); len = sqrt(n); for (int i = 1; i <= n; ++i) { scanf(\"%d\", &w[i]); sum[get(i)] += w[i]; } char op[2]; int l, r, d; while (m--) { scanf(\"%s%d%d\", op, &l, &r); if (*op == 'C') { scanf(\"%d\", &d); change(l, r, d); } else printf(\"%lld\\", query(l, r)); } return 0;} ``` **代码说明:** * `get(int i)` 函数用于计算元素在块中的位置。 * `change(int l, int r, int d)` 函数用于在块内对元素进行暴力操作。 * `query(int l, int r)` 函数用于计算块内所有元素的和。

正文

首先我们要了解一个问题:为什么要用分块
我们拿一道题目举例:

例题

给定一个长度为 \(N\) 的数列 \(A\),以及 \(M\) 条指令,每条指令可能是以下两种之一:

  1. C l r d,表示把 \(A[l],A[l+1],…,A[r]\) 都加上 \(d\)
  2. Q l r,表示询问数列中第 \(l \sim r\) 个数的和。

对于每个询问,输出一个整数表示答案。

输入格式

第一行两个整数 \(N,M\)

第二行 \(N\) 个整数 \(A[i]\)

接下来 \(M\) 行表示 \(M\) 条指令,每条指令的格式如题目描述所示。

输出格式

对于每个询问,输出一个整数表示答案。

每个答案占一行。

数据范围

\(1 \le N,M \le 10^5\),
\(|d| \le 10000\),
\(|A[i]| \le 10^9\)

输入样例:

10 5
1 2 3 4 5 6 7 8 9 10
Q 4 4
Q 1 10
Q 2 4
C 3 6 3
Q 2 4

输出样例:

4
55
9
15

大家可能会用线段树或树状数组解决这题,但是它们的码量和思路……
段错误树写一大堆递归(还要用懒标记),树状数组思路难想,理解很难。
然后就是调代码(都是血与泪的教训)
分块本质上就是一个优化的暴力算法,它的思路通俗易懂,虽然时间复杂度不及那些高级数据结构,但仔细想想,线段树和树状数组处理每一步的时的效率为\(\log n\),分块的效率则为\(\sqrt n\),显然\(\sqrt n\)\(\log n\)要处理的次数多上不少,但也是属于能过的状态,况且线段树因为巨大的常数有时还不如分块。

实测:

在这里插入图片描述

分块的思路

我们可以将一个序列分成\(\sqrt n\)块,如下图(当然最后一块不够的话就能分多少就分多少,不影响最终答案)
在这里插入图片描述
这里需要运用一点懒标记的思想,我们需要维护两个东西:
1、add:本段的所有数都要加上add
2、sum:本段的真实和是多少(算上add)
第一个操作:修改
在这里插入图片描述
修改的时候需要分两种情况:
1、完整段(蓝色部分):
\(add = add + d\)
\(sum = sum + d \times length\)(length为段的长度)
2、段内(红色部分):
直接暴力,枚举所有数
\(w_i = w_i + d\)
\(sum = sum + d\)
第二个操作:查询
同上,查询时也需要分两种情况:
1、完整段(蓝色部分):
直接累加sum即可
2、段内(红色部分):
直接暴力,枚举所有数,求和

其实如果在考试中出现这种有关序列的问题,如果实在想不出来正解,那就打个分块,一般能拿到70~80分

代码

#include <cmath>
#include <iostream>

using namespace std;

typedef long long ll;

const int N = 100010, M = 350;

int n, m, len;
ll add[M], sum[M];
int w[N];

int get(int i) { //	第i的位置映射到了哪一块
    return i / len;
}

void change(int l, int r, int d) {
    if (get(l) == get(r)) { //	段内,直接暴力
        for (int i = l; i <= r; ++i)
            w[i] += d, sum[get(i)] += d;
    } else {
        int i = l, j = r;
        while (get(i) == get(l))
            w[i] += d, sum[get(i)] += d, ++i;
        while (get(j) == get(r))
            w[j] += d, sum[get(j)] += d, --j;
        for (int k = get(i); k <= get(j); ++k) {
            sum[k] += len * d;
            /*
            这里有童鞋可能会问:如果最后一段有可能不足len长呢?
            我们可以先将它补全,不影响最后的答案,因为题目是不会询问到数列以外的部分的
            */
            add[k] += d;
        }
    }
}

ll query(int l, int r) {
    ll res = 0;
    if (get(l) == get(r)) { //	段内,直接暴力
        for (int i = l; i <= r; ++i)
            res += w[i] + add[get(i)];
    } else {
        int i = l, j = r;
        while (get(i) == get(l))
            res += w[i] + add[get(i)], ++i;
        while (get(j) == get(r))
            res += w[j] + add[get(j)], --j;
        for (int k = get(i); k <= get(j); ++k)
            res += sum[k];
    }

    return res;
}

int main() {
    scanf("%d%d", &n, &m);
    len = sqrt(n);
    for (int i = 1; i <= n; ++i) {
        scanf("%d", &w[i]);
        sum[get(i)] += w[i];
    }

    char op[2];
    int l, r, d;
    while (m--) {
        scanf("%s%d%d", op, &l, &r);
        if (*op == 'C') {
            scanf("%d", &d);
            change(l, r, d);
        } else
            printf("%lld\n", query(l, r));
    }

    return 0;
}

与分块相似的内容:

分块

首先我们要了解一个问题:为什么要用分块 我们拿一道题目举例: ### 例题 给定一个长度为 $N$ 的数列 $A$,以及 $M$ 条指令,每条指令可能是以下两种之一: 1. `C l r d`,表示把 $A[l],A[l+1],…,A[r]$ 都加上 $d$。 2. `Q l r`,表示询问数列中第

算法学习笔记(29):分块

分块 这是一种基于根号的算法,核心为大块标记,散块暴力,做到复杂度的平衡。 可能第一个想到于此相关的就是莫队吧,这是利用分块优化暴力的方法。 目录分块Rmq Problem / mex[国家集训队] 排队 - 洛谷[TJOI2009] 开关 - 洛谷[Violet] 蒲公英 - 洛谷小小总结 Rmq

14.3 Socket 字符串分块传输

首先为什么要实行分块传输字符串,一般而言`Socket`套接字最长发送的字节数为`8192`字节,如果发送的字节超出了此范围则后续部分会被自动截断,此时将字符串进行分块传输将显得格外重要,分块传输的关键在于封装实现一个字符串切割函数,将特定缓冲区内的字串动态切割成一个个小的子块,当切割结束后会得到该数据块的个数,此时通过套接字将个数发送至服务端此时服务端在依次循环接收数据包直到接收完所有数据包之后

fepk文件格式说明

1 卫星影像金字塔分块原理说明 通常我们在工作中使用的卫星影像数据,轻则几百M,重则几百个G甚至上TB级。影像数据太大,是大家经常会遇到的一个问题,尤其是想下载一个省以上数据的时候该问题尤为突出。那么该问题是否有一个比较好的解决方案呢? 以全球为例,我们以19级为例,共有2^18 * 2^17 张瓦

【学习笔记】基础算法:二次离线莫队/回滚莫队

【学习笔记】基础算法:二次离线莫队/回滚莫队 二次离线莫队 前置知识:莫队 前置知识:值域分块 值域分块,就是对 \(A\) 的值域进行分块,每个块维护该值域内数的个数 众所周知,莫队的复杂度是 \(O(n \sqrt m)\) 的,而在维护一些问题时左右端点移动一格并不是 \(\mathcal O

OI 数论中的上界估计与时间复杂度证明

渐进符号、约数函数、整除分块嵌套与杜教筛.

第129篇:JS模块化开发

好家伙,本篇为《JS高级程序设计》第二十六章“模块”学习笔记 JS开发会遇到代码量大和广泛使用第三方库的问题。 解决这个问题的方案通 常需要把代码拆分成很多部分,然后再通过某种方式将它们连接起来。 若代码量较大,我们使用模块化开发的模式,也能够使代码容易维护 我们需要模块 1.模块模式 把逻辑分块,

分页查询及其拓展应用案例

分页查询 分页查询是处理大量数据时常用的技术,通过分页可以将数据分成多个小部分,方便用户逐页查看。SQLAlchemy 提供了简单易用的方法来实现分页查询。 本篇我们也会在最终实现这样的分页效果: 1. 什么是分页查询 分页查询是将查询结果按照一定数量分成多页展示,每页显示固定数量的记录。分页查询通

分类模型的算法性能评价

一、概述 分类模型是机器学习中一种最常见的问题模型,在许多问题场景中有着广泛的运用,是模式识别问题中一种主要的实现手段。分类问题概况起来就是,对一堆高度抽象了的样本,由经验标定了每个样本所属的实际类别,由特定算法训练得到一个分类器,输入样本属性即自动计算出其所属类别,从而完成特定的识别任务。依实现原

半夜被慢查询告警吵醒,limit深度分页的坑

分享是最有效的学习方式。 博客:https://blog.ktdaddy.com/ 故事 梅雨季,闷热的夜,令人窒息,窗外一道道闪电划破漆黑的夜幕,小猫塞着耳机听着恐怖小说,辗转反侧,终于睡意来了,然而挨千刀的手机早不振晚不振,偏偏这个时候振动了一下,一个激灵,没有按捺住对内容的好奇,点开了短信,卧