主席树（可持久化线段树版）

求区间和模板

struct node{
    int l[maxn*20],r[maxn*20];  //区间大小    maxn = 1e5时为20倍，不够就开40
    int li[maxn*20],ri[maxn*20];//左孩子和右孩子的位置
    int sum[maxn*20],add[maxn*20]; //区间和，懒惰标记
    int tot;///数组大小，使用的时候先初始化为0

    void init(int x,int y){
        li[x] = li[y];  l[x] = l[y];
        ri[x] = ri[y];  r[x] = r[y];
        sum[x] = sum[y]; add[x] = add[y];
    }

    void build(int &pos,int L,int R){//建树,也可以建一个空树
        pos = ++tot;
        sum[pos] = 0;
        l[pos] = L;
        r[pos] = R;
        add[pos] = 0;
        if(L == R){
            scanf("%intd",&sum[pos]);
            return;
        }
        int mid = (L+R)>>1;
        build(li[pos],L,mid);
        build(ri[pos],mid+1,R);
        sum[pos] = sum[li[pos]] + sum[ri[pos]];
    }

    void update(int L,int R,int &x,int y,int k){//区间修改，x传入低k棵线段树的根节点，y为第k-1棵线段树的根节点
        x = ++tot;
        init(x,y);
        if(L == l[x] && R == r[x]){
            sum[x] += (R-L+1)*k;
            add[x] += k;
            return;
        }
        sum[x] += (R - L + 1)*k;
        int mid = (l[x]+r[x])>>1;
        if(R <= mid){
            update(L,R,li[x],li[y],k);
        }
        else if(L > mid){
            update(L,R,ri[x],ri[y],k);
        }
        else{
            update(L,mid,li[x],li[y],k);
            update(mid+1,R,ri[x],ri[y],k);
        }
    }

    int query(int L,int R,int pos,int add1){//查询第根节点为pos的树的l,r和，add1为父亲节点的add累加和，初始传入0
        if(L == l[pos] && R == r[pos]){
            return sum[pos] + add1*(R-L+1);
        }
        int mid = (l[pos]+r[pos])/2;
        if(R <= mid){
            return query(L,R,li[pos],add[pos]+add1);
        }
        else if(L > mid){
            return query(L,R,ri[pos],add[pos]+add1);
        }
        else{
            return query(L,mid,li[pos],add[pos]+add1)+query(mid+1,R,ri[pos],add[pos]+add1);
        }
    }
}tre;

 上面那个有点挫，补一个常数小一点的

int l[maxn*40],r[maxn*40],sum[maxn*40];
int mark[maxn],a[maxn],ans[maxn],root[maxn];
int tot;

void build(int &x,int s,int t){
    x = ++tot;
    sum[x] = 0;
    if(s == t) return;
    int mid = (s+t) >> 1;
    build(l[x],s,mid);
    build(r[x],mid+1,t);
}

void update(int &x,int y,int s,int t,int k,int v){
    x = ++tot;
    l[x] = l[y];
    r[x] = r[y];
    sum[x] = sum[y] + v;
    if(s == t) return;
    int mid = (s + t) >> 1;
    if(k <= mid) update(l[x],l[y],s,mid,k,v);
    else update(r[x],r[y],mid+1,t,k,v);
}

int query(int s,int t,int ll,int rr,int x){
    if(s <= ll && rr <= t) return sum[x];
    int mid = (ll + rr) >> 1;
    int res = 0;
    if(s <= mid) res += query(s,t,ll,mid,l[x]);
    if(t > mid) res += query(s,t,mid+1,rr,r[x]);
    return res;
}