To the moon 主席树 + 标记永久化

传送门

题目描述

给你一串序列，每次对它有四种操作， 1. 将l到r之间所有的数+d，同时时间戳+1。 2. 查询当前时间戳下l到r的所有数的总和。 3…查询时间戳d下的l到r的所有数的和。 4.将时间戳返回至t。

分析

主席树标记永久化的板子题了，今天来学一手

首先这个题如果没有历史版本的话，我们用线段树可以很轻松的解决这个问题，但是加上了可持久化
线段树在处理区间修改的时候是用懒标记，pushdown操作把标记一层一层往下传递，主席树的处理方法和懒标记有点类似，但是如果我们用pushdown传递标记的话，内存占用较大
标记永久化的写法就是当这个修改区间包含当前区间的时候，就给这个区间打上标记，但是不往下传递，最后查询的时候答案加上每个标记的贡献值就可以了

代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <queue>
#include <cstring>
#define debug(x) cout<<#x<<":"<<x<<endl;
#define _CRT_SECURE_NO_WARNINGS
#pragma GCC optimize("Ofast","unroll-loops","omit-frame-pointer","inline")
#pragma GCC option("arch=native","tune=native","no-zero-upper")
#pragma GCC target("avx2")
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> PII;
const int INF = 0x3f3f3f3f;
const int N = 1e5 + 10;
int n,m,idx;
int now;
int root[N];
ll a[N];

struct Node{
    int l,r;
    ll sum,add;
}tr[N * 30];

int build(int l,int r){
    int p = ++idx;
    if(l == r){
        tr[p].sum = a[l];
        return p;
    }
    int mid = l + r >> 1;
    tr[p].l = build(l,mid);
    tr[p].r = build(mid + 1,r);
    tr[p].sum = tr[tr[p].l].sum + tr[tr[p].r].sum;
    return p;
}

void up(int l,int r,int pre,int &pp,int L,int R,ll val){
    if(l > R || r < L) return ;
    tr[++idx] = tr[pre];
    pp = idx;
    if(l >= L && r <= R){
        tr[pp].add += val;
        tr[pp].sum += val * (r - l + 1);
    }
    else{
        int mid = l + r >> 1;
        tr[pp].sum += val * max(0,min(r,R) - max(l,L) + 1);
        up(l,mid,tr[pre].l,tr[pp].l,L,R,val);
        up(mid + 1,r,tr[pre].r,tr[pp].r,L,R,val);
    }
}

ll query(int l,int r,int x,int L,int R,ll sum){
    if(l > R || r < L) return 0;
    if(l >= L && r <= R) return tr[x].sum + sum * (r - l + 1);
    else{
        ll mid = l + r >> 1;
        return query(l,mid,tr[x].l,L,R,sum + tr[x].add) + query(mid + 1,r,tr[x].r,L,R,sum + tr[x].add);
    }
}

void init(){
    idx = 0;
    now = 0;
}

int main(){
    while(~scanf("%d%d",&n,&m)){
        init();
        for(int i = 1;i <= n;i++) scanf("%lld",&a[i]);
        root[0] = build(1,n);
        while(m--){
            char op;
            ll l,r,d;
            scanf(" %c %lld",&op,&l);
            if(op == 'C'){
                scanf("%lld%lld",&r,&d);
                now++;
                up(1,n,root[now - 1],root[now],l,r,d);
            }
            else if(op == 'Q'){
                scanf("%lld",&r);
                printf("%lld\n",query(1,n,root[now],l,r,0));
            }
            else if(op == 'H'){
                scanf("%lld%lld",&r,&d);
                printf("%lld\n",query(1,n,root[d],l,r,0));
            }
            else now = l;
        }
    }
}

/**
*　　┏┓　　　┏┓+ +
*　┏┛┻━━━┛┻┓ + +
*　┃　　　　　　　┃
*　┃　　　━　　　┃ ++ + + +
*  ████━████+
*  ◥██◤　◥██◤ +
*　┃　　　┻　　　┃
*　┃　　　　　　　┃ + +
*　┗━┓　　　┏━┛
*　　　┃　　　┃ + + + +Code is far away from 　
*　　　┃　　　┃ + bug with the animal protecting
*　　　┃　 　 ┗━━━┓ 神兽保佑,代码无bug　
*　　　┃ 　　　　　　 ┣┓
*　　  ┃ 　　　　　 　┏┛
*　    ┗┓┓┏━┳┓┏┛ + + + +
*　　　　┃┫┫　┃┫┫
*　　　　┗┻┛　┗┻┛+ + + +
*/