算法模板——线段树2（区间加+区间乘+区间求和）

该模板实现的功能——进行区间的乘法和加法，以及区间的求和（1：乘法 2：加法 3：求和）详见BZOJ1798

type
    vet=record
              a0,a1:int64;
    end;
var
   i,j,k,l,m,n,a2,a3,a4:longint;
   p:int64;
   d1,d2,d:vet;
   a,c:array[0..1000000] of int64;
   b:array[0..1000000] of vet;
function min(x,y:longint):longint;inline;
         begin
              if x<y then min:=x else min:=y;
         end;
function max(x,y:longint):longint;inline;
         begin
              if x>y then max:=x else max:=y;
         end;
function merge(d1,d2:vet):vet;inline;
         var d3:vet;
         begin
              d3:=d1;
              d3.a0:=d3.a0 mod p;
              d3.a1:=d3.a1 mod p;
              d2.a0:=d2.a0 mod p;
              d2.a1:=d2.a1 mod p;
              d3.a0:=(d3.a0*d2.a0) mod p;
              d3.a1:=((d3.a1*d2.a0) mod p+d2.a1) mod p;
              exit(d3);
         end;
procedure built(z,x,y:longint);inline;
          begin
               if (x=y) then
                  a[z]:=c[x] mod p
               else
                   begin
                        built(z*2,x,(x+y) div 2);
                        built(z*2+1,(x+y) div 2+1,y);
                        a[z]:=(a[z*2]+a[z*2+1]) mod p;
                   end;
               b[z].a0:=1;b[z].a1:=0;
          end;
procedure ext(z,x,y:longint);inline;
          begin
               a[z]:=((a[z]*b[z].a0) mod p+b[z].a1*(y-x+1)) mod p;
               b[z*2]:=merge(b[z*2],b[z]);
               b[z*2+1]:=merge(b[z*2+1],b[z]);
               b[z].a0:=1;b[z].a1:=0;
          end;
function op(z,x,y,l,r:longint;d:vet):int64;inline;
          var a2,a3,a4:int64;
          begin
               if l>r then exit(0);
               ext(z,x,y);a2:=a[z];
               if (x=l) and (y=r) then
                  begin
                       b[z]:=d;
                       exit((a2*(b[z].a0-1)) mod p+(b[z].a1*(r-l+1)) mod p);
                  end
               else
                   begin
                        a3:=op(z*2,x,(x+y) div 2,l,min((x+y) div 2,r),d) mod p;
                        a4:=op(z*2+1,(x+y) div 2+1,y,max((x+y) div 2+1,l),r,d) mod p;
                        a[z]:=(a[z]+(a3+a4) mod p) mod p;exit((a3+a4) mod p);
                   end;
          end;
function sum(z,x,y,l,r:longint;d:vet):int64;inline;
         var d1,d2:vet;
         begin
              if l>r then exit(0);
              d1:=b[z];d1:=merge(d1,d);;
              if (x=l) and (y=r) then
                 exit(((d1.a0*a[z]) mod p+(d1.a1*(r-l+1)) mod p) mod p)
              else
                  exit((sum(z*2,x,(x+y) div 2,l,min((x+y) div 2,r),d1)+sum(z*2+1,(x+y) div 2+1,y,max((x+y) div 2+1,l),r,d1)) mod p);
         end;
procedure showoff(z,x,y,l:longint);inline;
          begin
               writeln('':l*2,z,'(',x,',',y,') = Tag=(',b[z].a0,',',b[z].a1,') ',a[z]);
               if x<y then
                  begin
                       showoff(z*2,x,(x+y) div 2,l+1);
                       showoff(z*2+1,(x+y) div 2+1,y,l+1);
                  end;
          end;
begin
     readln(n,p);
     for i:=1 to n do read(c[i]);
     readln;
     built(1,1,n);
     readln(m);
     for i:=1 to m do
         begin
              read(j);
              case j of
                   1:begin
                          readln(a2,a3,a4);
                          d1.a0:=a4;d1.a1:=0;
                          op(1,1,n,a2,a3,d1);
                   end;
                   2:begin
                          readln(a2,a3,a4);
                          d1.a0:=1;d1.a1:=a4;
                          op(1,1,n,a2,a3,d1);
                   end;
                   3:begin
                          readln(a2,a3);
                          d1.a0:=1;d1.a1:=0;
                          writeln((sum(1,1,n,a2,a3,d1)+p) mod p);
                   end;
              end;
         end;
     readln;
end.