题目描述
如题,已知一个数列,你需要进行下面三种操作:
1.将某区间每一个数乘上x
2.将某区间每一个数加上x
3.求出某区间每一个数的和
输入输出格式
输入格式:第一行包含三个整数N、M、P,分别表示该数列数字的个数、操作的总个数和模数。
第二行包含N个用空格分隔的整数,其中第i个数字表示数列第i项的初始值。
接下来M行每行包含3或4个整数,表示一个操作,具体如下:
操作1: 格式:1 x y k 含义:将区间[x,y]内每个数乘上k
操作2: 格式:2 x y k 含义:将区间[x,y]内每个数加上k
操作3: 格式:3 x y 含义:输出区间[x,y]内每个数的和对P取模所得的结果
输出格式:输出包含若干行整数,即为所有操作3的结果。
输入输出样例
说明
时空限制:1000ms,128M
数据规模:
对于30%的数据:N<=8,M<=10
对于70%的数据:N<=1000,M<=10000
对于100%的数据:N<=100000,M<=100000
(数据已经过加强^_^)
样例说明:
故输出应为17、2(40 mod 38=2)
就是基本的线段树操作,但是有几个地方要多注意下,也就是我犯错的地方。
1 乘法优先级高于加法
2 初始tag数组时不能用memset把mul全变成1,要在build中赋初始值
3 long long记得最后要改好
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define ll long long
#define il inline
#define db double
using namespace std;
il int gi()
{
int x=0,y=1;
char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-')
y=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
x=x*10+ch-'0';
ch=getchar();
}
return x*y;
}
il ll gl()
{
ll x=0,y=1;
char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-')
y=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
x=x*10+ch-'0';
ch=getchar();
}
return x*y;
}
int n,mm;
ll p;
ll a[100045],sum[400045];
ll add[400045],mul[400045];
il void build(int rt,int l,int r)
{
mul[rt]=1;
add[rt]=0;
int m=(l+r)>>1;
if(l==r)
{
sum[rt]=a[l];
return;
}
build(rt<<1,l,m);
build(rt<<1|1,m+1,r);
sum[rt]=(sum[rt<<1]+sum[rt<<1|1])%p;
}
il void pushdown(int rt,int l,int r)
{
int m=(l+r)>>1;
sum[rt<<1]=(sum[rt<<1]*mul[rt]%p+add[rt]*(m-l+1)%p)%p;
sum[rt<<1|1]=(sum[rt<<1|1]*mul[rt]%p+add[rt]*(r-m)%p)%p;
mul[rt<<1]=mul[rt<<1]*mul[rt]%p;
mul[rt<<1|1]=mul[rt<<1|1]*mul[rt]%p;
add[rt<<1]=(add[rt<<1]*mul[rt]%p+add[rt])%p;
add[rt<<1|1]=(add[rt<<1|1]*mul[rt]%p+add[rt])%p;
mul[rt]=1;
add[rt]=0;
}
il void update1(int rt,int l,int r,int L,int R,ll k)
{
if(L<=l&&R>=r)
{
sum[rt]=sum[rt]*k%p;
mul[rt]=mul[rt]*k%p;
add[rt]=add[rt]*k%p;
return;
}
pushdown(rt,l,r);
int m=(l+r)>>1;
if(L<=m)
update1(rt<<1,l,m,L,R,k);
if(R>m)
update1(rt<<1|1,m+1,r,L,R,k);
sum[rt]=(sum[rt<<1]+sum[rt<<1|1])%p;
}
il void update2(int rt,int l,int r,int L,int R,ll k)
{
if(L<=l&&R>=r)
{
sum[rt]=(sum[rt]+k*(r-l+1))%p;
add[rt]=(k+add[rt])%p;
return;
}
pushdown(rt,l,r);
int m=(l+r)>>1;
if(L<=m)
update2(rt<<1,l,m,L,R,k);
if(R>m)
update2(rt<<1|1,m+1,r,L,R,k);
sum[rt]=(sum[rt<<1]+sum[rt<<1|1])%p;
}
il ll query(int rt,int l,int r,int L,int R)
{
if(L<=l&&R>=r)
return sum[rt];
pushdown(rt,l,r);
int m=(l+r)>>1;
ll s=0;
if(L<=m)
s=(s+query(rt<<1,l,m,L,R))%p;
if(R>m)
s=(s+query(rt<<1|1,m+1,r,L,R))%p;
return s%p;
}
int main()
{
n=gi(),mm=gi(),p=gl();
for(int i=1;i<=n;i++)
a[i]=gl();
build(1,1,n);
int x,y,k,r;
for(int i=1;i<=mm;i++)
{
r=gi();
if(r==1)
{
x=gi(),y=gi(),k=gl();
update1(1,1,n,x,y,k);
}
if(r==2)
{
x=gi(),y=gi(),k=gl();
update2(1,1,n,x,y,k);
}
if(r==3)
{
x=gi(),y=gi();
printf("%lld\n",query(1,1,n,x,y));
}
}
return 0;
}