Intro:

作为查询界的 \(O(1)\) 王者——前缀和的亲兄弟,差分,他可是修改界的 \(O(1)\) 王者


Prerequisite knowledge:

前缀和


Function: 仅单次询问的区间修改

模板题:洛谷P2367 语文成绩

先想一想朴素算法怎么做吧

对于输入的每一组 \((x,y)\) ,遍历序列 \(a_{x..y}\) ,每一项加上 \(z\) ,代码如下

while(p--)for(int i(x);i<=y;++i)a[i]+=z;

Time complexity: \(O(np)\)

Memory complexity: \(O(n)\)

这个时间复杂度明显不行,只要修改长度大了,修改时间就会变长

发现 查询只有一次

也就是说如果可以把数组变一下,使得修改 \(O(1)\) ,而查询 \(O(n)\) ,就可以接受了

我们构造一个数组 \(d_{1..n}\) ,其中

\(d_i=a_i-a_{i-1}, 1\leqslant i\leqslant n\)

也就是说 \(d\) 的每一项都是原来 \(a\) 的同位项与 \(a\) 的前一项的差

反过来看的话,\(a_i=a_{i-1}+d_i\)(是不是有点眼熟呢)

没错! \(a_i=\sum_{j=1}^id_j\) ,是前缀和!

考虑一下对 \(d_i\) 加上 \(v\) 对 \(a\) 会有什么效果

对于 \(a_{1..i-1}\) ,没效果

对于 \(a_{i..n}\) ,每一个值都加上了 \(v\) (认真思考下为什么)

所以对于每一次区间修改 \((l,r,v)\) (表示对 \(a_{l..r}\) 每一个元素加上 \(v\) ),对 \(d\) 来说,只不过是把 \(d_l\) 加上 \(v\) , \(d_{r+1}\) 减去 \(v\) 而已

这就是\(d_{1..n}\),差分数组

最后如果想要把\(d\)数组变回\(a\)数组,只要对\(d\)构造前缀和数组即可(这就是为什么我说差分是前缀和的逆运算了)


Code:

构造 \(d_{1..n}\) :
for(int i(n);i>=1;--i)d[i]=a[i]-a[i-1];

P.s: 其实顺序没有关系,但是如果原址构造必须倒序

修改
d[l]+=v,d[r+1]-=v;
更新原数组
for(int i(1);i<=n;++i)a[i]=a[i-1]+d[i]

Time complexity: \(O(n)\) 预处理 \(O(1)\) 修改 \(O(n)\) 查询(以后这三项简称为 \(O(n)-O(1)-O(n)\))

Memory complexity: \(O(n)\)

P.s 同样的,可以直接将原数组变成差分数组,则不需要额外空间,代码如下(必须倒序)

for(int i(n);i>=2;--i)a[i]-=a[i-1];

于是本题的答案就出来啦,具体见代码(记得开long long)

//This program is written by Brian Peng.
#pragma GCC optimize("Ofast","inline","-ffast-math")
#pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include<bits/stdc++.h>
using namespace std;
#define int long long
#define Rd(a) (a=read())
#define Gc(a) (a=getchar())
#define Pc(a) putchar(a)
inline int read(){
register int x;register char c(getchar());register bool k;
while(!isdigit(c)&&c^'-')if(Gc(c)==EOF)exit(0);
if(c^'-')k=1,x=c&15;else k=x=0;
while(isdigit(Gc(c)))x=(x<<1)+(x<<3)+(c&15);
return k?x:-x;
}
void wr(register int a){
if(a<0)Pc('-'),a=-a;
if(a<=9)Pc(a|'0');
else wr(a/10),Pc((a%10)|'0');
}
signed const INF(0x3f3f3f3f),NINF(0xc3c3c3c3);
long long const LINF(0x3f3f3f3f3f3f3f3fLL),LNINF(0xc3c3c3c3c3c3c3c3LL);
#define Ps Pc(' ')
#define Pe Pc('\n')
#define Frn0(i,a,b) for(register int i(a);i<(b);++i)
#define Frn1(i,a,b) for(register int i(a);i<=(b);++i)
#define Frn_(i,a,b) for(register int i(a);i>=(b);--i)
#define Mst(a,b) memset(a,b,sizeof(a))
#define File(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout)
#define N (5000010)
int n,p,a[N],x,y,z,ans(LINF);
signed main(){
Rd(n),Rd(p);
Frn1(i,1,n)Rd(a[i]);
Frn_(i,n,1)a[i]-=a[i-1];
while(p--)Rd(x),Rd(y),a[x]+=Rd(z),a[y+1]-=z;
Frn1(i,1,n)ans=min(ans,a[i]+=a[i-1]);
wr(ans),exit(0);
}

Example:

洛谷P3717 [AHOI2017初中组]cover

(这不是暴力就可以了吗)

当然正解就是暴力,所以假设 \(n,m,r\leqslant 1000\) ,再看看这道题

这道题就是典型的区间修改(将可以覆盖的地方 \(+1\) )与单次查询(最后问 \(>0\) 的有几个)

于是以行为单位差分,每一次修改相当于对每一行的元素做区间修改

但是每一行被修改的位置和长度互不相同,具体取决于半径 \(r\) 和行号 \(i\) 到圆心所在行号 \(y\) 的距离 \(|i-y|\)

假设存在 \(e_{0..r}\) 数组, \(e_i\) 表示斜边为 \(r\) ,对边为 \(i\) 的直角三角形的邻边长度的整数部分,即 \(e_i=\lfloor\sqrt{r^2-i^2}\rfloor\)

那么对于第 \(i\) 行,如果圆心是 \((x,y)\) ,半径为 \(r\) ,覆盖的区间就是 \([max(1,x-e_{|i-y|}),min(n,x+e_{|i-y|})]\) (可别忘了处理边界)

最后是如何计算 \(e\) 数组

(当然可以直接使用sqrt())

但是注意到\(e\)数组具有不上升单调性,可以使用一种类似递推的方法做,保证计算\(O(n)\)(虽然说其实完全没有必要啦)

具体做法就请大家自己想啦 (其实是本蒟蒻太懒了) ,上代码

//This program is written by Brian Peng.
#pragma GCC optimize("Ofast","inline","-ffast-math")
#pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include<bits/stdc++.h>
using namespace std;
#define Rd(a) (a=read())
#define Gc(a) (a=getchar())
#define Pc(a) putchar(a)
inline int read(){
register int x;register char c(getchar());register bool k;
while(!isdigit(c)&&c^'-')if(Gc(c)==EOF)exit(0);
if(c^'-')k=1,x=c&15;else k=x=0;
while(isdigit(Gc(c)))x=(x<<1)+(x<<3)+(c&15);
return k?x:-x;
}
void wr(register int a){
if(a<0)Pc('-'),a=-a;
if(a<=9)Pc(a|'0');
else wr(a/10),Pc((a%10)|'0');
}
signed const INF(0x3f3f3f3f),NINF(0xc3c3c3c3);
long long const LINF(0x3f3f3f3f3f3f3f3fLL),LNINF(0xc3c3c3c3c3c3c3c3LL);
#define Ps Pc(' ')
#define Pe Pc('\n')
#define Frn0(i,a,b) for(register int i(a);i<(b);++i)
#define Frn1(i,a,b) for(register int i(a);i<=(b);++i)
#define Frn_(i,a,b) for(register int i(a);i>=(b);--i)
#define Mst(a,b) memset(a,b,sizeof(a))
#define File(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout)
#define N (110)
int n,m,r,a[N][N],e[N],x,y,ans;
signed main(){
Rd(n),Rd(m),*e=Rd(r);
Frn1(i,1,r){e[i]=e[i-1];while(e[i]*e[i]+i*i>r*r)--e[i];}
while(m--){
Rd(x),Rd(y);
Frn1(i,max(1,y-r),min(n,y+r)){
++a[max(1,x-e[abs(i-y)])][i];
--a[min(n,x+e[abs(i-y)])+1][i];
}
}
Frn1(i,1,n)Frn1(j,1,n)if(a[i][j]+=a[i-1][j])++ans;
wr(ans),exit(0);
}

(个人认为如果\(n,m,r\leqslant 1000\),这道题至少上绿)

到此为止差分的所有基本操作都讲完啦!


Conclusion & Extension:

差分在区间修改和单词查询问题上,有着 \(O(n)-O(1)-O(n)\) 的优秀复杂度

它的基本思想就是用一个数字的变化表示一段区间的整体变化,达到优化的目的

其实,大家有木有感觉到这个差分有点像离散版的微分(把 \(\Delta x\) 设为 \(1\) )?

当差分遇上数据结构,一段美妙的旅程又将开始:树状数组(嘻嘻嘻又是我)

到此为止本篇文章就圆满结束啦,请各位奆佬们多多指教和支持,THX!

Basic Thought / Data Structure: 差分 Difference的更多相关文章

  1. Basic Thought / Data Structure: 前缀和 Prefix Sum

    Intro: 在OI中,前缀和是一种泛用性很高的数据结构,也是非常重要的优化思想 Function: 求静态区间和 模板题:输入序列\(a_{1..n}\),对于每一个输入的二元组\((l,r)\), ...

  2. hdu-5929 Basic Data Structure(双端队列+模拟)

    题目链接: Basic Data Structure Time Limit: 7000/3500 MS (Java/Others)    Memory Limit: 65536/65536 K (Ja ...

  3. HDU 5929 Basic Data Structure 模拟

    Basic Data Structure Time Limit: 7000/3500 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Oth ...

  4. HDU 5929 Basic Data Structure 【模拟】 (2016CCPC东北地区大学生程序设计竞赛)

    Basic Data Structure Time Limit: 7000/3500 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Oth ...

  5. Basic Data Structure

    Basic Data Structure Time Limit: 7000/3500 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Oth ...

  6. Basic Data Structure HDU - 5929 (这个模拟我要报警了)

    Mr. Frog learned a basic data structure recently, which is called stack.There are some basic operati ...

  7. Finger Trees: A Simple General-purpose Data Structure

    http://staff.city.ac.uk/~ross/papers/FingerTree.html Summary We present 2-3 finger trees, a function ...

  8. hdu 4217 Data Structure? 树状数组求第K小

    Data Structure? Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) ...

  9. Python: tree data structure

    # 树结构 from pythonds.basic.stack import Stack #pip install pythonds from pythonds.trees.binaryTree im ...

随机推荐

  1. Java 注解与单元测试

    注解 Java注解是在JDK1.5 之后出现的新特性,用来说明程序的,注解的主要作用体现在以下几个方面: 编译检查,例如 @Override 编写文档,java doc 会根据注解生成对应的文档 代码 ...

  2. 记录安装Python第三方包“tesserocr”的方法和遇到的坑

    1. 环境: 系统环境:Win7 32 位系统 Python版本: 3.6.5        虚拟环境为:Miniconda3 2. 共需要安装的模块: a. tesserocr b. tessera ...

  3. promethues安装

    prometheus 1. 下载安装 下载安装:https://github.com/prometheus/prometheus/releases/tag/v2.9.2 wget https://gi ...

  4. python接口自动化测试-unittest-生成测试报告

    用例的管理问题解决了后,接下来要考虑的就是报告我问题了,这里生成测试报告主要用到 HTMLTestRunner.py 这个模块,下面简单介绍一下如何使用: 一.下载HTMLTestRunner下载: ...

  5. linux入门基础指令大全(汇总)

    一.文件目录指令 1 pwd指令 pwd 显示当前所在的目录 2 ls指令 ls [选项] [目录或文件] 查看文件信息 ls -a 查看所有文件和目录,包括隐藏的 ls -l 以列表的方式显示 ll ...

  6. 为什么大家都说Java中只有值传递?

    最近跟Java中的值传递和引用传递杠上了,一度怀疑人生.查了很多资料,加上自己的理解,终于搞清楚了,什么是值传递和引用传递.也搞明白了,为什么大家都说Java只有值传递,没有引用传递.原来,我一直以来 ...

  7. 扫描器是如何判定有xss漏洞的

    这个问题,看似简单,实则触及到很多人的知识盲区 我们都知道,弹窗就能判定这个页面存在xss, 那么扫描器是怎么判断的呢,或者说扫描器是怎么判断是否弹窗的呢 测试发现 当响应的头中content-typ ...

  8. 【tf.keras】AdamW: Adam with Weight decay

    论文 Decoupled Weight Decay Regularization 中提到,Adam 在使用时,L2 与 weight decay 并不等价,并提出了 AdamW,在神经网络需要正则项时 ...

  9. Linux学习_菜鸟教程_2

    Linux 系统目录 /bin:  bin是Binary的缩写,这个目录存放着最经常使用的命令. /boot:  存放启动Linux时的一些核心文件,包括一些连接文件以及镜像文件. /dev : de ...

  10. hexo搭建博客系列(三)美化主题

    文章目录 其他搭建 1. 添加博客图标 2. 鼠标点击特效(二选一) 2.1 红心特效 2.2 爆炸烟花 3. 设置头像 4. 侧边栏社交小图标设置 5. 文章末尾的标签图标修改 6. 访问量统计 7 ...