前缀和与差分（分析与模板）

前缀和

处理数组公式：

​ s[i]=s[i-1]+num[i]

输出区间和公式:

​ s[r]-s[l-1]

#模板
#include<stdio.h>
int sum[100010]={0};
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    int num[100010];
    for(int i=1;i<=n;i++) scanf("%d",&num[i]);
    
    for(int i=1;i<=n;i++) sum[i]=sum[i-1]+num[i];//初始化数组
    while(m--){
        int l,r;
        scanf("%d %d",&l,&r);
        printf("%d\n",sum[r]-sum[l-1]);
    }
    return 0;
}


/*
3 5 6 7 8 1
sum[i]=sum[i-1]+num[i]

sum[1]=sum[0]+num[1]-----0+3-----3
sum[2]=sum[1]+num[2]-----3+5-----8
sum[3]=sum[2]+num[3]-----8+6-----14
sum[4]=sum[3]+num[4]-----14+7----21
sum[5]=sum[4]+num[5]-----21+8----29
...

2-4
sum[r]-sum[l-1]

sum[4]-sum[1]
21-3=19
*/

二维前缀和

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-kV5S8452-1676283406734)(C:\Users\zhaoh\Pictures\联想截图\联想截图_20221211230530.png)]

公式：
    初始化：s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+a[i][j]
    子矩阵和：s[x2][y2]-s[x1-1][y1]-s[x1][y1-1]+s[x1-1][y1-1]


#include <stdio.h>
int s[N][N];
int main()
{
    int n,m,q;
    scanf("%d%d%d",&n,&m,&q);
    int num[N][N];
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            scanf("%d",&num[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]+num[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;
}

差分

思路：
    A数组为原数组，构造B数组，使得A数组是B的前缀和数组即A[i]=B[1]+B[2]+...+B[i],此时将A数组的从l到r之间的所有数加上数c即可表示为B[l]+c,此时A[l]=B[1]+B[2]+...+B[l]+c,即是A[l]完成了加c，并且注意到A数组后面的数A[l+1]=B[1]+B[2]+...+B[l]+c+B[l+1],即A[l+1]也完成了加c，也就是A[l]之后的所有数都完成了加c操作，此时就需要将A[r+1]后的数据减c即可完成A[l]~A[r]之间的数加c，最后将B数组求其前缀和即能得到A数组经过操作之后的结果。
    
int A[N](前缀和数组),B[N](差分数组);
两者关系：A[i]=B[1]+B[2]+.....+B[i]
void insert(int l,int r,int c){
    B[l]+=c;
    B[r+1]-=c;
}
int main()
{
    int n,m;
    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]);
    //此处理解：B[1]=A[1],B[2]=A[2]-A[1]....
    while(m--){
        int l,r,c;
        scanf("%d%d%d",&l,&r,&c);
        insert(l,r,c);
    }
    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;
}

二维差分

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-taGjeSAf-1676283406736)(C:\Users\zhaoh\Pictures\联想截图\联想截图_20221213133447.png)]

思路：
    初始化即是在当前位置进行插入操作

公式：
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;
}

模板

#include<stdio.h>
const int N = 1010;
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()
{
    int n,m,q;
    scanf("%d%d%d",&n,&m,&q);
    int i,j;
    for(i=1;i<=n;i++){
        for(j=1;j<=m;j++){
            scanf("%d",&a[i][j]);
        }
    }
    for(i=1;i<=n;i++){
        for(j=1;j<=m;j++){
            insert(i,j,i,j,a[i][j]);
        }
    }
    while(q--){
        int x1,y1,x2,y2,c;
        scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&c);
        insert(x1,y1,x2,y2,c);
    }
    for(i=1;i<=n;i++){
        for(j=1;j<=m;j++){
            b[i][j]+=b[i][j-1]+b[i-1][j]-b[i-1][j-1];//求前缀和还原数组
        }
    }
    for(i=1;i<=n;i++){
        for(j=1;j<=m;j++){
            printf("%d ",b[i][j]);
        }
        printf("\n");
    }
    return 0;
}     

    }
}
for(i=1;i<=n;i++){
    for(j=1;j<=m;j++){
        printf("%d ",b[i][j]);
    }
    printf("\n");
}
return 0;

}