SP11470 TTM To the moon（线段树动态开点+时光倒流）

洛谷传送门

【题目分析】

洛谷第一道黑题唉。

题目也说的很清楚，维护一个数据结构，支持区间加，区间求和，区间历史求和，时光倒流。

很明显如果将标记下传，那么还会多个撤销操作，这不优秀，所以我们就考虑标记永久化。

然后其他就是主席树的基本操作了。。。。（240ms氵谷rank6）

【代码~】

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int MAXN=1e5+10;

int n,m,tot=0,rt[MAXN<<4],siz=0;
struct Node{
    int l,r;
    LL sum,add;
}tr[MAXN<<4];

LL Read()
{
    LL i=0,f=1;
    char c;
    for(c=getchar();(c>'9'||c<'0')&&c!='-';c=getchar());
    if(c=='-')
      f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar())
      i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

void push_up(int root)
{
    tr[root].sum=tr[tr[root].l].sum+tr[tr[root].r].sum;
}

void build(int &root,int l,int r)
{
    root=++tot;
    if(l==r)
    {
        tr[root].sum=Read();
        return;
    }
    int mid=l+r>>1;
    build(tr[root].l,l,mid);
    build(tr[root].r,mid+1,r);
    push_up(root);
}

void update(int &root,int las,int L,int R,int l,int r,LL v)
{
    tr[root=++tot]=tr[las];
    if(L<=l&&r<=R)
	{
		tr[root].add+=v;
		return;
	}
    tr[root].sum+=v*(min(R,r)-max(L,l)+1);
    int mid=l+r>>1;
    if(R<=mid)
	  update(tr[root].l,tr[las].l,L,R,l,mid,v);
    else 
	{
		if(L>mid)
		  update(tr[root].r,tr[las].r,L,R,mid+1,r,v);
    	else 
		  update(tr[root].l,tr[las].l,L,mid,l,mid,v),update(tr[root].r,tr[las].r,mid+1,R,mid+1,r,v);
    }
}

LL query(int root,int L,int R,int l,int r)
{
    LL ret=tr[root].add*(min(R,r)-max(L,l)+1);
    if(L<=l&&r<=R)
      return tr[root].sum+ret;
    int mid=l+r>>1;
    if(R<=mid)
      return ret+query(tr[root].l,L,R,l,mid);
    if(L>mid)
      return ret+query(tr[root].r,L,R,mid+1,r);
    return ret+query(tr[root].l,L,mid,l,mid)+query(tr[root].r,mid+1,R,mid+1,r);
}

int main()
{
    n=Read(),m=Read();
    build(rt[0],1,n);
    while(m--)
    {
        char cz[2];
        scanf("%s",cz);
        if(cz[0]=='C')
        {
            int l=Read(),r=Read();
            LL v=Read();
            ++siz,update(rt[siz],rt[siz-1],l,r,1,n,v);
        }
        if(cz[0]=='Q')
        {
            int l=Read(),r=Read();
            printf("%lld\n",query(rt[siz],l,r,1,n));
        }
        if(cz[0]=='H')
        {
            int l=Read(),r=Read(),p=Read();
            printf("%lld\n",query(rt[p],l,r,1,n));
        }
        if(cz[0]=='B')
          siz=Read(),tot=rt[siz+1]-1;
    }
    return 0;
}