[BZOJ1798][AHOI2009]Seq维护序列 线段树

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1798

一眼看过去线段树，事实上就是线段树。对于乘和加的两个标记，我们可以规定一个顺序，比如先乘后加。每次乘的时候利用乘法分配律就可以分别更新标记，而加就直接加上去就好了。

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 typedef long long ll;
 int inline readint(){
     int Num;char ch;
     while((ch=getchar())<''||ch>'');Num=ch-'';
     while((ch=getchar())>=''&&ch<='') Num=Num*+ch-'';
     return Num;
 }
 void outint(int x){
     if(x>=) outint(x/);
     putchar(x%+'');
 }
 int N,M,P;
 int sum[(<<)+],col1[(<<)+],col2[(<<)+];
 #define lson l,mid,rt<<1
 #define rson mid+1,r,rt<<1|1
 void Pushup(int &rt){
     sum[rt]=(sum[rt<<]+sum[rt<<|])%P;
 }
 void Build(int l,int r,int rt){
     col1[rt]=;
     if(l==r){
         sum[rt]=readint()%P;
         return;
     }
     int mid=l+r>>;
     if(l<=mid) Build(lson);
     if(r>mid) Build(rson);
     Pushup(rt);
 }
 void Pushdown(int &rt,int &l,int &r){
     if(col1[rt]!=){
         col1[rt<<]=(ll)col1[rt<<]*col1[rt]%P;
         col2[rt<<]=(ll)col2[rt<<]*col1[rt]%P;
         sum[rt<<]=(ll)sum[rt<<]*col1[rt]%P;
         col1[rt<<|]=(ll)col1[rt<<|]*col1[rt]%P;
         col2[rt<<|]=(ll)col2[rt<<|]*col1[rt]%P;
         sum[rt<<|]=(ll)sum[rt<<|]*col1[rt]%P;
         col1[rt]=;
     }
     if(col2[rt]){
         col2[rt<<]=(col2[rt<<]+col2[rt])%P;
         sum[rt<<]=((ll)col2[rt]*((l+r>>)-l+)%P+sum[rt<<])%P;
         col2[rt<<|]=(col2[rt<<|]+col2[rt])%P;
         sum[rt<<|]=((ll)col2[rt]*(r-(l+r>>))%P+sum[rt<<|])%P;
         col2[rt]=;
     }
 }
 void Mul(int l,int r,int rt,int L,int R,int x){
     if(l>=L&&r<=R){
         col1[rt]=(ll)col1[rt]*x%P;
         col2[rt]=(ll)col2[rt]*x%P;
         sum[rt]=(ll)sum[rt]*x%P;
         return;
     }
     Pushdown(rt,l,r);
     int mid=l+r>>;
     if(L<=mid) Mul(lson,L,R,x);
     if(R>mid) Mul(rson,L,R,x);
     Pushup(rt);
 }
 void Add(int l,int r,int rt,int L,int R,int x){
     if(l>=L&&r<=R){
         col2[rt]=(col2[rt]+x)%P;
         sum[rt]=((ll)(r-l+)*x%P+sum[rt])%P;
         return;
     }
     Pushdown(rt,l,r);
     int mid=l+r>>;
     if(L<=mid) Add(lson,L,R,x);
     if(R>mid) Add(rson,L,R,x);
     Pushup(rt);
 }
 int Qry(int l,int r,int rt,int L,int R){
     if(l>=L&&r<=R) return sum[rt];
     Pushdown(rt,l,r);
     int mid=l+r>>,ret=;
     if(L<=mid) ret=Qry(lson,L,R);
     if(R>mid) ret=(ret+Qry(rson,L,R))%P;
     return ret;
 }
 int main(){
     N=readint();
     P=readint();
     Build(,N,);
     M=readint();
     for(int i=;i<=M;i++){
         int opt=readint(),t,g,c;
         switch(opt){
             case :
                 t=readint();
                 g=readint();
                 c=readint()%P;
                 Mul(,N,,t,g,c);
                 break;
             case :
                 t=readint();
                 g=readint();
                 c=readint()%P;
                 Add(,N,,t,g,c);
                 break;
             case :
                 t=readint();
                 g=readint();
                 outint(Qry(,N,,t,g));
                 putchar('\n');
                 break;
         }
     }
     return ;
 }