#include<bits/stdc++.h>
using namespace std;
long long t[500000],a[500000],lan[500000],p;
void fl(long long x,long long l,long long r)
{
	if(lan[x]==0) return ;
	long long val=lan[x]%p,mid=(l+r)/2;
	lan[x*2]+=val%p,t[x*2]+=(mid-l+1)*val%p;
	lan[x*2+1]+=val%p,t[x*2+1]+=(r-mid)*val%p;
	lan[x]=0;
}
void f(long long x,long long l,long long r)
{
		if(l==r)
		{
			t[x]=a[l];
			return ;
		}
		long long mid=(l+r)/2;
		f(x*2,l,mid);
		f(x*2+1,mid+1,r);
		t[x]=(t[x*2]%p+t[x*2+1]%p)%p;
}
long long sum(long long x,long long l,long long r,long long L,long long R)
{
	if(L<=l&&r<=R)
	return t[x];
	long long ans=0;
	long long mid=(l+r)/2;
	fl(x,l,r);
	if(L<=mid)
	ans+=sum(x*2,l,mid,L,R)%p;
	if(mid<R)
	ans+=sum(x*2+1,mid+1,r,L,R)%p;
	t[x]=(t[x*2]%p+t[x*2+1]%p)%p;
	return ans%p;
}

void f3(long long x,long long l,long long r,long long L,long long R,long long v)
{
	if(L<=l&&r<=R){
		lan[x]+=v%p;
		t[x]+=(r-l+1)*v%p;
		return ;
	}
	fl(x,l,r);
	long long mid=(l+r)/2;
	if(L<=mid) f3(x*2,l,mid,L,R,v);
	if(mid<R) f3(x*2+1,mid+1,r,L,R,v);
	t[x]=(t[x*2]%p+t[x*2+1]%p)%p;
}
void f4(long long x,long long l,long long r,long long L,long long R,long long v)
{
	if(L<=l&&r<=R){
		lan[x]*=v%p;
		t[x]*=(r-l+1)*v%p;
		return ;
	}
	fl(x,l,r);
	long long mid=(l+r)/2;
	if(L<=mid) f3(x*2,l,mid,L,R,v);
	if(mid<R) f3(x*2+1,mid+1,r,L,R,v);
	t[x]=(t[x*2]%p+t[x*2+1]%p)%p;
}
long long n,m;
int main()
{
	cin>>n>>m>>p;
	for(long long i=1;i<=n;i++)
	scanf("%lld",&a[i]);
	f(1,1,n);
	for(long long i=1;i<=m;i++)
	{
		long long p2;
		scanf("%lld",&p2);
		if(p2==1)
		{
			long long u,v,w;
			cin>>w>>u>>v;
			f4(1,1,n,w,u,v);
		
		}
			else if(p2==3)
			{
				long long u,v;
				cin>>u>>v;
				cout<<sum(1,1,n,u,v)<<endl;
			}
			else if(p2==2)
			{
				long long u,v,w;
			cin>>w>>u>>v;
			f3(1,1,n,w,u,v);
			}
		
	}
	return 0;
}