维护区间加法和区间乘法

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

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

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

输入输出样例
输入 #1
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
输出 #1
2
35
8

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long ll;
const int maxn = 100005;
int a,b;
ll n,m,mod,x,ans;
struct node
{
    int l,r;
    ll fm,fs,sum;//fm是乘法的懒标记 fs是加法的懒标记
}tr[maxn*4+1];

void down(int k)
{
    tr[k*2].sum = (tr[k*2].sum*tr[k].fm%mod + (tr[k*2].r - tr[k*2].l + 1)*tr[k].fs%mod)%mod;
    tr[k*2+1].sum = (tr[k*2+1].sum*tr[k].fm%mod + (tr[k*2+1].r - tr[k*2+1].l + 1)*tr[k].fs%mod)%mod;
    tr[k*2].fm *= tr[k].fm; tr[k*2].fm %= mod;
    tr[k*2+1].fm *= tr[k].fm; tr[k*2+1].fm %= mod;
    tr[k*2].fs = tr[k*2].fs*tr[k].fm%mod + tr[k].fs; tr[k*2].fs %= mod;
    tr[k*2+1].fs = tr[k*2+1].fs*tr[k].fm%mod + tr[k].fs; tr[k*2+1].fs %= mod;;
    tr[k].fm = 1;
    tr[k].fs = 0;
}

void build(int k,int ll,int rr)
{
    tr[k].l = ll,tr[k].r = rr;
    tr[k].fm = 1,tr[k].fs = 0;
    if(tr[k].l == tr[k].r)
    {
        scanf("%lld",&tr[k].sum);
        tr[k].sum %= mod;
        return;
    }
    int mid = (ll+rr)/2;
    build(k*2,ll,mid);
    build(k*2+1,mid+1,rr);
    tr[k].sum = (tr[k*2].sum + tr[k*2+1].sum)%mod;
}

void upm(int k)
{
    if(tr[k].l >= a && tr[k].r <= b)
    {
        tr[k].sum *= x; tr[k].sum %= mod;
        tr[k].fm *= x; tr[k].fm %= mod;
        tr[k].fs *= x; tr[k].fs %= mod;
        return;
    }
    down(k);
    int m = (tr[k].l+tr[k].r)/2;
    if(a <= m)
        upm(k*2);
    if(b > m)
        upm(k*2+1);
    tr[k].sum = (tr[k*2].sum + tr[k*2+1].sum)%mod;
}

void ups(int k)
{
    if(tr[k].l >= a && tr[k].r <= b)
    {
        tr[k].sum += x*(tr[k].r - tr[k].l + 1);
        tr[k].sum %= mod;
        tr[k].fs += x; tr[k].fs %= mod;
        return;
    }
    down(k);
    int m = (tr[k].l+tr[k].r)/2;
    if(a <= m)
        ups(k*2);
    if(b > m)
        ups(k*2+1);
    tr[k].sum = (tr[k*2].sum + tr[k*2+1].sum)%mod;
}

void query(int k)
{
    if(tr[k].l >= a && tr[k].r <= b)
    {
        ans += tr[k].sum;
        ans %= mod;
        return;
    }
    down(k);
    int m = (tr[k].l+tr[k].r)/2;
    if(a <= m)
        query(k*2);
    if(b > m)
        query(k*2+1);
}

int main()
{
    scanf("%lld%lld%lld",&n,&m,&mod);
    build(1,1,n);
    while(m--)
    {
        int p;
        scanf("%d",&p);
        ans = 0;
        if(p == 1)
        {
            scanf("%d%d%lld",&a,&b,&x);
            upm(1);
        }
        if(p == 2)
        {
            scanf("%d%d%lld",&a,&b,&x);
            ups(1);
        }
        if(p == 3)
        {
            scanf("%d%d",&a,&b);
            query(1);
            printf("%lld\n",ans%mod);
        }
    }
}