「线段树」[AHOI2009]维护序列

双倍经验，还是蓝题，岂不美哉

题目描述

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

输入格式

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

输出格式

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

输入输出样例

输入

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

输出

2
35
8

说明/提示

【样例说明】

初始时数列为(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

Source: Ahoi 2009

题目题解

还是三个操作，线段树基本题

//#define fre yes

#include <cstdio>
#define int long long

const int N = 100005;
struct Node {
    int l, r;
    long long sum, addv, mul;
} tree[N << 2];
int MOD;

long long ans;

void build(int k, int l, int r) {
    tree[k].l = l; tree[k].r = r; tree[k].addv = 0; tree[k].mul = 1;
    if(l == r) {
        scanf("%lld", &tree[k].sum);
        tree[k].sum %= MOD;
        return ;
    }

    int mid = (l + r) >> 1;
    build(k * 2, l, mid);
    build(k * 2 + 1, mid + 1, r);
    tree[k].sum = (tree[k * 2].sum + tree[k * 2 + 1].sum) % MOD;
}

void update(int k) {
    tree[k * 2].addv = (tree[k].mul * tree[k * 2].addv + tree[k].addv) % MOD;
    tree[k * 2 + 1].addv = (tree[k].mul * tree[k * 2 + 1].addv + tree[k].addv) % MOD;
    tree[k * 2].mul = (tree[k].mul * tree[k * 2].mul) % MOD;
    tree[k * 2 + 1].mul = (tree[k].mul * tree[k * 2 + 1].mul) % MOD;
    tree[k * 2].sum = (tree[k * 2].sum * tree[k].mul + (tree[k * 2].r - tree[k * 2].l + 1) * tree[k].addv) % MOD;
    tree[k * 2 + 1].sum = (tree[k * 2 + 1].sum * tree[k].mul + (tree[k * 2 + 1].r - tree[k * 2 + 1].l + 1) * tree[k].addv) % MOD;
    tree[k].addv = 0;
    tree[k].mul = 1;
}

void mul_interval(int k, int l, int r, int x) {
    if(tree[k].l >= l && tree[k].r <= r) {
        tree[k].sum = (tree[k].sum * x) % MOD;
        tree[k].addv = (tree[k].addv * x) % MOD;
        tree[k].mul = (tree[k].mul * x) % MOD;
        return ;
    }

    if(tree[k].addv || tree[k].mul != 1) update(k);
    int mid = (tree[k].l + tree[k].r) >> 1;
    if(mid >= l) mul_interval(k * 2, l, r, x);
    if(mid < r) mul_interval(k * 2 + 1, l, r, x);
    tree[k].sum = (tree[k * 2].sum + tree[k * 2 + 1].sum) % MOD;
}

void add_interval(int k, int l, int r, int x) {
    if(tree[k].l >= l && tree[k].r <= r) {
        tree[k].sum = (tree[k].sum + (tree[k].r - tree[k].l + 1) * x) % MOD;
        tree[k].addv = (tree[k].addv + x) % MOD;
        return ;
    }

    if(tree[k].addv || tree[k].mul != 1) update(k);
    int mid = (tree[k].l + tree[k].r) >> 1;
    if(mid >= l) add_interval(k * 2, l, r, x);
    if(mid < r) add_interval(k * 2 + 1, l, r, x);
    tree[k].sum = (tree[k * 2].sum + tree[k * 2 + 1].sum) % MOD;
}

void query(int k, int l, int r) {
    if(tree[k].l >= l && tree[k].r <= r) {
        ans = (ans + tree[k].sum) % MOD;
        return ;
    }

    if(tree[k].addv || tree[k].mul != 1) update(k);
    int mid = (tree[k].l + tree[k].r) >> 1;
    if(mid >= l) query(k * 2, l, r);
    if(mid < r) query(k * 2 + 1, l, r);
}

signed main() {
    static int n, m;
    scanf("%lld %lld", &n, &MOD);
    build(1, 1, n);
    scanf("%lld", &m);
    for (int i = 1; i <= m; i++) {
        int k;
        scanf("%lld", &k);
        if(k == 1) {
            int l, r, x;
            scanf("%lld %lld %lld", &l, &r, &x);
            mul_interval(1, l, r, x);
        }
        if(k == 2) {
            int l, r, x;
            scanf("%lld %lld %lld", &l, &r, &x);
            add_interval(1, l, r, x);
        }
        if(k == 3) {
            int l, r; ans = 0;
            scanf("%lld %lld", &l, &r);
            query(1, l, r);
            printf("%lld\n", ans % MOD);
        }
    } return 0;
}

原文地址：https://www.cnblogs.com/Nicoppa/p/11445895.html