求调！！！

#include <bits/stdc++.h>
using namespace std;
long long n,m,a[100001],ls,l,r,z,k,q;
struct ab{
    int l,r;
    long long z,lb,cb=1;
} t[1000001];
long long j(int x,int l,int r){
    t[x].l=l;
    t[x].r=r;
    if(l==r){
        t[x].z=a[l];
    }else{
        int mid=l+r>>1;
        t[x].z=j(x*2,l,mid)+j(x*2+1,mid+1,r);
    }
    t[x].z%=q;
    return t[x].z;
}
void gx(int x){
	t[x*2].lb+=t[x].lb;
	t[x*2].cb*=t[x].cb;
	t[x*2].z*=t[x].cb;
	t[x*2].z+=t[x].lb*(t[x*2].r-t[x*2].l+1);
	t[x*2].z%=q;
	t[x*2+1].lb+=t[x].lb;
	t[x*2+1].cb*=t[x].cb;
	t[x*2+1].z*=t[x].cb;
	t[x*2+1].z+=t[x].lb*(t[x*2+1].r-t[x*2+1].l+1);
	t[x*2+1].z%=q;
	t[x].lb=0;
	t[x].cb=1;
}
long long qj(int x,int l,int r,long long z){
    if ((l==t[x].l)&&(r==t[x].r)){
        t[x].lb+=z;
        t[x].lb%=q;
        t[x].z+=z*(r-l+1);
        t[x].z%=q;
        return z*(r-l+1)%q;
    }if(t[x].lb||(t[x].cb!=1)){
        gx(x);
	}if (t[x*2].r>=r){
        long long ls=qj(x*2,l,r,z);
        t[x].z+=ls;
        t[x].z%=q;
        return ls%q;
    }else if (l>=t[x*2+1].l){
        long long ls=qj(x*2+ 1, l, r, z);
        t[x].z += ls;
        t[x].z%=q;
        return ls%q;
    }else{
        long long ls = qj(x * 2, l, t[x * 2].r, z), lls = qj(x * 2 + 1, t[x * 2 + 1].l, r, z);
        t[x].z += ls;
        t[x].z += lls;
        t[x].z%=q;
        return (ls + lls)%q;
    }
}
long long qc(int x, int l, int r, long long z){
    if ((l == t[x].l)&& (r == t[x].r)){
    	int ls=t[x].lb,lls=t[x].z;
        t[x].lb *= z;
        t[x].z*=z;
        t[x].cb*=z;
        t[x].z%=q;
        t[x].lb%=q;
        t[x].cb%=q;
        return (ls+lls)*(z-1)%q;
    }if(t[x].lb||(t[x].cb!=1)){
        gx(x);
	}if (t[x * 2].r >= r){
        long long ls = qc(x*2,l,r,z);
        t[x].z += ls;
        t[x].z%=q;
        return ls%q;
    }else if (l >= t[x * 2 + 1].l){
        long long ls = qc(x * 2 + 1, l, r, z);
        t[x].z += ls;
        t[x].z%=q;
        return ls%q;
    }else{
        long long ls = qc(x * 2, l, t[x * 2].r, z), lls = qc(x * 2 + 1, t[x * 2 + 1].l, r, z);
        t[x].z += ls;
        t[x].z += lls;
        t[x].z%=q;
        return (ls + lls)%q;
    }
}	
void c(int x,int l,int r){
    if ((t[x].l==l) && (t[x].r==r)){
        k+=t[x].z;
        k%=q;
        return;
    }if(t[x].lb||t[x].cb){
        gx(x);
	}if (t[x*2].r>=r){
        c(x*2,l,r);
    }else if (l>=t[x*2+1].l){
        c(x*2+1,l,r);
    }else{
        c(x*2,l,t[x*2].r);
        c(x*2+1,t[x*2+1].l,r);
    }
}
int main(){
    cin>>n>>m>>q;
    for(int i=1;i<=n;i++)
        cin>>a[i];
    j(1,1,n);
    while (m--){
        cin>>ls>>l>>r;
        if (ls==1){
            cin>>z;
            qc(1,l,r,z);
        }else if(ls==2){
            cin>>z;
            qj(1,l,r,z);
		}else{
            k=0;
            c(1,l,r);
            cout<<k<<"\n";
        }
    }
    return 0;
}