Loj 6283. 数列分块入门 7

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int M = 1e5+10;
int mtag[M],atag[M],a[M],bl[M],n,block;
const int mod = 10007;
void resert(int x){
    for(int i = (x-1)*block+1;i <= min(x*block,n);i ++)
        a[i] = (a[i]*mtag[x] + atag[x])%mod;
    atag[x] = 0; mtag[x] = 1;
}

void update(int f,int l,int r,int c){
     resert(bl[l]);
     for(int i = l;i <= min(bl[l]*block,r);i ++){
        if(f == 1) a[i] = (a[i]*c)%mod;
        else a[i] = (a[i]+c)%mod;
     }
     if(bl[l] != bl[r]){
        resert(bl[r]);
        for(int i = (bl[r]-1)*block+1;i <= r;i ++){
            if(f == 1) a[i] = (a[i]*c)%mod;
            else a[i] = (a[i] + c)%mod;
         }
     }
     for(int i = bl[l]+1;i <= bl[r]-1;i ++){
        if(f == 1) {
            atag[i] = (atag[i]*c)%mod; mtag[i] = (mtag[i]*c)%mod;
        }
        else atag[i] = (atag[i]+c)%mod;
     }
}

int main()
{
    int f,l,r,c;
    scanf("%d",&n);
    block = sqrt(n);
    for(int i = 1;i <= n;i ++)  scanf("%d",&a[i]);
    for(int i = 1;i <= n;i ++)  bl[i] = (i-1)/block + 1;
    for(int i = 1;i <= bl[n];i ++)  mtag[i] = 1,atag[i] = 0;
    for(int i = 1;i <= n;i ++){
        scanf("%d%d%d%d",&f,&l,&r,&c);
        if(f == 2) printf("%d\n",(a[r]*mtag[bl[r]] + atag[bl[r]])%mod);
        else update(f,l,r,c);
    }
    return 0;
}