#include<bits/stdc++.h>
using namespace std;
#define int long long
int n,m,MOD,a[100030],cnt[400030],f1[400030],f2[400030];
inline void build(int x,int l,int r)
{
	f1[x]=1;
	if(l==r)
	{
		cnt[x]=cnt[l];
		return;
	}
	int mid=(l+r)/2;
	build(x*2,l,mid);
	build(x*2+1,mid+1,r);
	cnt[x]=cnt[x*2]+cnt[x*2+1];
}
inline void change(int x,int l,int r,int k1,int k2)
{
	cnt[x]=(k1*cnt[x]+(r-l+1)*k2%MOD)%MOD;
	f1[x]*=k1;
	f1[x]%=MOD;
	f2[x]=(k1*f2[x]+k2)%MOD;
}
inline void spread(int x,int l,int r)
{
	int mid=(l+r)/2;
	change(x*2,l,mid,f1[x],f2[x]);
	change(x*2+1,mid+1,r,f1[x],f2[x]);
	f1[x]=1;
	f2[x]=0;
}
inline void update1(int t,int x,int y,int l,int r,int k)
{
	if(l<=x&&y<=r)
	{
		cnt[t]*=k;cnt[t]%=MOD;
		f1[t]*=k;f1[t]%=MOD;
		f2[t]*=k;f2[t]%=MOD;
		return;
	}
	spread(t,x,y);
	int mid=(x+y)/2;
	if(l<=mid) update1(t*2,x,mid,l,r,k);
	if(mid<r) update1(t*2+1,mid+1,y,l,r,k);
	cnt[t]=cnt[t*2]+cnt[t*2+1];
} 
inline void update2(int t,int x,int y,int l,int r,int k)
{
	if(l<=x&&y<=r)
	{
		cnt[t]=cnt[t]+k*(y-x+1)%MOD;
		cnt[t]%=MOD;
		f2[t]+=k;
		f2[t]%=MOD;
		return;
	} 
	spread(t,x,y);
	int mid=(x+y)/2;
	if(l<=mid) update2(t*2,x,mid,l,r,k);
	if(mid<r) update2(t*2+1,mid+1,y,l,r,k);
	cnt[t]=cnt[t*2]+cnt[t*2+1]; 
}
inline int ask(int t,int x,int y,int l,int r)
{
	if(l<=x&&y<=r) return cnt[t];
	spread(t,x,y);
	int res=0,mid=(x+y)/2;
	if(l<=mid)
	{
		res+=ask(t*2,x,mid,l,r);
		res%=MOD;
	}
	if(mid<r)
	{
		res+=ask(t*2+1,mid+1,y,l,r);
		res%=MOD;
	}
	return res%MOD;
}
signed main()
{
	cin>>n>>m>>MOD;
	for(int i=1;i<=n;i++) cin>>a[i],a[i]%=MOD;
	build(1,1,n);
	for(int i=1;i<=m;i++)
	{
		int cnt,l,r,k;
		cin>>cnt;
		if(cnt==1)
		{
			cin>>l>>r>>k;
			update1(1,1,n,l,r,k);
		}
		else if(cnt==2)
		{
			cin>>l>>r>>k;
			update2(1,1,n,l,r,k);
		}
		else
		{
			cin>>l>>r;
			cout<<ask(1,1,n,l,r)<<endl; 
		}
	}
	return 0;
}