可持久化线段树——区间更新hdu4348

和线段树类似，每个结点也要打lazy标记

但是lazy标记和线段树不一样

具体区别在于可持久化后lazy-tag不用往下传递，而是固定在这个区间并不断累加，变成了这个区间固有的性质(有点像分块的标记了)

update就按照这么来

int update(int last,int L,int R,int c,int l,int r){
    int now=++size;
    T[now]=T[last];

    if(L<=l && R>=r){
        T[now].sum+=(r-l+)*c;
        T[now].add+=c;
        return now;
    }

    int mid=l+r>>;
    if(L<=mid)T[now].lc=update(T[last].lc,L,R,c,l,mid);
    if(R>mid)T[now].rc=update(T[last].rc,L,R,c,mid+,r);
    pushup(l,r,now);
    return now;
}

查询时由于lazytag固定在区间上。所以向下查询的时候要把上层的lazytag的影响都算上，即递归时传递一个上层区间的  影响值（例如add）

ll query(int now,int L,int R,int add,int l,int r){
    if(L<=l && R>=r) return T[now].sum+(ll)add*(r-l+);
    int mid=l+r>>;
    ll res=;add+=T[now].add;
    if(L<=mid)res+=query(T[now].lc,L,R,add,l,mid);
    if(R>mid)res+=query(T[now].rc,L,R,add,mid+,r);
    return res;
}

此外还有合并维护时，由于子区间没有收到父区间的影响，所以合并时还要算父区间的lazytag

void pushup(int l,int r,int rt){T[rt].sum=T[T[rt].lc].sum+T[T[rt].rc].sum+T[rt].add*(r-l+);}

最后是完整代码，其实本题版本回滚时还可以吧size往回滚，以此节省内存

/*
主席树区间更新
*/
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define maxn 100005
ll n,m,a[maxn];
struct Node{int lc,rc;ll sum,add;}T[maxn*];
int size,rt[maxn];
void pushup(int l,int r,int rt){T[rt].sum=T[T[rt].lc].sum+T[T[rt].rc].sum+T[rt].add*(r-l+);}
int build(int l,int r){
    int now=++size;
    if(l==r){
        T[now].lc=T[now].rc=;
        T[now].sum=a[l];
        return now;
    }
    int mid=l+r>>;
    T[now].lc=build(l,mid);
    T[now].rc=build(mid+,r);
    pushup(l,r,now);
    return now;
}
int update(int last,int L,int R,int c,int l,int r){
    int now=++size;
    T[now]=T[last];

    if(L<=l && R>=r){
        T[now].sum+=(r-l+)*c;
        T[now].add+=c;
        return now;
    }

    int mid=l+r>>;
    if(L<=mid)T[now].lc=update(T[last].lc,L,R,c,l,mid);
    if(R>mid)T[now].rc=update(T[last].rc,L,R,c,mid+,r);
    pushup(l,r,now);
    return now;
}
ll query(int now,int L,int R,int add,int l,int r){
    if(L<=l && R>=r) return T[now].sum+(ll)add*(r-l+);
    int mid=l+r>>;
    ll res=;add+=T[now].add;
    if(L<=mid)res+=query(T[now].lc,L,R,add,l,mid);
    if(R>mid)res+=query(T[now].rc,L,R,add,mid+,r);
    return res;
}
void init(){
    size=;
    memset(rt,,sizeof rt);
    memset(T,,sizeof T);
}
int main(){
    while(scanf("%lld%lld",&n,&m)==){
        init();
        for(int i=;i<=n;i++)scanf("%lld",&a[i]);

        int cur=,l,r,c;char op[];
        rt[cur]=build(,n);
        while(m--){
            scanf("%s",op);
            if(op[]=='C'){scanf("%d%d%d",&l,&r,&c);rt[++cur]=update(rt[cur-],l,r,c,,n);}
            if(op[]=='Q'){scanf("%d%d",&l,&r);cout<<query(rt[cur],l,r,,,n)<<'\n';}
            if(op[]=='H'){
                scanf("%d%d%d",&l,&r,&c);
                cout<<query(rt[c],l,r,,,n)<<'\n';
            }
            if(op[]=='B'){scanf("%d",&c);cur=c;}
        }
    //    puts("");
    }
}