前缀和与差分小结

前缀和

一维数组

前缀和其实就相当于数列中的前n项和，在m次询问求某个子序列的和时，如果每次都遍历序列，则时间复杂度为O(n*m)，而若先统计前缀和，则时间复杂度为O(n+m)。关键步骤如下所示：

    //输入
    for(int i = 1; i <= n; i ++) scanf("%d", &a[i]);
    //计算前缀和
    for(int i = 1; i <= n; i ++) s[i] = s[i - 1] + a[i];
    while(m -- ){//m次询问
        int l, r;
        cin >> l >> r ;
        cout << s[r] - s[l - 1] << endl;//输出[l,r]区间的和
    }

二维数组

如上为对于一维数组，而对于二维数组/矩阵，同样可以采用前缀和方法，可进行类比。考虑这样一个问题，对于一个n*m的矩阵，给出q个询问，每次计算(x1, y1), (x2, y2)内的子矩阵和。

#include <iostream>
using namespace std;
const int N = 1010;
int a[N][N], s[N][N];
int main(){
    int n, m, q;
    cin >> n >> m >> q;
    for(int i = 1; i <= n; i ++ )
       for(int j = 1; j <= m; j ++ ){
           //边输入边计算前缀和
           scanf("%d", &a[i][j]);
           s[i][j] = a[i][j] + s[i][j - 1] + s[i - 1][j] - s[i - 1][j - 1];
       }
    while(q -- ){
        int x1, y1, x2, y2;
        scanf("%d%d%d%d",&x1, &y1, &x2, &y2);
        //输出子矩阵和
        printf("%d\n", s[x2][y2] - s[x2][y1 - 1] - s[x1 - 1][y2] + s[x1 - 1][y1 - 1]);
    }
    return 0;
}

如上，关键在于计算前缀和与输出子矩阵和，这里体现了容斥原理的思想，也类似概率论课程中求概率分布的过程（只是这里求的是离散值）。
tips：在用前缀和或差分法求解问题时，矩阵下标从1开始，这是为了方便计算，下标为0的元素自动初始化为0。

差分

一维数组

差分其实类似于前缀和的逆运算，可以由差分数组的前缀和求出原数组的元素。什么时候运用差分的思想与性质可以减少时间复杂度呢，考虑这样一个问题，对一个长度为n的整数序列，输入m个操作，每个操作包含三个整数l, r, c，表示将序列中[l, r]之间的每个数加上c，然后输出进行完所有操作后的序列。
对数组[l, r]区间上的所有数加上C，其实等价于对其差分数组b进行如下操作：

因为差分数组的前缀和为原数组，在位置l上加C相当于对l之后的所有数都加C，因此还需要在r+1位置减去C，从而对数组a的多个数的操作转化为了对其差分数组b中两个数的操作。具体代码如下：

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

这样，最后只需对差分数组b求和后，即可得原数组a。
而在刚开始初始化差分数组b时，可以设想刚开始a,b均为0数组，此时一一插入a数组的每个数，即可得到a数组的差分数组b。而该插入操作相当于对特殊区间[i, i]加a[i]，如下所示：

for(int i = 1; i <= n; i ++ ) insert(i, i, a[i]);

完整代码如下：

#include <iostream>
using namespace std;
const int N = 1e5 + 10;
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);
    for(int i = 1; i <= n; i ++ ) scanf("%d", &a[i]);
    //得到a的差分数组
    for(int i = 1; i <= n; i ++ ) insert(i, i, a[i]);
    while(m -- ){//进行m次操作，转化为对差分数组b的操作
        int l, r, c;
        scanf("%d%d%d", &l, &r, &c);
        insert(l, r, c);
    }
    //由差分数组b求和得到原数组a
    for(int i = 1; i <= n; i ++ ) b[i] = b[i - 1] + b[i];
    for(int i = 1; i <= n; i ++ ) printf("%d ", b[i]);
    return 0;
}

总结：初始时要得到a的差分数组b也可以利用差分的原始定义，即通过b[1]=a[1], b[2]=a[2]-a[1], b[3]=a[3]-a[2],…得到。通过这种差分的方法，时间复杂度显然减小，n与m的关系由乘积变为加法。

二维数组

有了一维数组的基础，二维数组通过类比与画图也可以较好地解决。在对子矩阵的加C操作中，等价于对其差分矩阵做如下操作：

代码表示如下：

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;
}

同样在对原始矩阵a求其差分时，设想a矩阵每个元素一一插入：

insert(i, j, i, j, a[i][j]);

最后，对差分矩阵b进行求前缀和操作，然后输出，完整代码及注释如下：

#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(){
    scanf("%d%d%d", &n, &m, &q);
    //边输入a矩阵边计算其差分矩阵b
    for(int i = 1; i <= n; i ++)
      for(int j = 1; j <= m; j ++ ){
          scanf("%d", &a[i][j]);
          insert(i, j, i, j, a[i][j]);//计算差分矩阵b
      }
    //q次操作
    while(q -- ){
        int x1, y1, x2, y2, c;
        cin >> x1 >> y1 >> x2 >> y2 >> c;
        insert(x1, y1, x2, y2, c);
    }
    //对差分矩阵b求前缀和得到q次操作后的a矩阵
    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;
}

总结：差分思想解题，虽然多开辟了一个数组b，但其减小了时间复杂度，以空间换取了时间。

寒假复习数据结构与算法，同时刷一些题目，然后刚开始试着写一点题解，请各位大佬指点，共同学习。