题意:
luogu
\text{luogu}
luogu的智能推荐是个好东西。
题目传送门:
https://www.luogu.org/problem/P3934
思路:
弱化版题目详见:
CF906D Power Tower
\text{CF906D Power Tower}
CF906D Power Tower。
我们发现这道题只需要满足区间修改单点查询即可。
于是用差分树状数组即可。
这个树状数组是不能取模的,而题目保证用
long long
\text{long long}
long long可以存下。
这题数据较强,因此我在这道题中发现了上一道题代码的
bug
\text{bug}
bug——快速幂开始的时候就要将
x
x
x搞一下,防止
x
x
x过大,相乘爆
long long
\text{long long}
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);
}
}