差分:指给定一个数组a,对a的每个元素a[i]都存在b[i] = a[i] - a[i - 1]。
例如:
b1 = a1;
b2 = a2 - a1;
b3 = a3 - a2;
…
bn = an-1 - an;
我们可以发现对任何 b1 一直加到 bi 都存在相加的结果等于ai,所以当 b 是 a 的差分数组时,a 也是 b 的前缀和数组。 由于a
是 b 的前缀数组,所以想要对 a 某一区间的值进行多次修改的话,可以用O(1)的复杂度修改 b 数组再求 b 数组的前缀和得到 a 数组。
例如:对于数组 a : 1 2 3 4 5 ,可以得到差分数组 b :1 1 1 1 1; 如果想要对第3位到第4位进行加2操作,l =
3,r = 4,b[l] += 2,b[r + 1] -= 2; 这时的 b 数组:1 1 3 1 -1,进行前缀和,a = 1 2 5 6 5。正好对应之前操作。(进行乘除操作可以类推到前缀乘)
差分的修改操作可以理解为当 b[l] += c 时,数组 a 从 l 开始的每个数都加了 c ,要想 l 到 r 区间 +c,还需要在 b[r + 1] -= c,从而使 r 后面的数不受影响。
来道模板题:
输入一个长度为 n 的整数序列。
接下来输入 m 个操作,每个操作包含三个整数 l,r,c ,表示将序列中 [l,r] 之间的每个数加上 c 。
请你输出进行完所有操作后的序列。
输入格式
第一行包含两个整数 n 和 m 。
第二行包含 n 个整数,表示整数序列。
接下来 m 行,每行包含三个整数 l,r,c ,表示一个操作。
输出格式
共一行,包含 n 个整数,表示最终序列。
数据范围
1≤n,m≤100000 ,
1≤l≤r≤n ,
−1000≤c≤1000 ,
−1000≤整数序列中元素的值≤1000
输入样例:
6 3
1 2 2 1 2 1
1 3 1
3 5 1
1 6 1
输出样例:
3 4 5 3 4 2
#include<bits/stdc++.h>
using namespace std;
const int N = 100010;
int a[N],b[N];
int n,m,l,r,c,x;
int main()
{
cin>>n>>m;
for(int i = 1; i <= n; i ++ )
{
cin>>a[i];
b[i] = a[i] - a[i - 1];//得到差分数组
}
while(m -- )
{
cin>>l>>r>>c;
b[l] += c;
b[r + 1] -= c;//对l到r区间进行区间加c操作
}
for(int i = 1; i <= n ; i ++ )
{
x += b[i];//对差分数组求前缀和得到操作后的原始数组
cout<<x<<" ";
}
return 0;
}
类似的同样有差分矩阵:
输入一个 n 行 m 列的整数矩阵,再输入 q 个操作,每个操作包含五个整数 x1,y1,x2,y2,c ,其中 (x1,y1) 和 (x2,y2) 表示一个子矩阵的左上角坐标和右下角坐标。
每个操作都要将选中的子矩阵中的每个元素的值加上 c 。
请你将进行完所有操作后的矩阵输出。
输入格式
第一行包含整数 n,m,q 。
接下来 n 行,每行包含 m 个整数,表示整数矩阵。
接下来 q 行,每行包含 5 个整数 x1,y1,x2,y2,c ,表示一个操作。
输出格式
共 n 行,每行 m 个整数,表示所有操作进行完毕后的最终矩阵。
数据范围
1≤n,m≤1000 ,
1≤q≤100000 ,
1≤x1≤x2≤n ,
1≤y1≤y2≤m ,
−1000≤c≤1000 ,
−1000≤矩阵内元素的值≤1000
输入样例:
3 4 3
1 2 2 1
3 2 2 1
1 1 1 1
1 1 2 2 1
1 3 2 3 2
3 1 3 4 1
输出样例:
2 3 4 1
4 3 4 1
2 2 2 2
#include<bits/stdc++.h>
using namespace std;
const int N = 1010;
int a[N][N],b[N][N],n,m,q,x1,x2,y1,y2,c;
int main()
{
cin>>n>>m>>q;
for(int i = 1; i <= n; i ++ )
for(int j = 1; j <= m; j ++ )
{
cin>>a[i][j];
b[i][j] = a[i][j] + a[i - 1][j - 1] - a[i - 1][j] - a[i][j - 1];//与求前缀和矩阵符号相反
}
while(q -- )
{
cin>>x1>>y1>>x2>>y2>>c;
b[x1][y1] += c;//表示从a矩阵从x1 y1 往后的每个数都加了c
b[x1][y2 + 1] -= c;
b[x2 + 1][y1] -= c;
b[x2 + 1][y2 + 1] += 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];//前缀和
cout<<b[i][j]<<" ";
}
cout<<endl;
}
return 0;
}