（模板）差分&&前缀和

差分数组在acm中有很广泛的应用。对于原数组a：1 2 3 5 5，其差分数组就是sub：1 1 1 2 0，即每一项与前一项的差。其性质有：

• 差分数组求前缀和能得到原数组
• 对区间[l,r]上都加上d在差分数组上表现为sub[l]+=d , sub[r+1]-=d
• 进一步的，在[l,r]上加上首项为k，公差为d的等差数列，在差分数组的表现为sub[l]+=k , sub[i]+=d (l+1<=i<=r) , sub[r+1]-=k+(r-l)*d

void add(int l,int r,int a,int k){
    d2[l]+=a;
    d2[l+1]+=k-a;
    d2[r+1]-=(r-l+1)*k+a;
    d2[r+2]-=(l-r)*k-a;
}

求前缀和的操作如下：

void pre_sum(){
    for(int i=1;i<=n;++i)
    {
        d2[i]+=d2[i-1];
    }
}

因此还原数组的操作：

pre_sum();
pre_sum();

离线操作的话复杂度即O(m)。（m为操作次数）

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;

const int maxn=1e5+5;
int n,m,a[maxn],sub[maxn];
int tr[maxn<<2],lz[maxn<<2];

void pushup(int u){
    tr[u]=tr[u<<1]+tr[u<<1|1];
}

void pushdown(int u,int l,int r){
    int mid=(l+r)>>1;
    tr[u<<1]+=lz[u]*(mid-l+1);
    lz[u<<1]+=lz[u];
    tr[u<<1|1]+=lz[u]*(r-mid);
    lz[u<<1|1]+=lz[u];
    lz[u]=0;
}

void build(int u,int l,int r){
    if(l==r){
        tr[u]=sub[l];
        return;
    }
    int mid=(l+r)>>1;
    build(u<<1,l,mid);
    build(u<<1|1,mid+1,r);
    pushup(u);
}

void update(int u,int l,int r,int L,int R,int v){
    if(l>=L&&r<=R){
        tr[u]+=v*(r-l+1);
        lz[u]+=v;
        return;
    }
    if(lz[u]) pushdown(u,l,r);
    int mid=(l+r)>>1;
    if(L<=mid) update(u<<1,l,mid,L,R,v);
    if(R>mid) update(u<<1|1,mid+1,r,L,R,v);
    pushup(u);
}

int query(int u,int l,int r,int L,int R){
    if(l>=L&&r<=R){
        return tr[u];
    }
    if(lz[u]) pushdown(u,l,r);
    int ret=0,mid=(l+r)>>1;
    if(L<=mid) ret+=query(u<<1,l,mid,L,R);
    if(R>mid) ret+=query(u<<1|1,mid+1,r,L,R);
    pushup(u);
    return ret;
}

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)
        sub[i]=a[i]-a[i-1];
    build(1,1,n);
    while(m--){
        int op,l,r,k,d,p;
        scanf("%d",&op);
        if(op==1){
            scanf("%d%d%d%d",&l,&r,&k,&d);
            update(1,1,n,l,l,k);
            if(r>l) update(1,1,n,l+1,r,d);
            if(r!=n) update(1,1,n,r+1,r+1,-k-(r-l)*d);
        }
        else{
            scanf("%d",&p);
            printf("%d\n",query(1,1,n,1,p));
        }
    }
    return 0;
}

原文地址：https://www.cnblogs.com/FrankChen831X/p/12427963.html