P3373 【模板】线段树 2 区间求和 区间乘 区间加

题目描述

如题，已知一个数列，你需要进行下面两种操作：

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的结果。

输入输出样例

输入样例#1：

5 5 38
1 5 4 2 3
2 1 4 1
3 2 5
1 2 4 2
2 3 5 5
3 1 4

输出样例#1：

17
2

说明

时空限制：1000ms,128M

数据规模：

对于30%的数据：N<=8，M<=10

对于70%的数据：N<=1000，M<=10000

对于100%的数据：N<=100000，M<=100000

（数据已经过加强^_^）

样例说明：

故输出应为17、2（40 mod 38=2）

根据加减法原理，，

好像只能这么解释，

先放乘法标记

再放加法标记

注意查询的时候ll和rr是不变的

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<cmath>
 #define LLI long long
 using namespace std;
 const LLI MAXN=;
 LLI read(LLI & n)
 {
     char p='+';LLI x=;
     while(p<''||p>'')
         p=getchar();
     while(p>=''&&p<='')
     x=x*+p-,p=getchar();
     n=x;
 }
 LLI n,m,mod,wl,wr,wv,ans;
 struct node
 {
     LLI l,r,w,fc,fj;
 }a[MAXN];
 void update(LLI k)
 {
     a[k].w=(a[k<<].w+a[k<<|].w)%mod;
 }
 void build_tree(LLI k,LLI ll,LLI rr)
 {
     a[k].l=ll;a[k].r=rr;
     a[k].fc=;
     a[k].fj=;
     if(a[k].l==a[k].r)
     {
         read(a[k].w);
         return ;
     }
     LLI mid=(ll+rr)/;
     build_tree(k<<,ll,mid);
     build_tree(k<<|,mid+,rr);
     update(k);
 }
 void pushdown(LLI k,LLI ll,LLI rr,LLI mid)
 {
     a[k<<].w*=a[k].fc;a[k<<|].w*=a[k].fc;
     a[k<<].w+=a[k].fj*(mid-ll+);a[k<<|].w+=a[k].fj*(rr-mid);
     a[k<<].fc*=a[k].fc;a[k<<|].fc*=a[k].fc;
     a[k<<].fj*=a[k].fc;a[k<<|].fj*=a[k].fc;
     a[k<<].fj+=a[k].fj;a[k<<|].fj+=a[k].fj;
     a[k].fc=;a[k].fj=;
     a[k<<].w%=mod;a[k<<].fj%=mod;a[k<<].fc%=mod;
     a[k<<|].w%=mod;a[k<<|].fj%=mod;a[k<<|].fc%=mod;
 }
 void interval_add(LLI k,LLI ll,LLI rr,LLI v)
 {
     if(a[k].l>rr||a[k].r<ll)
         return ;
     if(ll<=a[k].l&&rr>=a[k].r)
     {
         a[k].w=(a[k].w+v*(a[k].r-a[k].l+))%mod;
         a[k].fj=(a[k].fj+v)%mod;
         return ;
     }
     LLI mid=(a[k].l+a[k].r)/;
     pushdown(k,a[k].l,a[k].r,mid);
     //if(ll<=mid)
     interval_add(k<<,ll,rr,v);
     //if(rr>mid)
     interval_add(k<<|,ll,rr,v);
     update(k);
 }
 void interval_mul(LLI k,LLI ll,LLI rr,LLI v)
 {
     if(a[k].l>rr||a[k].r<ll)
         return ;
     if(ll<=a[k].l&&rr>=a[k].r)
     {
         a[k].w*=v%mod;
         a[k].fc*=v%mod;
         a[k].fj*=v%mod;
         return ;
     }
     LLI mid=(a[k].l+a[k].r)/;
     pushdown(k,a[k].l,a[k].r,mid);
     //if(ll<=mid)
     interval_mul(k<<,ll,rr,v);
     //if(rr>mid)
     interval_mul(k<<|,ll,rr,v);
     update(k);
 }
 void interval_sum(LLI k,LLI ll,LLI rr)
 {
     if(a[k].l>rr||a[k].r<ll)
         return ;
     if(ll<=a[k].l&&rr>=a[k].r)
     {
         ans=(ans+a[k].w)%mod;
         return ;
     }
     LLI mid=(a[k].l+a[k].r)/;
     pushdown(k,a[k].l,a[k].r,mid);
     //if(ll<=mid)
         interval_sum(k<<,ll,rr);
     //if(rr>mid)
         interval_sum(k<<|,ll,rr);
 }
 int main()
 {
     read(n);read(m);read(mod);
     build_tree(,,n);
     for(LLI i=;i<=m;i++)
     {
         LLI p;
         read(p);
         if(p==)
         {
             read(wl);read(wr);read(wv);
             interval_mul(,wl,wr,wv);
         }
         else if(p==)
         {
             read(wl);read(wr);read(wv);
             interval_add(,wl,wr,wv);
         }
         else if(p==)
         {
             ans=;
             read(wl);read(wr);
             interval_sum(,wl,wr);
             //cout<<ans%mod<<endl;
             printf("%lld\n",ans%mod);
         }
     }
     return ;
 }