题目描述

如题，已知一个数列，你需要进行下面三种操作：

  • 将某区间每一个数乘上 xx

  • 将某区间每一个数加上 xx

  • 求出某区间每一个数的和

输入格式

第一行包含三个整数 n,m,pn,m,p，分别表示该数列数字的个数、操作的总个数和模数。

第二行包含 nn 个用空格分隔的整数，其中第 ii 个数字表示数列第 ii 项的初始值。

接下来 mm 行每行包含若干个整数，表示一个操作，具体如下：

操作 11： 格式：1 x y k 含义：将区间 [x,y][x,y] 内每个数乘上 kk

操作 22： 格式：2 x y k 含义：将区间 [x,y][x,y] 内每个数加上 kk

操作 33： 格式：3 x y 含义：输出区间 [x,y][x,y] 内每个数的和对 pp 取模所得的结果

输出格式

输出包含若干行整数，即为所有操作 33 的结果。

输入输出样例

输入 #1

5 5 38
1 5 4 2 3
2 1 4 1
3 2 5
1 2 4 2
2 3 5 5
3 1 4

输出 #1复制

17
2

说明/提示

【数据范围】

对于 30\%30% 的数据：n \le 8n≤8，m \le 10m≤10
对于 70\%70% 的数据：n \le 10^3n≤103，m \le 10^4m≤104
对于 100\%100% 的数据：n \le 10^5n≤105，m \le 10^5m≤105

除样例外，p = 571373p=571373

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#pragma warning(disable:4996)
struct Node
{
	ll l, r;
	ll lazy = 0, sum = 0, mlz = 1;
}tree[500005];
ll input[500005];
int n, m, p;
void pushup(ll k)
{
	ll k1 = tree[k].mlz, k2 = tree[k].lazy;
	tree[k << 1].sum = (tree[k << 1].sum*k1 + k2 * (tree[k << 1].r - tree[k << 1].l + 1)) % p;
	tree[k << 1 | 1].sum = (tree[k << 1 | 1].sum*k1 + k2 * (tree[k << 1 | 1].r - tree[k << 1 | 1].l + 1)) % p;
	tree[k << 1].mlz = (tree[k << 1].mlz*k1) % p;
	tree[k << 1 | 1].mlz = (tree[k << 1 | 1].mlz*k1) % p;
	tree[k << 1].lazy = (tree[k << 1].lazy*k1 + k2) % p;
	tree[k << 1 | 1].lazy = (tree[k << 1 | 1].lazy*k1 + k2) % p;
	tree[k].lazy = 0;
	tree[k].mlz = 1;
	return;
}
void build(int k, int l, int r)
{
	tree[k].l = l, tree[k].r = r;
	if (l == r)
	{
		tree[k].sum = input[l]%p; return;
	}
	int mid = (l + r) >> 1;
	build(2 * k, l, mid);
	build(2 * k + 1, mid + 1, r);
	tree[k].sum = (tree[2 * k].sum + tree[2 * k + 1].sum) % p;
}
void add1(int k, int dis, int x)//单点修改
{
	if (tree[k].l == tree[k].r)
	{
		//cout << "!!!!!   " << tree[k].l<<" "<< tree[k].sum<<"  << endl;
		tree[k].sum += x;
		return;
	}
	int mid = (tree[k].l + tree[k].r) >> 1;
	if (mid >= dis)add1(2 * k, dis, x);
	else add1(2 * k + 1, dis, x);
	tree[k].sum = tree[2 * k].sum + tree[2 * k + 1].sum;
}
void add2(int k, ll l, ll r, ll x)
{
	if (tree[k].l == l && tree[k].r == r)
	{
		tree[k].sum += ((r - l + 1)*x) % p;
		tree[k].lazy = (tree[k].lazy + x) % p;
		return;
	}
	pushup(k);
	int mid = (tree[k].l + tree[k].r) >> 1;
	if (mid >= r)add2(2 * k, l, r, x);
	else if (mid < l)add2(2 * k + 1, l, r, x);
	else add2(2 * k, l, mid, x), add2(2 * k + 1, mid + 1, r, x);
	tree[k].sum = (tree[2 * k].sum + tree[2 * k + 1].sum) % p;
}
void mul(int k, ll l, ll r, ll x)
{
	if (tree[k].l == l && tree[k].r == r)
	{
		tree[k].sum = (tree[k].sum*x) % p;
		tree[k].mlz = (tree[k].mlz*x) % p;
		tree[k].lazy = (tree[k].lazy*x) % p;
		return;
	}
	pushup(k);
	int mid = (tree[k].l + tree[k].r) >> 1;
	if (mid >= r)mul(2 * k, l, r, x);
	else if (mid < l)mul(2 * k + 1, l, r, x);
	else mul(2 * k, l, mid, x), mul(2 * k + 1, mid + 1, r, x);
	tree[k].sum = (tree[2 * k].sum + tree[2 * k + 1].sum) % p;
}
ll query(int k, ll l, ll r)
{
	if (tree[k].l == l && tree[k].r == r)
	{
		return tree[k].sum%p;
	}
	pushup(k);
	int mid = (tree[k].l + tree[k].r) >> 1;
	if (mid >= r)return query(2 * k, l, r) % p;
	if (mid < l)return query(2 * k + 1, l, r) % p;
	return (query(2 * k, l, mid) % p + query(2 * k + 1, mid + 1, r) % p) % p;
}
int main()
{

	cin >> n >> m >> p;
	for (int i = 1; i <= n; i++)scanf("%d", &input[i]);
	build(1, 1, n);
	//for (int i = 1; i <= 7; i++)cout << tree[i].sum << " ";
	while (m--)
	{
		int a;
		scanf("%d", &a);
		if (a == 1)
		{
			int x, y, k;
			
			scanf("%d%d%d", &x, &y, &k);k %= p;
			mul(1, x, y, k);
		}
		else if (a == 2)
		{
			int x, y, k;
			
			scanf("%d%d%d", &x, &y, &k);k %= p;
			add2(1, x, y, k);
		}
		else
		{
			int x, y;
			scanf("%d%d", &x, &y);
			cout << query(1, x, y)%p << endl;
		}
	}



}