【BZOJ】【1798】【AHOI2009】Seq维护序列

线段树

　　属于线段树中级应用吧……

　　要打两种标记：乘法和加法标记。一开始我想着可以像只有加法标记那样，永不下传，查询的时候依次累加就好了。后来发现不会写……只好每次update的时候……遇到标记！下传！query的时候遇到标记！下传！暴力地来搞……

　　然后说下下传的细节：先传乘法，后传加法。因为传乘法标记的时候要连sum带add都要一起乘，如果先传add就会多乘……

 /**************************************************************
     Problem: 1798
     User: Tunix
     Language: C++
     Result: Accepted
     Time:3868 ms
     Memory:6748 kb
 ****************************************************************/

 //BZOJ 1798
 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<iostream>
 #include<algorithm>
 #define rep(i,n) for(int i=0;i<n;++i)
 #define F(i,j,n) for(int i=j;i<=n;++i)
 #define D(i,j,n) for(int i=j;i>=n;--i)
 //using namespace std;
 const int N=;
 typedef long long LL;
 //#define debug
 int n,m,P,a[N];
 int mul[N<<],add[N<<],sum[N<<];

 void read(int &v){
     v=; int sign=; char ch=getchar();
     while(ch<'' || ch>''){ if (ch=='-') sign=-; ch=getchar();}
     while(ch>='' && ch<=''){ v=v*+ch-''; ch=getchar();}
     v*=sign;
 }
 #define mid (l+r>>1)
 #define L (o<<1)
 #define R (o<<1|1)

 void Push_down(int o,int l,int r){
     if (mul[o]!=){
         add[L]=((LL)add[L]*mul[o])%P; add[R]=((LL)add[R]*mul[o])%P;
         mul[L]=((LL)mul[L]*mul[o])%P; mul[R]=((LL)mul[R]*mul[o])%P;
         mul[o]=;
     }
     if (add[o]!=){
         add[L]=((LL)add[L]+add[o])%P; add[R]=((LL)add[R]+add[o])%P;
         add[o]=;
     }
 }
 void Push_up(int o,int l,int r){
     if (l<r) sum[o]=(sum[L]+sum[R])%P; else sum[o]=;
     sum[o]=((LL)sum[o]*mul[o]+(LL)add[o]*(r-l+))%P;
 }
 void build(int o,int l,int r){
     if (l==r) add[o]=sum[o]=a[l],mul[o]=;
     else{
         build(L,l,mid);
         build(R,mid+,r);
         add[o]=; mul[o]=;
         Push_up(o,l,r);
     }
 }
 int ql,qr;
 void update_add(int o,int l,int r,int v){
     if (ql<=l && qr>=r)
         add[o]=((LL)add[o]+v)%P;
     else{
         if (mul[o]!= || add[o]) Push_down(o,l,r);
         if (ql<=mid) update_add(L,l,mid,v); else Push_up(L,l,mid);
         if (qr>mid) update_add(R,mid+,r,v); else Push_up(R,mid+,r);
     }
     Push_up(o,l,r);
 }
 void update_mul(int o,int l,int r,int v){
     if (ql<=l && qr>=r)
         mul[o]=((LL)mul[o]*v)%P,add[o]=(LL)add[o]*v%P;
     else{
         if (mul[o]!= || add[o]) Push_down(o,l,r);
         if (ql<=mid) update_mul(L,l,mid,v); else Push_up(L,l,mid);
         if (qr>mid) update_mul(R,mid+,r,v); else Push_up(R,mid+,r);
     }
     Push_up(o,l,r);
 }
 int sumv;
 void query(int o,int l,int r){
     if (ql<=l && qr>=r)
         sumv=((LL)sumv+sum[o])%P;
     else{
         if (mul[o]!= || add[o]) {Push_down(o,l,r); Push_up(L,l,mid); Push_up(R,mid+,r);}
         if (ql<=mid) query(L,l,mid);
         if (qr>mid) query(R,mid+,r);
     }
 }
 int ans[N],cnt=;
 int main(){
     #ifndef ONLINE_JUDGE
     freopen("input.txt","r",stdin);
 //  freopen("output.txt","w",stdout);
     #endif
     read(n); read(P);
     F(i,,n) read(a[i]);
 //  F(i,1,n) printf("%d ",a[i]);
 //  printf("\n");
     build(,,n);

     read(m);
 //  printf("m=%d\n",m);
     int c,v;
     F(i,,m){
         read(c); read(ql); read(qr);
 //      printf("c=%d\n",c);
         switch(c){
             case : read(v); update_mul(,,n,v); break;
             case : read(v); update_add(,,n,v); break;
             case : sumv=; query(,,n); ans[cnt++]=sumv;break;
         }
     }
     rep(i,cnt-) printf("%d\n",ans[i]);
     printf("%d",ans[cnt-]);
     return ;
 }