值得一做》关于双标记线段树两三事BZOJ 1798 （NORMAL-）

　　这是一道双标记线段树的题，很让人很好的预习/学习/复习线段树，我不知道它能让别人学习什么，反正让我对线段树的了解更加深刻。

　　题目没什么好讲的，程序也没什么好讲的，所以也没有什么题解，但是值得一做

　　给出题目&代码

Description

老师交给小可可一个维护数列的任务，现在小可可希望你来帮他完成。 有长为N的数列，不妨设为a1,a2,…,aN 。有如下三种操作形式： (1)把数列中的一段数全部乘一个值; (2)把数列中的一段数全部加一个值; (3)询问数列中的一段数的和，由于答案可能很大，你只需输出这个数模P的值。

Input

第一行两个整数N和P(1≤P≤1000000000）。第二行含有N个非负整数,从左到右依次为a1,a2,…,aN, (0≤ai≤1000000000,1≤i≤N)。第三行有一个整数M，表示操作总数。从第四行开始每行描述一个操作，输入的操作有以下三种形式： 操作1：“1 t g c”(不含双引号)。表示把所有满足t≤i≤g的ai改为ai×c (1≤t≤g≤N,0≤c≤1000000000)。 操作2：“2 t g c”(不含双引号)。表示把所有满足t≤i≤g的ai改为ai+c (1≤t≤g≤N,0≤c≤1000000000)。 操作3：“3 t g”(不含双引号)。询问所有满足t≤i≤g的ai的和模P的值 (1≤t≤g≤N)。 同一行相邻两数之间用一个空格隔开，每行开头和末尾没有多余空格。

Output

对每个操作3，按照它在输入中出现的顺序，依次输出一行一个整数表示询问结果。

Sample Input

7 43
1 2 3 4 5 6 7
5
1 2 5 5
3 2 4
2 3 7 9
3 1 3
3 4 7

Sample Output

2
35
8

HINT

【样例说明】

初始时数列为(1,2,3,4,5,6,7)。
经过第1次操作后，数列为(1,10,15,20,25,6,7)。
对第2次操作，和为10+15+20=45，模43的结果是2。
经过第3次操作后，数列为(1,10,24,29,34,15,16}
对第4次操作，和为1+10+24=35，模43的结果是35。
对第5次操作，和为29+34+15+16=94,模43的结果是8。

测试数据规模如下表所示

数据编号 1 2 3 4 5 6 7 8 9 10
N= 10 1000 1000 10000 60000 70000 80000 90000 100000 100000
M= 10 1000 1000 10000 60000 70000 80000 90000 100000 100000

 /**************************************************************
     Problem: 1798
     User: PencilWang
     Language: C++
     Result: Accepted
     Time:4320 ms
     Memory:27388 kb
 ****************************************************************/

 #include<stdio.h>
 long long n,m;
 long long MOD;
 struct shit{
 int L,R;
 long long num,j,c;
 }s[];
 long long w[];
 long long C,J;
 void fuck(int p)
 {
     s[p].num=(s[p<<].num+s[p<<|].num)%MOD;
     return ;
 }
 void suck(int p)
 {
     int mid=(s[p].L+s[p].R)>>;
     int LL=p<<,RR=p<<|;
     s[LL].num=(s[LL].num*s[p].c+(s[LL].R-s[LL].L+)*s[p].j)%MOD;
     s[LL].c=(s[LL].c*s[p].c)%MOD;
     s[LL].j=(s[LL].j*s[p].c+s[p].j)%MOD;
     s[RR].num=(s[RR].num*s[p].c+(s[RR].R-s[RR].L+)*s[p].j)%MOD;
     s[RR].c=(s[RR].c*s[p].c)%MOD;
     s[RR].j=(s[RR].j*s[p].c+s[p].j)%MOD;
     s[p].c=,s[p].j=;
     return ;
 }
 void build(int p,int l,int r)
 {
     s[p].L=l;
     s[p].R=r;
     s[p].c=;
     s[p].j=;
     if(l==r)
     {
     s[p].num=w[l];
     return ;
     }
     int mid=(l+r)>>;
     build(p<<,l,mid);
     build(p<<|,mid+,r);
     fuck(p);
     return ;
 }
 void cc(int a,int b,int x,int p)
 {
     if(a<=s[p].L&&s[p].R<=b)
     {
     s[p].num=(x*s[p].num)%MOD;
     s[p].c=(x*s[p].c)%MOD;
     s[p].j=(x*s[p].j)%MOD;
     return ;
     }
     suck(p);
     int mid=(s[p].L+s[p].R)>>;
     if(mid>=a)cc(a,b,x,p<<);
     if(mid<b)cc(a,b,x,p<<|);
     fuck(p);
     return ;
 }
 void jj(int a,int b,int x,int p)
 {
     if(a<=s[p].L&&s[p].R<=b)
     {
     s[p].num=(s[p].num+x*(s[p].R-s[p].L+))%MOD;
     s[p].j=(x+s[p].j)%MOD;
     return ;
     }
     suck(p);
     int mid=(s[p].L+s[p].R)>>;
     if(mid>=a)jj(a,b,x,p<<);
     if(mid<b)jj(a,b,x,p<<|);
     fuck(p);
     return ;
 }
 long long Q(int p,int a,int b)
 {
     if(a<=s[p].L&&s[p].R<=b)
     return s[p].num%MOD;
     int mid=(s[p].L+s[p].R)>>;
     suck(p);
     long long ANS=;
     if(a<=mid)
     {
     ANS+=Q(p<<,a,b);
     ANS%=MOD;
     }
     if(b>mid)
     {
     ANS+=Q(p<<|,a,b);
     ANS%=MOD;
     }
     fuck(p);
     return ANS;
 }
 int main()
 {
     long long a,b,c,f;
     scanf("%lld%lld",&n,&MOD);
     for(int i=;i<=n;i++)scanf("%lld",w+i);
     build(,,n);
     scanf("%lld",&m);
     while(m--)
     {
         scanf("%lld%lld%lld",&f,&a,&b);
         if(f==)
         {
         scanf("%lld",&c);
         c%=MOD;
         cc(a,b,c,);
         }
         else if(f==)
         {
         scanf("%lld",&c);
         c%=MOD;
         jj(a,b,c,);
         }
         else
         printf("%lld\n",Q(,a,b)%MOD);
     }
     return ;
 }

1798