前缀和与差分

1、前缀和

(1)一维前缀和

<1>思想

下标从1开始。

<2>代码

#include <iostream>

using namespace std;

const int N = 1e5+10;

int n,m;
int a[N];// 数组
int s[N];// 前缀和

int main(){
    scanf("%d%d",&n,&m);
    
    // 下标从1开始
    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--){
        int l,r;
        scanf("%d%d",&l,&r);
        printf("%d\n",s[r]-s[l-1]); // 区间和的计算
    }
    return 0;
}

(2)二维前缀和

<1>思想

<2>代码

#include <iostream>

using namespace std;

const int N = 1010;

int n,m,q;
int a[N][N];
int s[N][N];

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]);
        }
    }
    
    for(int i=1;i<=n;++i){
        for(int j=1;j<=m;++j){
            s[i][j] = s[i-1][j] + s[i][j-1] - s[i-1][j-1] + a[i][j];// 算前缀和
        }
    }
    
    while(q--){
        int x1,y1,x2,y2;
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        printf("%d\n",s[x2][y2]-s[x1-1][y2]-s[x2][y1-1]+s[x1-1][y1-1]); // 算子矩阵的和
    }
    return 0;
}

2、差分

可以通过差分方程来求原数组（就是对差分进行求和）。

(1) 一维差分

作用：对原数组中的某个连续区间加上某个固定值。

<1>思想

<2>代码

#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]);
    
    // 由原数组构建差分方程
    for(int i=1;i<=n;++i) insert(i,i,a[i]);
    
    while(m--){
        int r,l,c;
        scanf("%d%d%d",&l,&r,&c);
        insert(l,r,c);
    }
    
    // 由差分方程求前缀和就可以得到原数组
    // a[i]=b[i]+b[i-1]
    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;
}

(1) 二维差分

<1>思路

<2>代码

#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[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]);
        }
    }
    
    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]);
        }
        printf("\n");
    }
}