分块——入门7

如果只是加法操作的话就和入门1没啥区别了，现在最大的问题就是乘法咋搞？？
首先我们要知道的是，乘法的优先级是高于加法，那么我们在做乘法的时候就要十分的小心了，所以我们在修改乘法标记的时候要把加法标记也作修改。
例如：
若当前的一个块乘以m1后加上a1，这时进行一个乘m2的操作，则原来的标记变成m1*m2，a1*m2
若当前的一个块乘以m1后加上a1，这时进行一个加a2的操作，则原来的标记变成m1，a1+a2

代码：

//By Bibi
///                 .-~~~~~~~~~-._       _.-~~~~~~~~~-.
///             __.'              ~.   .~              `.__
///           .'//                  \./                  \\`.
///        .'//                     |                     \\`.
///       .'// .-~"""""""~~~~-._     |     _,-~~~~"""""""~-. \\`.
///     .'//.-"                 `-.  |  .-'                 "-.\\`.
///   .'//______.============-..   \ | /   ..-============.______\\`.
/// .'______________________________\|/______________________________`.
#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define dep(i,a,b) for(int i=a;i>=b;--i)
using namespace std;
const int MAXN=100005;
const int mod=10007;
int read()
{
    int sum=0,flag=1;
    char c;
    for(;c<'0'||c>'9';c=getchar())if(c=='-') flag=-1;
    for(;c>='0'&&c<='9';c=getchar())sum=(sum<<1)+(sum<<3)+c-'0';
    return sum*flag;
}
int n,block,num,l[MAXN],r[MAXN],L,R,c,op,a[MAXN],p[MAXN],be[MAXN],pp[MAXN];
void build()
{
    block=sqrt(n);
    num=n/block;
    if(n%block) num++;
    rep(i,1,n) be[i]=(i-1)/block+1;
    rep(i,1,num)
    {
        l[i]=(i-1)*block+1;
        r[i]=i*block;
    }
    r[num]=n;
}
void push(int ll,int rr)
{
    rep(i,ll,rr) a[i]=(a[i]*pp[be[ll]]+p[be[ll]])%mod;
    p[be[ll]]=0,pp[be[ll]]=1;
} 

void add()
{
    push(l[be[L]],r[be[L]]);
    if(be[L]==be[R]) rep(i,L,R) a[i]+=c,a[i]%=mod;
    else
    {
        rep(i,L,r[be[L]]) a[i]+=c,a[i]%=mod;
        push(l[be[R]],r[be[R]]);
        rep(i,l[be[R]],R) a[i]+=c,a[i]%=mod;
        rep(i,be[L]+1,be[R]-1) p[i]+=c,p[i]%=mod;
    } 
}
void cheng()
{
    push(l[be[L]],r[be[L]]);
    if(be[L]==be[R]) rep(i,L,R) a[i]*=c,a[i]%=mod;
    else
    {
        rep(i,L,r[be[L]]) a[i]*=c,a[i]%=mod;
        push(l[be[R]],r[be[R]]);
        rep(i,l[be[R]],R) a[i]*=c,a[i]%=mod;
        rep(i,be[L]+1,be[R]-1) pp[i]*=c,p[i]*=c,pp[i]%=mod,p[i]%=mod;
    }
}
void init()
{
    n=read();
    rep(i,1,n) a[i]=read(),pp[i]=1;
    build();
    rep(i,1,n)
    {
        op=read();L=read();R=read();c=read();
        if(!op) add();
        else if(op==1) cheng(); 
        else printf("%d\n",(a[R]*pp[be[R]]+p[be[R]])%mod);
    }
}
int main()
{
    init();
    return 0;
}