##Segment Tree

馅断术模板啦~

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

1. 将某区间每一个数乘上x

2. 将某区间每一个数加上x

3. 求出某区间每一个数的和

输入格式

第一行包含三个整数N、M、P，分别表示该数列数字的个数、操作的总个数和模数。

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

接下来M行每行包含3或4个整数，表示一个操作，具体如下：

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

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

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

输出格式

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

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <memory.h>

#define ll long long 
using namespace std;

const int maxn = 100005;
const int maxm = 1000005;
ll a[maxn];
ll sum[maxm];
ll addv[maxm];
ll mulv[maxm];
int p;
inline ll read()
{
	char c = getchar(); ll x = 0;
	for (; !isdigit(c); c = getchar());
	for (; isdigit(c); c = getchar())
		x = x * 10 + (c - '0');
	return x;
}


void  cons(int o, int L, int R)
{
	addv[o] = 0;
	mulv[o] = 1;
	if (L == R)
		sum[o] = a[L];
	else
	{
		int mid = (L + R) / 2;
		cons(o * 2, L, mid);
		cons(o * 2 + 1, mid + 1, R);
		sum[o] = sum[o * 2] + sum[o * 2 + 1];
	}
	sum[o] %= p;
	return;
}


inline void maintain(int o, int L, int R)
{
	int mid = (L + R) / 2;
	sum[o * 2] = (sum[o * 2] * mulv[o] + addv[o] * (mid - L + 1)) % p;
	sum[o * 2 + 1] = (sum[o * 2 + 1] * mulv[o] + addv[o] * (R - mid)) % p;

	mulv[o * 2] = (mulv[o] * mulv[o * 2]) % p;
	mulv[o * 2 + 1] = (mulv[o] * mulv[o * 2 + 1]) % p;
	addv[o * 2] = (addv[o * 2] * mulv[o] + addv[o]) % p;
	addv[o * 2 + 1] = (addv[o * 2 + 1] * mulv[o] + addv[o]) % p;

	addv[o] = 0;
	mulv[o] = 1;
}


ll ql, qr;
void u_add(int o, int L, int R, int d)
{
	if (R<ql || L>qr)return;
	if (ql <= L && R <= qr)
	{
		addv[o] = (addv[o] + d) % p;
		sum[o] = (sum[o] + d * (R - L + 1)) % p;
		return;
	}

	maintain(o, L, R);
	int mid = (L + R) / 2;
	u_add(2 * o, L, mid, d);
	u_add(2 * o + 1, mid + 1, R, d);
	sum[o] = (sum[o * 2] + sum[o * 2 + 1]) % p;
	return;



}

void u_mul(int o, int L, int R, int m)
{
	if (R<ql || L>qr)return;
	if (ql <= L && R <= qr)
	{
		sum[o] = (sum[o] * m) % p;
		addv[o] = (addv[o] * m) % p;
		mulv[o] = (mulv[o] * m) % p;
		return;
	}
	maintain(o, L, R);
	int mid = (L + R) / 2;
	u_mul(2 * o, L, mid, m);
	u_mul(2 * o + 1, mid + 1, R, m);
	sum[o] = (sum[o * 2] + sum[o * 2 + 1]) % p;
	return;



}


ll query(int o, int L, int R)//intial d=0,m=1;
{
	if (R<ql || L>qr)return 0;
	if (ql <= L && R <= qr)
		return sum[o];

	ll s1 = 0, s2 = 0;
	maintain(o, L, R);
	int mid = (L + R) / 2;
	s1 = query(2 * o, L, mid);
	s2 = query(2 * o + 1, mid + 1, R);
	return (s1 + s2) % p;
}

int main()
{
	memset(mulv, 1, sizeof(mulv));
	int n, m;
	cin >> n >> m >> p;
	for (int i = 1; i <= n; i++)
		a[i] = read();
	cons(1, 1, n);
	for (int i = 1; i <= m; i++)
	{
		ll g = 0;
		g = read();
		ql = read();
		qr = read();
		if (g == 1)
		{
			int k = read();
			u_mul(1, 1, n, k);
			continue;
		}
		if (g == 2)
		{
			int k = read();
			u_add(1, 1, n, k);
			continue;
		}
		if (g == 3)
		{
			ll s = query(1, 1, n);
			printf("%d\n", s);
			continue;
		}
	}


	return 0;
}