[BZOJ 3821]玄学

一、题目

UOJ的传送门

二、解法

离线的做法

这道题强制在线，但是没有关系，我们讲一下离线的做法，会很有启发。
考靠把询问和修改都放在$[1,n]$上面，我们从$1$扫到$n$，设当前点为$x$，我们维护包含当前点的操作为一颗线段树（树上节点表示连续的操作编号），当$x$为某一个操作的左端点时，加入这个操作的$a,b$，如果是某一个操作的右端点，删去这个操作的$a,b$。然后我们回答所有位置为$x$的问题，直接在线段树上询问即可。至于操作本质上是一个一次函数，合并的话很简单（$ax+b,cx+d$合并成为$acx+(bc+d)$），时间复杂度$O(n\log n)$，由于这不是正解，代码自然没有，讲解仅供参考。

在线：离散+归并的运用

在线就没办法了，只能够暴力存下所有信息。我们维护一个线段树，每一个节点表示$[l,r]$的操作之后，原数列所有$a,b$不同段的状态（注意，这里我们很暴力地把一个数组状态存了下来，因为一个点的状态数是$O(L)$的，所以总的状态是$O(n\log n)$的）。

定义如下，考虑修改。我们先改底层的状态，然后其他的状态能修改的条件是$r=id$（其中$id$是当前修改的操作编号），因为我们要用到一个点的状态必需要它操作区间里面都凑齐了，所以每个点都只会更新一次，那么复杂度是$O(n\log n)$。至于更新的话，由于原来用来合并的两个段都可以看做有序的，所以我们可以用类似于归并排序的方法更新。

考虑询问，我们找到线段树上$\log n$个点，找到询问的$x$对应的段的$a,b$，把这些$a,b$合并之后得到总的$a,b$，即可算出答案，这里需要套二分找位置，时间复杂度$O(n{\log }^{2}n)$。

#include <cstdio>
#include <iostream>
using namespace std;
const int M = 100005;
int read()
{
	int x=0,flag=1;char c;
	while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
	while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
	return x*flag;
}
int n,m,q,type,cnt,ans,ql,qr,qa,qb,ra,rb,a[M];
int L[4*M],R[4*M],dl[50*M],dr[50*M],da[50*M],db[50*M];
void up(int x)
{
	int i=L[x<<1],j=L[x<<1|1],cur=1;
	L[x]=cnt+1;
	while(i<=R[x<<1] || j<=R[x<<1|1])
	{
		dl[++cnt]=cur;
		dr[cnt]=min(dr[i],dr[j]);
		da[cnt]=1ll*da[i]*da[j]%m;
		db[cnt]=(1ll*da[j]*db[i]+db[j])%m;
		cur=dr[cnt]+1;
		if(dr[i]<dr[j]) i++;
		else if(dr[i]>dr[j]) j++;
		else i++,j++;
	}
	R[x]=cnt;
}
void modify(int i,int l,int r,int id)
{
	if(l==r)
	{
		L[i]=cnt+1;
		if(ql>1) dl[++cnt]=1,dr[cnt]=ql-1,da[cnt]=1;
		dl[++cnt]=ql,dr[cnt]=qr,da[cnt]=qa,db[cnt]=qb;
		if(qr<n) dl[++cnt]=qr+1,dr[cnt]=n,da[cnt]=1;
		R[i]=cnt; 
		return ;
	}
	int mid=(l+r)>>1;
	if(mid>=id) modify(i<<1,l,mid,id);
	else modify(i<<1|1,mid+1,r,id);
	if(r==id) up(i);
}
int Find(int l,int r)
{
	int ans=0;
	while(l<=r)
	{
		int mid=(l+r)>>1;
		if(dr[mid]<qa) l=mid+1;
		else r=mid-1,ans=mid;
	}
	return ans;
}
void ask(int i,int l,int r)
{
	if(ql>r || l>qr) return ;
	if(ql<=l && r<=qr)
	{
		int k=Find(L[i],R[i]);
		ra=1ll*ra*da[k]%m;rb=(1ll*da[k]*rb+db[k])%m;
		return ;
	}
	int mid=(l+r)>>1;
	ask(i<<1,l,mid);
	ask(i<<1|1,mid+1,r);
}
int main()
{
	type=read();n=read();m=read();
	for(int i=1;i<=n;i++)
		a[i]=read();
	q=read();int t=0;
	for(int i=1;i<=q;i++)
	{
		int op=read();ql=read();qr=read();qa=read();
		if(type&1) ql^=ans,qr^=ans;
		if(op==1)
		{
			qb=read();
			modify(1,1,1e5,++t);
		}
		else
		{
			if(type&1) qa^=ans;
			ra=1;rb=0;
			ask(1,1,1e5);
			ans=(1ll*ra*a[qa]+rb)%m;
			printf("%d\n",ans);
		}
	}
}