rt,只有124对了,其他全WA
#include<bits/stdc++.h>
#define ll long long
#define lid id<<1
#define rid (lid)+1
using namespace std;
const int N=100005;
ll MOD,k,a[N];
int n,m,op,x,y;
struct tree
{
int l,r;
ll sum,add,mul;
}tr[N<<2];
void pu(int id)
{
tr[id].sum=(tr[lid].sum*tr[lid].mul+tr[lid].add*(tr[lid].r-tr[lid].l+1)+\
tr[rid].sum*tr[rid].mul+tr[rid].add*(tr[rid].r-tr[rid].l+1))%MOD;
}
void pd(int id)
{
ll &add=tr[id].add,&mul=tr[id].mul;
tr[id].sum=(tr[id].sum*mul+add)%MOD;
tr[lid].add=(tr[lid].add*mul+add)%MOD,tr[lid].mul=tr[lid].mul*mul%MOD;
tr[rid].add=(tr[rid].add*mul+add)%MOD,tr[rid].mul=tr[rid].mul*mul%MOD;
add=0,mul=1;
}
void build(int l,int r,int id)
{
tr[id].l=l,tr[id].r=r,tr[id].add=0,tr[id].mul=1;
if(l==r)
{
tr[id].sum=a[l];
return;
}
int mid=l+r>>1;
build(l,mid,lid);
build(mid+1,r,rid);
pu(id);
}
void add(int l,int r,int id)
{
if(tr[id].l==l&&tr[id].r==r)
{
tr[id].add=(tr[id].add+k)%MOD;
return;
}
pd(id);
if(tr[lid].r>=l)
{
if(tr[rid].l<=r)
{
add(l,tr[lid].r,lid);
add(tr[rid].l,r,rid);
}
else add(l,r,lid);
}
else add(l,r,rid);
pu(id);
}
void mul(int l,int r,int id)
{
if(tr[id].l==l&&tr[id].r==r)
{
tr[id].add=tr[id].add*k%MOD,tr[id].mul=tr[id].mul*k%MOD;
return;
}
pd(id);
if(tr[lid].r>=l)
{
if(tr[rid].l<=r)
{
mul(l,tr[lid].r,lid);
mul(tr[rid].l,r,rid);
}
else mul(l,r,lid);
}
else mul(l,r,rid);
pu(id);
}
ll query(int l,int r,int id)
{
if(tr[id].l==l&&tr[id].r==r)
return (tr[id].sum*tr[id].mul+tr[id].add*(r-l+1))%MOD;
pd(id);
if(tr[lid].r>=l)
if(tr[rid].l<=r)
return (query(l,tr[lid].r,lid)+query(tr[rid].l,r,rid))%MOD;
else return query(l,r,lid);
else query(l,r,rid);
}
int main()
{
cin>>n>>m>>MOD;
for(int i=1;i<=n;i++)scanf("%lld",&a[i]);
build(1,n,1);
while(m--)
{
scanf("%d%d%d",&op,&x,&y);
if(op==3)printf("%lld\n",query(x,y,1)%MOD);
else
{
scanf("%lld",&k);
if(op==1)mul(x,y,1);
else add(x,y,1);
}
}
return 0;
}