思路：扩展欧拉定理
提交：(geq5)次
错因：快速幂时刚开始没有判断(a)是否大于(p)
题解：
用树状数组维护差分，查询时暴力从左端点的第一个数向右端点递归，若递归时发现指数变为(1)，则指数返回(1)；若递归出右端点，指数也返回(1)；
#pragma GCC optimize (3)
#include<cstdio>
#include<iostream>
#define ll long long
#define R register ll
using namespace std;
namespace Luitaryi {
static char B[1<<15],*S=B,*T=B;
#define getchar() (S==T&&(T=(S=B)+fread(B,1,1<<15,stdin),S==T)?EOF:*S++)
template<class I> inline I g(I& x) { x=0; register I f=1;
register char ch; while(!isdigit(ch=getchar())) f=ch=='-'?-1:f;
do x=x*10+(ch^48); while(isdigit(ch=getchar())); return x*=f;
} const int N=500010,M=2e7;
int n,m,cnt,p[M/2],phi[M+10],a[N];
ll c[N]; bool v[M+10];
inline void PRE() { phi[1]=1;
for(register int i=2;i<=M;++i) {
if(!v[i]) p[++cnt]=i,phi[i]=i-1;
for(register int j=1;j<=cnt&&i*p[j]<=M;++j) {
v[i*p[j]]=true;
if(i%p[j]==0) {
phi[i*p[j]]=phi[i]*p[j]; break;
} phi[i*p[j]]=phi[i]*(p[j]-1);
}
}
}
inline void add(int x,int d) {for(;x<=n+4;x+=x&-x) c[x]+=d;}
inline ll query(int x) { R ret=0; for(;x;x-=x&-x) ret+=c[x]; return ret;}
inline ll qpow(ll a,ll p,int M) { R ret=1;
register bool flg=false,flg1=false;
if(a>=M) flg1=true,a%=M;//先判a
while(p) { if(p&1) {
ret*=a; flg|=flg1;
if(ret>=M) flg=true,ret%=M;
} a*=a; if(a>=M) flg1=true,a%=M; p>>=1;
} return ret+(flg?M:0);//指数是否大于模数
}
inline int solve(int l,int r,int x) {
if(x==1) return 1; //上面是1指数返回1
if(l>r) return 1; //到区间的最后，还是指数返回1
R cur=query(l)+a[l],p=solve(l+1,r,phi[x]);
return qpow(cur,p,x);
}
inline void main() { freopen("in.in","r",stdin); freopen("out.out","w",stdout);
PRE(); g(n),g(m); for(register int i=1;i<=n;++i) g(a[i]);
for(register int i=1,x,l,r,d;i<=m;++i) {
g(x),g(l),g(r),g(d); if(x&1) add(l,d),add(r+1,-d);
if(!(x&1)) printf("%d
",solve(l,r,d)%d);
}
}
} signed main() {Luitaryi::main(); return 0;}
2019.08.23
77