luogu P3934 Nephren Ruq Insania

題意：

$\text{luogu}$的智能推薦是個好東西。

題目傳送門：

https://www.luogu.org/problem/P3934

思路：

弱化版題目詳見：$\text{CF906D Power Tower}$。
我們發現這道題只需要滿足區間修改單點查詢即可。
於是用差分樹狀數組即可。
這個樹狀數組是不能取模的，而題目保證用$\text{long long}$可以存下。
這題數據較強，因此我在這道題中發現了上一道題代碼的$\text{bug}$——快速冪開始的時候就要將$x$搞一下，防止$x$過大，相乘爆$\text{long long}$。

代碼：

#pragma GCC optimize("Ofast")
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
#define lowbit(x) ((x)&(-(x)))
#define I inline
#define R register
using namespace std;
	LL tr[1000010];
	int phi[20000010],prime[20000010];
	int n,m;
I void init(int ma)
{
	int t=0;
	phi[1]=1;
	for(R int i=1;i<=ma;i++)
	{
		if(!phi[i]) phi[i]=i-1,prime[++t]=i;
		for(R int j=1;j<=t&&i*prime[j]<=ma;j++)
		{
			if(!(i%prime[j]))
			{
				phi[i*prime[j]]=phi[i]*prime[j];
				break;
			}
			phi[i*prime[j]]=phi[i]*phi[prime[j]];
		}
	}
}
I void add(int x,int y)
{
	for(;x<=n;x+=lowbit(x))
		tr[x]+=y;
}
I LL getsum(int x)
{
	LL sum=0;
	for(;x;x-=lowbit(x))
		sum+=tr[x];
	return sum;
}
I LL calc(LL x,int p)
{
	return x<p?x:x%p+p;
}
I LL ksm(LL x,int k,int p)
{
	LL tot=1;
	x=calc(x,p);
	for(;k;k>>=1)
	{
		if(k&1) tot=calc(tot*x,p);
		x=calc(x*x,p);
	}
	return tot;
}
LL solve(int x,int y,int p)
{
	if(x==y||p==1) return calc(getsum(x),p);
	return ksm(getsum(x),solve(x+1,y,phi[p]),p);
}
int main()
{
	int t,x,y=0,z;
	init(20000000);
	scanf("%d %d",&n,&m);
	for(R int i=1;i<=n;i++)
	{
		scanf("%d",&x);
		add(i,x-y);
		y=x;
	}
	for(R int i=1;i<=m;i++)
	{
		scanf("%d %d %d %d",&t,&x,&y,&z);
		if(t==1) add(x,z),add(y+1,-z); else printf("%lld\n",solve(x,y,z)%z);
	}
}