分块【乘法和加法更新模板】

这次我的分块终于跑的比线段树慢了，好开心！

那我为啥要写怎么慢的算法啊！

因为分块可以维护一些乱七八糟更新的东西。

对于乘法和加法来说，分块也和线段树一样，需要mul和add两个标记。

每一次update都需要pushdown(其实不应该叫pushdown，应该是reset，目的是把标记更新上去)。

对于散块来说，每一次都需要pushdown, 只需要更新a和sum,不用打标记。

mul的时候要更新add，因为乘法的优先级比较高。

下面是代码：mul一定要设置成1，belong我搞错了好几次，应该是从1到n更新belong。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 1e5+7;
ll block,num,l[maxn],r[maxn],belong[maxn],mul[maxn],add[maxn];
ll a[maxn],sum[maxn];
int n,m;
ll p;
void build()
{
	block = sqrt(n);
	num = n/block;
	if(n%block) num++;
	for(int i=1;i<=num;i++)
	{
		l[i] = (i-1)*block+1;
		r[i] = i*block;
		mul[i] = 1;
	}
	r[num] = n;
	for(int i=1;i<=n;i++)
	{
		belong[i] = (i-1)/block+1;
	}
	for(int i=1;i<=num;i++)
	{
		for(int j=l[i];j<=r[i];j++)
		{
			sum[i] = (sum[i]+a[j])%p;
		}	
	}
}
void pushdown(int x)
{
	for(int i=l[x];i<=r[x];i++)
	{
		a[i] = (a[i]*mul[x]+add[x])%p;
	}
	mul[x] = 1;
	add[x] = 0;
}
void Mul(int x,int y,int v)
{
	if(belong[x]==belong[y])
	{
		pushdown(belong[x]);
		for(int i=x;i<=y;i++)
		{
			sum[belong[x]] += (v-1)*a[i]%p;   //orginally, sum[belong[x]] == sigma(ai), then now, sum[belong[x]] += (v-1)*a[i];
			a[i] = (a[i]*v)%p;
		}
		return;
	}
	pushdown(belong[x]);
	for(int i=x;i<=r[belong[x]];i++)
	{
		sum[belong[x]] += (v-1)*a[i]%p;
		a[i] = (a[i]*v)%p;
	}
	for(int i=belong[x]+1;i<belong[y];i++)
	{
		mul[i] = (mul[i]*v)%p;
		add[i] = (add[i]*v)%p;
		sum[i] = (sum[i]*v)%p;
	}
	pushdown(belong[y]);
	for(int i=l[belong[y]];i<=y;i++)
	{
		sum[belong[y]] += (v-1)*a[i]%p;
		a[i] = (a[i]*v)%p;
	}
}
void Add(int x,int y,int v)
{
	if(belong[x]==belong[y])
	{
		pushdown(belong[x]);
		for(int i=x;i<=y;i++)
		{
			a[i] = (a[i]+v)%p;
		}
		sum[belong[x]] = (sum[belong[x]]+(y-x+1)*v)%p;
		return;
	}
	pushdown(belong[x]);
	for(int i=x;i<=r[belong[x]];i++)
	{
		a[i] = (a[i]+v)%p;
	}
	sum[belong[x]] = (sum[belong[x]]+(r[belong[x]]-x+1)*v)%p;
	for(int i=belong[x]+1;i<belong[y];i++)
	{
		sum[i] = (sum[i]+block*v)%p;
		add[i] = (add[i]+v)%p;
	}
	pushdown(belong[y]);
	for(int i=l[belong[y]];i<=y;i++)
	{
		a[i] = (a[i]+v)%p;
	}
	sum[belong[y]] = (sum[belong[y]]+(y-l[belong[y]]+1)*v)%p;
}
ll query(int x,int y)
{
	if(belong[x]==belong[y])
	{
		ll ans = 0;
		for(int i=x;i<=y;i++)
		{
			ans = (ans+a[i]*mul[belong[x]]+add[belong[x]])%p;
		}
		return ans;
	}
	ll ans = 0;
	for(int i=x;i<=r[belong[x]];i++)
	{
		ans = (ans+a[i]*mul[belong[x]]+add[belong[x]])%p;
	}
	for(int i=belong[x]+1;i<belong[y];i++)
	{
		ans = (ans+sum[i])%p;
	}
	for(int i=l[belong[y]];i<=y;i++)
	{
		ans = (ans+a[i]*mul[belong[y]]+add[belong[y]])%p;
	}
	return ans;
}
int main()
{
	scanf("%d%d%lld",&n,&m,&p);
	for(int i=1;i<=n;i++)
	{
		scanf("%lld",a+i);
	}
	build();
	int op,x,y,z;
	for(int i=0;i<m;i++)
	{
		scanf("%d",&op);
		if(op==1)
		{
			scanf("%d%d%d",&x,&y,&z);
			Mul(x,y,z);
		}
		else if(op==2)
		{
			scanf("%d%d%d",&x,&y,&z);
			Add(x,y,z);
		}
		else
		{
			scanf("%d%d",&x,&y);
			printf("%lld\n",query(x,y));
		}
	}
	return 0;
}