前缀和与差分

前缀和的概念

前缀和，顾名思义就是前n项和，这一点对于没有接触过算法的人也十分容易理解。

至于如何计算前缀和，我们这里有一个约定，计算前缀和的数组我们要求第零项等于0，不用于存储数据，这是由于计算前缀和的时候我们使用的是递推公式，以下面的例子来说明这一点

一维数组的前缀和计算方式

int q[N];
int s[N];
for(int i = 1 ; i<=N ; i++)
   s[i] = s[i-1] + q[i];

 可以看到，为了使第一项递推满足，我们需要从第一项开始。

前缀和的意义以及基本应用

前缀和为我们以后要进行的算法学习提供有力的支持。

实际意义上，前缀和也有它的用处，它主要用于处理涉及到数组中某一区域的值的和等问题，比如求区间内[l，r]的值，或者求解有多少个区间满足和为k的倍数等，如果我们用了前缀和处理数据，就不再需要每次都去暴力遍历，而是直接用O(1)的时间复杂度求解[l，r]区间内的和

//区间[l，r]内的和
Sum=s[r]-s[l-1]；

一维的前缀和是很容易理解和计算的，下面我们来介绍复杂一点的二维前缀和。

二维前缀和

定义

S(i,j)等于矩阵中以q[i][j]为右下角的矩形范围内的所有数字的和。

递推的计算前缀和

将s(i，j)拆乘s(i-1，j)+s(i，j-1)+q[i][j]-s(i-1，j-1)就找到了递推公式

计算子矩阵的和 

如图，与一维前缀和类似的计算子矩阵的和

//我们可以用"+="运算来省去额外的数组空间用于存储前缀和
//int q[n][m];
int s[n][m];

for(int i = 1; i<=n; i++)
    for(int j = 1; i<=m; j++)
        s[i][j]+=s[i-1][j]+s[i][j-1]-s[i-1][j-1];
//如果不使用"+="运算的话,只需要将最后一项s[i-1][j-1]改成q[i-1][j-1]

//输出子矩阵以(x1,y1)为左上角(x2,y2)为右下角的和,包含端点
for(int i = 1; i<=n; i++)
    for(int j = 1; i<=m; j++)
        printf("%d", s[x2][y2] - s[x1-1][y2] - s[x2][y1-1] + s[x1-1][y1-1]);

差分

差分的定义

如果我们有如下数组  ，如果 ，那么就称b数组是a数组的差分，可以发现还有的关系，实际上，差分和前缀和是互逆的

差分的意义和实际作用

与前缀和一样，差分也具有它的为后续算法铺垫的意义。

差分处理往往出现在需要使某个数组某个区间[l，r]内的所有数字做相同增减处理的情况下，比如我们要让  改变成，如果不用差分的话，我们就需要遍历数组[l，r]，但是由于差分的存在，我们可以巧妙的完成这样的工作，由于数组a是b的前缀和，我们只需要用O(1)的时间复杂度去改变数组b的值便可实现；差分的构造我们只需要循环一遍insert函数，如下代码：

int a[N];
int b[N];

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

int main()
{
    int l, r, c;
    for(int i = i; i<=N; i++)
        insert(i, i, a[i]);    //这一步就构造了数组b，使得b[n]=a[n]-a[n-1]
    insert(l, r, c); //使下标属于[l，r]的a[]整体+c
    for(int i = 1; i<=N; i++) b[i]+=b[i-1];  //改变后的a[]等于b[]的前缀和
    for（int i =1； i<=N; i++）printf("%d ",b[i]);
}

要想输出a[]的值，我们只需要输出b[]的前缀和即可。

二维差分

由于差分和前缀和是一对互逆的运算，我们不必去思考如何构造差分，我们假设有二维数组a[][]，构造a的差分b[][]，也就是说数组b的前缀和是a[][]，回忆一维差分的构造，我们首先是写insert函数，再类似在区间[i，i]中插入a[i]，那么二维的话，就是类似的在小矩形(i，j)→(i，j)中插入a[i][j]，所以重点是insert函数（实现改变区间内a[][]的值的操作的函数）,根据下面这个图看一下

可以看到，改变b[i][j]的影响跟前缀和是对应着的，所以有以下insert函数

int a[n][m];
int b[n][m];
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;
}
//构造b数组
for(int i=1 ; i<=n; i++)
    for(int j=1; j<=m; j++)
        insert(i,j,i,j,a[i][j]);

同理，如果我们想要输出改变后的子矩阵中的a[i][j]，只需计算b数组的前缀和即可，这里不再赘述二维前缀和了。