数据结构：分块-区间加法与区间求和

块儿外的元素还是暴力

为了快速求答案，要提前预处理每个块儿的元素和

区间修改的时候，不完整的块儿直接修改，顺便更新元素和

完整的块儿还是打标记

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
long long read()
{
    long long x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
const int maxn=50005;
int n,blo;
int bl[maxn];
long long v[maxn],atag[maxn],sum[maxn];
void add(int a,int b,int c)
{
    for(int i=a;i<=min(bl[a]*blo,b);i++)
        v[i]+=c,sum[bl[a]]+=c;
    if(bl[a]!=bl[b])
        for(int i=(bl[b]-1)*blo+1;i<=b;i++)
            v[i]+=c,sum[bl[b]]+=c;
    for(int i=bl[a]+1;i<=bl[b]-1;i++) atag[i]+=c;
}
long long query(int a,int b)
{
    long long ans=0;
    for(int i=a;i<=min(bl[a]*blo,b);i++)
        ans+=v[i]+atag[bl[a]];
    if(bl[a]!=bl[b])
        for(int i=(bl[b]-1)*blo+1;i<=b;i++)
            ans+=v[i]+atag[bl[b]];
    for(int i=bl[a]+1;i<=bl[b]-1;i++)
        ans+=sum[i]+blo*atag[i];
    return ans;
}
int main()
{
    n=read();blo=sqrt(n);
    for(int i=1;i<=n;i++) v[i]=read();
    for(int i=1;i<=n;i++)
    {
        bl[i]=(i-1)/blo+1;
        sum[bl[i]]+=v[i];  //预处理 
    }
    for(int i=1;i<=n;i++)
    {
        int f=read(),a=read(),b=read(),c=read();
        if(f==0) add(a,b,c);
        if(f==1) printf("%d
",query(a,b)%(c+1));
    } 
    return 0;
}