LOJ#6283. 数列分块入门 7

对于每个区间先乘在加，如果我修改的是部分的块，我就需要把现这个块的add和mul标记全部放下去，然后再更新。

#include<map>
#include<set>
#include<ctime>
#include<cmath>
#include<stack>
#include<queue>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define first fi
#define second se
#define lowbit(x) (x & (-x))

typedef unsigned long long int ull;
typedef long long int ll;
const double pi = 4.0*atan(1.0);
const int inf = 0x3f3f3f3f;
const int maxn = ;
const int maxm = ;
const int mod = ;
using namespace std;

int n, m, tol, T;
int block;
int a[maxn];
int add[maxn];
int mul[maxn];
int belong[maxn];

void init() {
    memset(a, , sizeof a);
    memset(add, , sizeof add);
    memset(mul, , sizeof mul);
    memset(belong, , sizeof belong);
}

int L(int x) {
    return (x-)*block + ;
}

int R(int x) {
    return min(n, x*block);
}

void update1(int l, int r, int c) {
    for(int i=L(belong[l]); i<=R(belong[l]); i++) {
        a[i] = a[i] * mul[belong[i]] % mod + add[belong[i]];
        a[i] %= mod;
    }
    add[belong[l]] = ;
    mul[belong[l]] = ;
    for(int i=l; i<=min(r, R(belong[l])); i++) {
        a[i] += c;
        a[i] %= mod;
    }
    if(belong[l] == belong[r])    return ;
    if(belong[l] != belong[r]) {
        for(int i=L(belong[r]); i<=R(belong[r]); i++) {
            a[i] = a[i] * mul[belong[i]] % mod + add[belong[i]];
            a[i] %= mod;
        }
        add[belong[r]] = ;
        mul[belong[r]] = ;
        for(int i=L(belong[r]); i<=r; i++) {
            a[i] += c;
            a[i] %= mod;
        }
    }
    for(int i=belong[l]+; i<belong[r]; i++) {
        add[i] += c;
        add[i] %= mod;
    }
}

void update2(int l, int r, int c) {
    for(int i=L(belong[l]); i<=R(belong[l]); i++) {
        a[i] = a[i] * mul[belong[i]] % mod + add[belong[i]];
        a[i] %= mod;
    }
    add[belong[l]] = ;
    mul[belong[l]] = ;
    for(int i=l; i<=min(r, R(belong[l])); i++) {
        a[i] *= c;
        a[i] %= mod;
    }
    if(belong[l] == belong[r])    return ;
    if(belong[l] != belong[r]) {
        for(int i=L(belong[r]); i<=R(belong[r]); i++) {
            a[i] = a[i] * mul[belong[i]] % mod + add[belong[i]];
            a[i] %= mod;
        }
        add[belong[r]] = ;
        mul[belong[r]] = ;
        for(int i=L(belong[r]); i<=r; i++) {
            a[i] *= c;
            a[i] %= mod;
        }
    }
    for(int i=belong[l]+; i<belong[r]; i++) {
        add[i] *= c;
        mul[i] *= c;
        add[i] %= mod;
        mul[i] %= mod;
    }
}

int main() {
    while(~scanf("%d", &n)) {
        init();
        block = sqrt(n);
        for(int i=; i<=n; i++) {
            scanf("%d", &a[i]);
            belong[i] = (i-) / block + ;
            mul[i] = ;
        }
        m = n;
        while(m--) {
            int op, l, r, c;
            scanf("%d%d%d%d", &op, &l, &r, &c);
            if(op == ) {
                update1(l, r, c%mod);
            } else if(op == ) {
                update2(l, r, c%mod);
            } else {
                int ans = a[r]*mul[belong[r]]%mod + add[belong[r]]%mod;
                printf("%d\n", ans%mod);
            }
           // for(int i=1; i<=n; i++)    printf("%d%c", a[i]*mul[belong[i]]%mod + add[belong[i]]%mod, i==n ? '\n' : ' ');
        }
    }
    return ;
}