洛谷P3373 【模板】线段树 2

题目描述

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

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的结果。

 

输入输出样例

输入样例#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

说明

时空限制：1000ms,128M

数据规模：

对于30%的数据：N<=8，M<=10

对于70%的数据：N<=1000，M<=10000

对于100%的数据：N<=100000，M<=100000

（数据已经过加强^_^）

样例说明：

故输出应为17、2（40 mod 38=2）

---

这题做了两天，，仍然找不到哪错了，好伤心，先挂着，以后回头看，走过路过的coder，顺便帮我看一下吧，生活还要继续，擦干眼泪继续干

#include<bits/stdc++.h>
#define lom long long 
#define maxn 100000
using namespace std;
struct point
{ 
	int l,r; 
	lom sum,mul,pul;
};
lom mod;
int a[maxn+5];
point tree[maxn*4+5];
void bulid(int x,int l,int r)
{
	tree[x].l=l; tree[x].r=r;tree[x].mul=1;
	if(l==r)
	{
		tree[x].sum=a[l];
		tree[x].mul=1;
		return ;
	}
	
	int mid=l+r>>1;
	bulid(2*x,l,mid);
	bulid(2*x+1,mid+1,r);
	
	tree[x].sum=tree[x<<1].sum+tree[x<<1|1].sum%mod;
}
void spread(int x)
{
	if(tree[x].mul!=1)
	{
		tree[2*x].mul=(tree[x].mul*tree[2*x].mul)%mod;
		tree[2*x+1].mul=(tree[x].mul*tree[2*x+1].mul)%mod;
		tree[2*x].pul=(tree[2*x].pul+tree[x].pul)%mod;
		tree[2*x+1].pul=(tree[2*x+1].pul+tree[x].pul)%mod;
		tree[2*x].sum=(tree[2*x].sum*tree[x].mul)%mod;
		tree[2*x+1].sum=(tree[2*x+1].sum*tree[x].mul)%mod;
		tree[x].mul=1;
	}
	if(tree[x].pul!=0)
	{
		tree[2*x].pul=(tree[2*x].pul+tree[x].pul)%mod;
		tree[2*x+1].pul=(tree[2*x+1].pul+tree[x].pul)%mod;
		tree[2*x+1].sum=(tree[2*x+1].sum+tree[x].pul*(tree[2*x+1].r-tree[2*x+1].l+1))%mod;
		tree[2*x].sum=(tree[2*x].sum+tree[x].pul*(tree[2*x].r-tree[2*x].l+1))%mod;
		tree[x].pul=0;
	}
}
void update_mul(int x,int l,int r,int m)
{
	if(l<=tree[x].l&&tree[x].r<=r)
	{
		tree[x].pul=(m*tree[x].pul)%mod;
		tree[x].mul=(m*tree[x].mul)%mod;
		tree[x].sum=(m*tree[x].sum)%mod;
		return ;
	}
	
	spread(x);
	int mid=tree[x].l+tree[x].r>>1;
	if(l<=mid) update_mul(2*x,l,r,m);
	if(r>mid) update_mul(2*x+1,l,r,m);
	tree[x].sum=tree[x<<1].sum+tree[x<<1|1].sum%mod;
}
void update_pul(int x,int l,int r,int m)
{
	if(l<=tree[x].l&&tree[x].r<=r)
	{
		tree[x].sum=(tree[x].sum+(lom)(tree[x].r-tree[x].l+1)*m)%mod;
		tree[x].pul=(tree[x].pul+m)%mod;
		return;
	}
	spread(x);
	int mid=tree[x].l+tree[x].r>>1;
	if(l<=mid) update_pul(2*x,l,r,m);
	if(r>mid) update_pul(2*x+1,l,r,m);
	tree[x].sum=tree[x<<1].sum+tree[x<<1|1].sum%mod;
}

lom find(int x,int l,int r)
{
	if(l<=tree[x].l&&tree[x].r<=r)
	return tree[x].sum%mod;
	spread(x);
	lom ans=0;
	int mid=tree[x].l+tree[x].r>>1;
	if(l<=mid) ans+=find(2*x,l,r)%mod;
	if(r>mid)  ans+=find(2*x+1,l,r)%mod;
	return ans%mod;
}
void print(int n)
{
	for(int i=1;i<=n*4;i++)
	cout<<tree[i].sum<<' ';
	cout<<endl<<endl;
 } 
int main()
{
	int n,m;
	cin>>n>>m>>mod;
	for(int i=1;i<=n;i++)
	cin>>a[i];
	bulid(1,1,n);
	print(n);
	while(m--)
	{
		int order,x,y,k;
		cin>>order>>x>>y;
		if(order==1) 
		{
			cin>>k;
			update_mul(1,x,y,k); 
			print(n);
		}
		else if(order==2)
		{
			cin>>k;
			update_pul(1,x,y,k);
			print(n);
		}
		else cout<<find(1,x,y)<<endl;
	}
	return 0;
}