【1.基础算法-模板】5.前缀和与差分 (模板)

前缀和一维

#include<bits/stdc++.h>
using namespace std;
const int N=100000+10;
long long s[N];//全局变量s[0]=0;
int main()
{
    int a[N],m,n,l,r;
    cin>>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--)
    {
        scanf("%d%d",&l,&r);
        printf("%lld\n",s[r]-s[l-1]);//区间求和
    }
}

前缀和二维

#include<bits/stdc++.h>
using namespace std;
const int maxn=1000+10;
int a[maxn][maxn];
long long sum[maxn][maxn];
int main()
{
    int n,m,k;
    cin>>n>>m>>k;
    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++)
    sum[i][j]=sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1]+a[i][j];//前缀和
    
    while (k--)
    {
        int x1,x2,y1,y2;
        cin>>x1>>y1>>x2>>y2;
        printf("%lld\n",sum[x2][y2]-sum[x2][y1-1]-sum[x1-1][y2]+sum[x1-1][y1-1]);//子矩阵求和    
        
    }
    
    
}

差分一维

#include<bits/stdc++.h>
using namespace std;
const int N=100000+10;
int b[N],a[N],s[N];
void inster(int l,int r,int c)//插入操作
{
    b[l]=b[l]+c;
    b[r+1]=b[r+1]-c;
}
int main()
{
    int n,m;
    cin>>n>>m;
    for(int i=1;i<=n;i++)//输入原数组
    scanf("%d",&a[i]);
    
    for(int i=1;i<=n;i++)//更新差分数组
    inster(i,i,a[i]);
    
    while(m--)
    {
     int l,r,c; 
     cin>>l>>r>>c;
     inster(l,r,c);//更新差分数组
    }
    
    for(int i=1;i<=n;i++)
    s[i]=s[i-1]+b[i];//差分数组的前缀和=原数组
    
    for(int i=1;i<=n;i++)
    printf("%d ",s[i]);//更新后的原数组
}

差分二维

#include<bits/stdc++.h>
using namespace std;
const int N=1000+10;
long long b[N][N],a[N][N];

void inster(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,k;
    cin>>n>>m>>k;
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    scanf("%lld ",&a[i][j]);
    
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    inster(i,j,i,j,a[i][j]);
    
    while(k--)
    {
        int x1,x2,y1,y2,c;
        scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&c);
        inster(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++)
    {
    printf("%lld ",b[i][j]);
    }
    printf("\n");
    }
    return 0;
}