ACM-差分算法

欢迎来到@一夜看尽长安花 博客，您的点赞和收藏是我持续发文的动力

对于文章中出现的任何错误请大家批评指出，一定及时修改。有任何想要讨论的问题可联系我：3329759426@qq.com 。发布文章的风格因专栏而异，均自成体系，不足之处请大家指正。

    专栏：

• java全栈
• C&C++
• PythonAI
• PCB设计
• Linux云计算&运维

文章概述：ACM算法 ——差分

关键词：ACM 差分 一维差分 二维差分

本文目录：

差分

一维差分

原理

运算逻辑（目的）

示例：

二维差分

原理

运算步骤

特点分析

示例：

差分

相当于前缀和的逆运算

 

一维差分

原理

构造前一个b数组，使得b数组的前n项和为b

也就是：b为差分，a为b的前缀和

 

运算逻辑（目的）

原本如果想要使 ai (i=1,2,3,4....) 每个数都增加个常数C,需要计算每个ai，计算时间复杂度为O(n)

但如果用差分的思想，将a转化为b的前缀和，在某个b上加个C，对所有的b求前缀和就相当于a加C了，这时候由于只在一个数上加了个值，时间复杂度就转化为了O(1)

具体逻辑如下：

区间[ l , r ], 在第l个数上加上C 那么前 l 项和 到前r项和这段区间上的 ai 均加上了C

在第r+1 个数上减去C就能够和前面第l个数上加上C抵消，也就是说从第r+1项和往后的a都没变

这样就实现了区间化的增值变化

假设初始化 ai = 0 ,那么易得bi=0；

若是ai != 0 ,即ai 有非零值，就相当于 0 + ai (即 insert[i,i]+ai ),相当于与在bi加上数ai , bi+1上减去ai, 求和前bi项和即为ai，从bi+1项往后即为0，因为一加一减正好抵消

 

示例：

#include <iostream>

using namespace std;
const int N = 100010;

//输入一个长度为 n 的整数序列, 输入 m次操作
int n,m;
int a[N],b[N];

void insert(int l,int r,int c)
{
	b[l] +=c;
	b[r+1]-=c;
}

int main()
{
	scanf("%d%d",&n,&m);
	//获取a数组
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	//初始化b数组
	for(int i=1;i<=n;i++) insert(i,i,a[i]);
	
	while(m--)
	{
		int l,r,c;
		scanf("%d%d%d",&l,&r,&c);
		insert(l,r,c);
	}
	//输出进行完所有操作后的序列。
	for(int i=1;i<=n;i++) b[i]+=b[i-1];
	for(int i=1;i<=n;i++) printf("%d ",b[i]);
	
	return 0;
}

二维差分

原理

与一维差分类似，逻辑上就是构造一个差分矩阵b使得 矩阵中部分区间的前缀和为a

运算步骤

在（x1,y1）插入一个值，那么以其为左上顶点的矩阵前缀和都会加上c,

在（x2+1,y1）(x1,y2+1)分别减去c ,那么以（x2+1,y1）(x1,y2+1)为左上顶点的矩阵前缀和就会减去c

在只有蓝色或绿色的区域相当于没变，因为+c和-c抵消了，在红色区间+c-c-c 最终是-c,

在（x2+1,y2+1）加c,就可以和-c抵消啦！

这样下来就只有黑色区域子差分矩阵的前缀和都加上了c

即

 

特点分析

和一位差分类似，将时间复杂度由O(n)转换到了O(1)

 

示例：

#include <iostream>

using namespace std;

const int N=1010;

//n*m的矩阵，查询q次
int n,m,q;
int a[N][N] ,b[N][N];

void insert(int x1,int y1,int x2,int y2,int c)
{
	b[x1][y1] +=c;
	b[x2+1][y1] -=c;
	b[x1][y2+1] -=c;
	b[x2+1][y2+1] +=c;
}

int main()
{
	scanf("%d%d%d",&n,&m,&q);
	
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
				scanf("%d",&a[i][j]);
	
	//初始化b差分矩阵
	for(int i =1;i<=n;i++)
		for(int j=1;j<=m;j++)
			insert(i,j,i,j,a[i][j]);
	
	while(q--)
	{
		int x1, y1 ,x2 ,y2 ,c;
		cin >> x1 >> y1 >> x2 >> y2 >> c;
		insert(x1,y1,x2,y2,c);
	}
	
	//输出操作后的矩阵
	for(int i=1;i<=n; i++)
		for(int j=1;j<=m;j++)
			b[i][j] +=b[i-1][j]+b[i][j-1]-b[i-1][j-1];
		
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++) printf("%d ",b[i][j]);
        //换行
		puts("");
	}
	
	return 0;
}