#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);

        }

    }

}