线段树

题目大意：一个线段树支持下面4中操作。

1.对于一个区间[l,r]把区间中的每个元素加c；

2.对于一个区间[l,r]把区间中的每个元素减c;

3.对于一个区间[l,r]把区间中的每个元素置c；

4.对于一个区间[l,r]求区间的和，并把这个和加到列表的每个元素上。

#include <bits/stdc++.h>
using namespace std;

const int imax_n = 10005;
const long long PP = 1000000009;

struct Node
{
    int l, r;
    long long sum;
    long long val;
    long long c;
    long long s;
}f[8 * imax_n];

long long a[imax_n];
int n, m;

void build(int t, int l, int r, long long *a)
{
    // printf("build %d %d %d
", t, l, r);
    f[t].l = l;
    f[t].r = r;
    f[t].val = 0LL;
    f[t].c = 0LL;
    f[t].s = 0LL;
    if (l == r)
    {
        // printf("add = %d
", a[l-1]);
        f[t].val = a[l-1];
        f[t].sum = a[l-1];
        // printf("%lld
", f[t].sum);
        return ;
    }
    int mid = (l + r) / 2;
    build(2*t, l , mid, a);
    build(2*t+1, mid+1, r, a);
    f[t].sum = (f[2*t].sum + f[2*t+1].sum) % PP;
}

void pushDown(int t)
{
    if (f[t].s)
    {
        f[2*t].s = f[t].s;
        f[2*t+1].s = f[t].s;
        f[2*t].val = f[t].val;
        f[2*t].sum = f[t].val * (f[2*t].r - f[2*t].l + 1) % PP;
        f[2*t + 1].val = f[t].val;
        f[2*t + 1].sum = f[t].val * (f[2*t + 1].r - f[2*t+1].l + 1) % PP;
        f[t].s = 0LL;
        f[t].c = 0LL;
        f[t].val = 0LL;
    }

    if (f[t].c)
    {
        f[2*t].c = (f[2*t].c + f[t].c) % PP;
        f[2*t+1].c  = (f[2*t+1].c + f[t].c) %PP;
        f[2*t].sum = (f[2*t].sum + f[t].c * (f[2*t].r - f[2*t].l + 1)) % PP;
        f[2*t + 1].sum = (f[2*t+1].sum + f[t].c * (f[2*t+1].r - f[2*t+1].l + 1)) % PP;
        f[t].c = 0;
    }
}

void pushUp(int t)
{
    f[t].sum = (f[2*t].sum + f[2*t+1].sum) % PP;
}

void modify(int t, int l, int r, int c)
{
    if (f[t].l == l && f[t].r == r)
    {
        if (f[t].s)
            f[t].val = (f[t].val + c) % PP;
        else
            f[t].c = (f[t].c + c) % PP;
        f[t].sum = (f[t].sum + (long long )c * (f[t].r - f[t].l + 1)) % PP;
        return ;
    }

    pushDown(t);
    int mid = (f[t].l + f[t].r) / 2;
    if (r<=mid) modify(2*t, l, r, c);
    else if (l > mid) modify(2*t+1, l, r, c);
    else
    {
        modify(2*t, l, mid, c);
        modify(2*t+1, mid + 1, r, c);
    }
    pushUp(t);
}

void set_(int t, int l, int r, int val)
{
    // printf("set_ %d %d %d
", t, l, r);
    if (f[t].l == l && f[t].r == r)
    {
        f[t].s = 1LL;
        f[t].val = (long long )val;
        f[t].sum = f[t].val* (f[t].r - f[t].l + 1) % PP;
        return;
    }
    pushDown(t);
    int mid = (f[t].l + f[t].r) / 2;
    if (r <= mid) set_(2*t, l, r, val);
    else if (l > mid) set_(2*t+1, l, r, val);
    else
    {
        set_(2*t, l, mid, val);
        set_(2*t+1, mid+1, r, val);
    }
    pushUp(t);
}

long long query(int t, int l, int r)
{
    if (f[t].l == l && f[t].r == r)
    {
        return f[t].sum % PP;
    }
    pushDown(t);
    int mid = (f[t].l + f[t].r) / 2;
    if (r<=mid) return query(2*t, l, r);
    else if (l > mid) return query(2*t + 1, l, r);
    else
    {
        return (query(2*t, l, mid) + query(2*t+1, mid+1, r)) % PP;
    }
}

void print_tree(int t)
{
    printf("node: %d l=%d r=%d s=%lld c=%lld sum=%lld
", t, f[t].l, f[t].r, f[t].s, f[t].c, f[t].sum);
    if (f[t].l == f[t].r)
    {
        return;
    }
    print_tree(2*t);
    print_tree(2*t+1);
}

void dfs(int t, vector<int> &ans)
{

    if (f[t].l == f[t].r)
    {
        ans.push_back((int)(f[t].sum%PP));
        return ;
    }
    pushDown(t);
    dfs(2*t, ans);
    dfs(2*t+1, ans);
    pushUp(t);
}

int main()
{
    scanf("%d%d", &n, &m);
    for (int i = 0; i < n; ++i)
    {
        scanf("%lld", &a[i]);
    }
    build(1, 1, n, a);
    // print_tree(1);
    long long op4 = 0LL;
    for (int i = 0; i < m; ++i)
    {
        int op, x, y, c;
        scanf("%d%d%d%d", &op, &x, &y, &c);
        if (op==1)
        {
            modify(1, x, y, c);
            // print_tree(1);
        }
        else if (op==2)
        {
            modify(1, x, y, -c);
            // print_tree(1);
        }
        else if (op==3)
        {
            set_(1, x, y, c);
            // print_tree(1);
        }
        else
        {
            op4 = query(1, x, y) % PP;
            // printf("qeury %d %d op4 = %lld
", x, y, op4);
            modify(1, 1, n, op4);
        }
    }
    vector<int> ans;
    dfs(1, ans);
    for (int i = 0; i < ans.size(); ++i)
    {
        printf("%d ", ans[i]);
    }
    printf("
");
    return 0;
}