uoj#46. 【清华集训2014】玄学

难得的1A。。。

今天上午红太阳Lych大讲了归并，分块，推和可并推，马上练了这题。

题意很简单对于一个数列有操作使x~y区间内的数*a+b，多次询问有x~y操作时z的值（对m取模），很显然*a+b的操作满足结合律且多次*a+b可合并为一个*c+d，不满足区间减法（m不为质数）。

对于询问一段操作，可以考虑对时间建线段树，l~r的时间内对数列的修改是O（r-l）段，于是归并，然后就好了。。。

#include<iostream>
#include<cstdio>
#define N 400005
#define ll long long
using namespace std;
int Hz,n,m,x,y,z,a,b,q,s;
int v[N],cnt;
struct type{int l,r;}c[N];
struct node{int x,y,a,b;}p[N*10];
void gb(int k)
{
	int l1=c[k<<1].l;
	int l2=c[k<<1|1].l;
	int l=1;
	for (c[k].l=cnt+1;l<=n;)
	{
		if (p[l1].y<p[l2].y)
		{
			p[++cnt]=(node){l,p[l1].y,(ll)p[l1].a*p[l2].a%m,((ll)p[l1].b*p[l2].a+p[l2].b)%m};
			l=p[l1].y+1;l1++;continue;
		}
		if (p[l1].y>p[l2].y)
		{
			p[++cnt]=(node){l,p[l2].y,(ll)p[l1].a*p[l2].a%m,((ll)p[l1].b*p[l2].a+p[l2].b)%m};
			l=p[l2].y+1;l2++;continue;
		}
		if (p[l1].y==p[l2].y)
		{
			p[++cnt]=(node){l,p[l1].y,(ll)p[l1].a*p[l2].a%m,((ll)p[l1].b*p[l2].a+p[l2].b)%m};
			l=p[l1].y+1;l1++;l2++;continue;
		}
	}
	c[k].r=cnt;
}
void update(int l,int r,int z)
{
	int mid;
	while (l<r)
	{
		int mid=(l+r)>>1;
		if (p[mid].y<z) l=mid+1;else r=mid;
	}
	a=(ll)a*p[l].a%m;
	b=((ll)b*p[l].a+p[l].b)%m;
}
void mdy(int k,int l,int r,int t)
{
	if (l==r)
	{
		c[k].l=++cnt;
		if (x!=1) p[cnt]=(node){1,x-1,1,0};
		p[++cnt]=(node){x,y,a,b};
		if (y!=n) p[++cnt]=(node){y+1,n,1,0};
		c[k].r=cnt;
		return;
	}
	int mid=(l+r)>>1;
	if (t<=mid) mdy(k<<1,l,mid,t);
	else mdy(k<<1|1,mid+1,r,t);
	if (t==r) gb(k);
}
void qry(int k,int l,int r,int x,int y)
{
	if (x<=l&&y>=r)
	{
		update(c[k].l,c[k].r,z);
		return;
	}
	int mid=(l+r)>>1;
	if (x<=mid) qry(k<<1,l,mid,x,y);
	if (y>mid) qry(k<<1|1,mid+1,r,x,y);
}
int main()
{
	scanf("%d",&Hz);Hz=Hz%2;
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;i++)
		scanf("%d",&v[i]);
	scanf("%d",&q);int sum=0,LAns=0,Mx=min(100000,q);
	for (int i=1;i<=q;i++)
	{
		scanf("%d",&s);
		if (s==1)
		{
			scanf("%d%d%d%d",&x,&y,&a,&b);
			x^=LAns*Hz;y^=LAns*Hz;
			mdy(1,1,Mx,++sum);
		}
		else
		{
			scanf("%d%d%d",&x,&y,&z);
			x^=LAns*Hz;y^=LAns*Hz;z^=LAns*Hz;
			a=1;b=0;qry(1,1,Mx,x,y);
			LAns=((ll)a*v[z]+b)%m;
			printf("%d\n",LAns);
		}
	}
	return 0;
}