suoi37 清点更多船只 (卡空间线段树)

sbw巨佬的卡空间方法，把线段树的叶节点只记到长度为16的区间，然后在叶节点上暴力修改查询，这样点数是$frac{N}{8}$的，可以过...

orz

#include<bits/stdc++.h>
#define pa pair<int,int>
#define CLR(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
const int maxn=(1<<20),maxp=(1<<18)+10;

inline ll rd(){
    ll x=0;char c=getchar();int neg=1;
    while(c<'0'||c>'9'){if(c=='-') neg=-1;c=getchar();}
    while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*neg;
}

ll sum[maxp],v[maxn],laz[maxp];
int N,M;

inline void update(int p,int l,int r){
    if(r-l+1<=16){    
        sum[p]=0;
        for(int i=l;i<=r;i++){
            sum[p]+=v[i];
        }
    }else{
        sum[p]=sum[p<<1]+sum[p<<1|1];
    }
}
inline void deal(int p,int l,int r,ll z){
    if(r-l+1<=16){
        for(int i=l;i<=r;i++){
            v[i]+=z,sum[p]+=z;
        }
    }else sum[p]+=z*(r-l+1);
}

inline void pushdown(int p,int l,int r){
    if(!laz[p]||r-l+1<=16) return;
    int a=p<<1,b=p<<1|1,m=l+r>>1;
    deal(a,l,m,laz[p]),deal(b,m+1,r,laz[p]);
    laz[a]+=laz[p],laz[b]+=laz[p];
    laz[p]=0;
}

void build(int p,int l,int r){
    if(r-l+1<=16) update(p,l,r);
    else{
        int m=l+r>>1;
        build(p<<1,l,m);
        build(p<<1|1,m+1,r);
        // printf("!%d %d %d
",p,l,r);
        update(p,l,r);
    }
}

ll query(int p,int l,int r,int x,int y){
    // printf("~~%d %d %d %d %d
",p,l,r,x,y);
    pushdown(p,l,r);
    if(x<=l&&r<=y) return sum[p];
    if(r-l+1<=16){
        ll re=0;
        for(int i=x;i<=y;i++) re+=v[i];
        return re;
    }
    int m=l+r>>1;ll re=0;
    if(x<=m) re=query(p<<1,l,m,x,min(m,y));
    if(y>=m+1) re+=query(p<<1|1,m+1,r,max(x,m+1),y);
    return re;
}

void add(int p,int l,int r,int x,int y,ll z){
    // printf("!!%d %d %d %d %d
",p,l,r,x,y);
    pushdown(p,l,r);
    if(r-l+1<=16){
        for(int i=x;i<=y;i++)
            v[i]+=z,sum[p]+=z;
    }else if(x<=l&&r<=y){
        deal(p,l,r,z);
        laz[p]+=z;
    }else{
        int m=l+r>>1;
        if(x<=m) add(p<<1,l,m,x,min(m,y),z);
        if(y>=m+1) add(p<<1|1,m+1,r,max(m+1,x),y,z);
        update(p,l,r);
    }
}

int main(){
    //freopen("","r",stdin);
    int i;
    N=rd(),M=rd();
    for(i=1;i<=N;i++)
        v[i]=rd();
    N=1<<((int)log2(N)+1);
    build(1,1,N);
    for(i=1;i<=M;i++){
        char s[10];
        scanf("%s",s);
        int l=rd(),r=rd();
        if(s[1]=='d'){
            int d=rd();
            add(1,1,N,l,r,d);
        }else printf("%lld
",query(1,1,N,l,r));
    }
    return 0;
}