P3373线段树2

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const ll N=;
ll sum[N<<],lazy1[N<<],lazy2[N<<];
ll a[N];
ll n,m,p;

inline ll read()
{
    ll ans=;
    char last=' ',ch=getchar();
    while(ch<''||ch>'') last=ch,ch=getchar();
    while(ch>=''&&ch<='') ans=ans*+ch-'',ch=getchar();
    if(last=='-') ans=-ans;
    return ans;
}

void build(ll cnt,ll l,ll r)
{
    if(l==r)
    {
        sum[cnt]=a[l]%p;
        return;
    }
    else
    {
        ll mid=(l+r)>>;
        build(cnt<<,l,mid);
        build((cnt<<)|,mid+,r);
        sum[cnt]=(sum[cnt<<]+sum[(cnt<<)|])%p;
    }
}

inline bool cover(ll nl,ll nr,ll l,ll r)
{
    return l<=nl&&r>=nr;
}

inline bool intersection(ll nl,ll nr,ll l,ll r)
{
    return l<=nr&&r>=nl;
}

void pushdown(ll cnt,ll l,ll r)
{
    ll mid=(l+r)>>;
    lazy1[cnt<<]=(lazy1[cnt<<]*lazy2[cnt]+lazy1[cnt])%p;
    lazy1[(cnt<<)|]=(lazy1[(cnt<<)|]*lazy2[cnt]+lazy1[cnt])%p;
    lazy2[cnt<<]=(lazy2[cnt<<]*lazy2[cnt]%p)%p;
    lazy2[(cnt<<)|]=(lazy2[(cnt<<)|]*lazy2[cnt]%p)%p;
    sum[cnt<<]=((lazy1[cnt]*(mid-l+)%p)+(sum[cnt<<]*lazy2[cnt]%p))%p;
    sum[cnt<<|]=((lazy1[cnt]*(r-mid)%p)+(sum[cnt<<|]*lazy2[cnt]%p))%p;
    lazy2[cnt]=;
    lazy1[cnt]=;
}

void add1(ll cnt,ll nl,ll nr,ll l,ll r,ll a)
{
    if(cover(nl,nr,l,r))
    {
        sum[cnt]=(sum[cnt]+(nr-nl+)*a%p)%p;
        lazy1[cnt]=(lazy1[cnt]+a%p)%p;
        return ;
    }
    pushdown(cnt,nl,nr);
    ll mid=(nl+nr)>>;
    if(intersection(nl,mid,l,r)) add1(cnt<<,nl,mid,l,r,a);
    if(intersection(mid+,nr,l,r)) add1(cnt<<|,mid+,nr,l,r,a);
    sum[cnt]=(sum[cnt<<]+sum[cnt<<|])%p;
}

void add2(ll cnt,ll nl,ll nr,ll l,ll r,ll a)
{
    if(cover(nl,nr,l,r))
    {
        sum[cnt]=(sum[cnt]*a)%p;
        lazy1[cnt]=(lazy1[cnt]*a)%p;
        lazy2[cnt]=(lazy2[cnt]*a)%p;
        return ;
    }
    pushdown(cnt,nl,nr);
    ll mid=(nl+nr)>>;
    if(intersection(nl,mid,l,r)) add2(cnt<<,nl,mid,l,r,a);
    if(intersection(mid+,nr,l,r)) add2(cnt<<|,mid+,nr,l,r,a);
    sum[cnt]=(sum[cnt<<]+sum[cnt<<|])%p;
}

ll query(ll cnt,ll nl,ll nr,ll l,ll r)
{
    if(cover(nl,nr,l,r))
    {
        return sum[cnt]%p;
    }
    pushdown(cnt,nl,nr);
    ll mid=(nl+nr)>>;
    ll ans=;
    if(intersection(nl,mid,l,r)) ans+=query(cnt<<,nl,mid,l,r)%p;
    if(intersection(mid+,nr,l,r)) ans+=query(cnt<<|,mid+,nr,l,r)%p;
    return ans;
}

int main()
{
    n=read(),m=read(),p=read();
    for(ll i=;i<=n;i++)
        a[i]=read();
    for(int i=;i<=N<<;i++) lazy2[i]=;
    build(,,n);
    while(m--)
    {
        ll k;
        k=read();
        if(k==)
        {
            ll l=read(),r=read(),t=read();
            add2(,,n,l,r,t);
        }
        if(k==)
        {
            ll l=read(),r=read(),t=read();
            add1(,,n,l,r,t);
        }
        if(k==)
        {
            ll l=read(),r=read();
            printf("%lld\n",query(,,n,l,r)%p);
        }
    }
}