#include <cstdio> #include <iostream> using namespace std; const int N=1000010; int n,m,mod; long long a[N],addv[N],mulv[N]; void pushdown(int o,int len){ if(mulv[o]!=1){ mulv[o*2]=mulv[o*2]*mulv[o]%mod; mulv[o*2+1]=mulv[o*2+1]*mulv[o]%mod; a[o*2]=a[o*2]*mulv[o]%mod; a[o*2+1]=a[o*2+1]*mulv[o]%mod; addv[o*2]=addv[o*2]*mulv[o]%mod; addv[o*2+1]=addv[o*2+1]*mulv[o]%mod; mulv[o]=1; } if(addv[o]){ addv[o*2]=(addv[o*2]+addv[o])%mod; addv[o*2+1]=(addv[o*2+1]+addv[o])%mod; a[o*2]=(a[o*2]+addv[o]*(len-len/2))%mod; a[o*2+1]=(a[o*2+1]+addv[o]*(len-len/2))%mod; addv[o]=0; } return; } void pushup(int o){ a[o]=(a[o*2]+a[o*2+1])%mod; return; } void update(int l,int r,int L,int R,long long v,int o){ if(l>R||r<L)return; if(L<=l && r<=R){addv[o]=(addv[o]+v)%mod; a[o]=(a[o]+v*(r-l+1))%mod; return;} pushdown(o,r-l+1); int mid=(l+r)>>1; update(l,mid,L,R,v,o*2); update(mid+1,r,L,R,v,o*2+1); pushup(o); } int query(int l,int r,int L,int R,int o){ int ans=0; if(l>R||r<L)return 0; if(L<=l && r<=R)return a[o]; pushdown(o,r-l+1); int mid=(l+r)>>1; ans+=query(l,mid,L,R,o*2)%mod; ans+=query(mid+1,r,L,R,o*2+1)%mod; return ans%mod; } void mul(int l,int r,int L,int R,int v,int o){ if(l>R||r<L)return; if(L<=l && r<=R){mulv[o]=mulv[o]*v%mod; a[o]=(a[o]*v)%mod; addv[o]=addv[o]*v%mod; return;} pushdown(o,r-l+1); int mid=(l+r)>>1; mul(l,mid,L,R,v,o*2); mul(mid+1,r,L,R,v,o*2+1); pushup(o); } int main(){ int k=0; scanf("%d%d%d",&n,&m,&mod); while((1<<k)<n)k++; k=1<<k; for(int i=k;i<k+n;i++)scanf("%lld",&a[i]); for(int i=k-1;i;i--)a[i]=a[i*2]+a[i*2+1]; for(int i=1;i<k+n;i++)mulv[i]=1; while(m--){ int x; scanf("%d",&x); if(x==2){ int l,r; long long v; scanf("%d%d%lld",&l,&r,&v); v%=mod; update(1,k,l,r,v,1); } else if(x==3){ int l,r; scanf("%d%d",&l,&r); printf("%d ",query(1,k,l,r,1)); } else { int l,r; long long v; scanf("%d%d%lld",&l,&r,&v); v%=mod; mul(1,k,l,r,v,1); } } return 0; }