差分矩阵
原题链接
https://www.acwing.com/problem/content/800/
初始化的公式
差分跟前缀和的关系,差分是前缀和的逆运算
前缀和是将原本所有的数列都当成是一堆差分,然后对其进行求前缀和
差分是将原本所有的数列每一个数都当成一个前缀和,然后求他们的差分
而他们所使用的公式是一样的,只是所求的变量是不一样的,前缀和是根据原数列的数求前缀和,而现在是根据前缀和求原数列的数
差分矩阵初始化公式
由上图的前缀和公式可得,s[x][y] = (s[x][y-1] - s[x-1][y-1]) + s[x-1][y-1] + (s[x-1][y] - s[x-1][y-1]) + a[x][y]
数组s表示各个坐标的前缀和,而差分中原数列的每个坐标的值都表示差分的前缀和,所以,可以将原数列的数当成前缀和公式中的s数组,而a数组的数表示前缀和的数
所以a[x][y] = s[x][y] -(s[x][y-1] - s[x-1][y-1]) - s[x-1][y-1] - (s[x-1][y] - s[x-1][y-1])
从运算的式子上面来看,就是互逆运算。
差分数列初始化公式
应用的公式
快速对子矩阵进行增加或者减去相同的数
差分矩阵公式
给以(x1, y1)为左上角,(x2, y2)为右下角的子矩阵中的所有元素加上c:
b[x1][y1] += c;
b[x1][y2+1] -= c;
b[x2+1][y1] -= c;
b[x2+1][y2+1] += c;
数列差分公式
给区间[l, r]中的每个数加上c:
B[l] += c, B[r + 1] -= c
代码
#include<iostream>
using namespace std;
const int N = 1010;
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[x1][y2+1] -= c;
b[x2+1][y1] -= c;
b[x2+1][y2+1] += c;
}
int main()
{
cin>>n>>m>>q;
for(int i = 1; i <= n; i++)
for(int j = 1; j <= m; j++)
scanf("%d",&a[i][j]);
//求差分矩阵
//对当前位置加上a[i][j] 然后再对其后面的数的减去a[i][j] 相当于只对当前位置增加
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][j-1] + b[i-1][j] - b[i-1][j-1];
//输出答案
for(int i = 1; i <= n; i++)
{
for(int j = 1; j <= m; j++)cout<<b[i][j]<<' ';
cout<<endl;
}
return 0;
}
for(int j = 1; j <= m; j++)cout<<b[i][j]<<' ';
cout<<endl;
}
return 0;
}