題目大意：給你一個序列，需要支持區間修改，以及查詢一段區間$a_{i}^{a_{i+1}^{a_{i+2}...}}mod\;p$的值，每次詢問的$p$的值不同

對于區間修改，由線段樹完成，沒什么好說的

對于查詢，利用"上帝與集合的正確用法"那道題的方法，不斷取$\phi(p)$降冪，那么最多迭代$log$層

由于$ai$不一定和$p$互質，需要使用拓展歐拉定理

$ans=ai^{Ans_{i+1}\;mod\;\phi(p)+Ans_{i+1}>=\phi(p)?\phi(p):0}$

每層都記錄$Ans_{i}$是否對這一層的$p$取過模，用于判斷Ans_{i+1}>=\phi(p)

?

特判比較多，注意特判的優先級

1.如果$ai mod p==0$，那么不論接下來是幾次冪都返回0

2.如果$ai==1$，那么不論$ai$的多少次冪都是1，同時清空已經取過模的標記

3.如果$i+1$層取過模，那么這一層的運算即使沒取模也要標記為取過模，因為實際值一定取模了

4.如果$i+1$層取過模，那么快速冪里需要加上$\phi(p)$，反之$i+1$沒取模一定不要加，不然答案是錯的

5.接下來在快速冪里更新第i層計算的取模標記

1 #include <cstdio> 2 #include <cstring> 3 #include <algorithm> 4 #define NN 501000 5 #define MM 20001000 6 #define maxn 20000000 7 #define ll long long 8 #define uint unsigned int 9 #define ull unsigned long long 10 using namespace std; 11 12 int n,m,g,r,root; 13 int gint() 14 { 15 int ret=0,fh=1;char c=getchar(); 16 while(c<'0'||c>'9'){if(c=='-')fh=-1;c=getchar();} 17 while(c>='0'&&c<='9'){ret=ret*10+c-'0';c=getchar();} 18 return ret*fh; 19 } 20 int use[MM],pr[2000000],phi[MM],cnt; 21 void Pre() 22 { 23 phi[1]=1; 24 for(int i=2;i<=maxn;i++) 25 { 26 if(!use[i]) pr[++cnt]=i,phi[i]=i-1; 27 for(int j=1;j<=cnt&&i*pr[j]<=maxn;j++){ 28 use[i*pr[j]]=1; 29 if(i%pr[j]) phi[i*pr[j]]=phi[i]*phi[pr[j]]; 30 else {phi[i*pr[j]]=phi[i]*pr[j];break;} 31 } 32 } 33 } 34 struct Seg{ 35 ull val[NN<<2],tag[NN<<2]; 36 void pushdown(int rt){ 37 if(!tag[rt]) return; 38 val[rt<<1]+=tag[rt],val[rt<<1|1]+=tag[rt]; 39 tag[rt<<1]+=tag[rt],tag[rt<<1|1]+=tag[rt]; 40 tag[rt]=0; 41 } 42 void update(int L,int R,int l,int r,int rt,int w) 43 { 44 if(L<=l&&r<=R){tag[rt]+=w,val[rt]+=w;return;} 45 int mid=(l+r)>>1;pushdown(rt); 46 if(L<=mid) update(L,R,l,mid,rt<<1,w); 47 if(R>mid) update(L,R,mid+1,r,rt<<1|1,w); 48 //pushup(rt); 49 } 50 ull query(int x,int l,int r,int rt) 51 { 52 if(l==r) return val[rt]; 53 int mid=(l+r)>>1;pushdown(rt); 54 if(x<=mid) return query(x,l,mid,rt<<1); 55 else return query(x,mid+1,r,rt<<1|1); 56 //pushup(rt); 57 } 58 }s; 59 ull qpow(ull x,ull y,int p,int &fl) 60 { 61 ull ans=1; 62 if(y==0) return 1; 63 if(x>=p) x%=p,fl=1; 64 while(y){ 65 if(y&1){ 66 if(ans*x>=p) 67 ans=ans*x%p,fl=1; 68 else ans=ans*x; 69 } 70 if(x*x>=p) 71 x=x*x%p,fl=1; 72 else x=x*x; 73 y>>=1; 74 }return ans; 75 } 76 ull euler(int i,int ma,int p,int &fl) 77 { 78 ull ai=s.query(i,1,n,1); 79 if(p==1) {fl=1;return 0;} 80 if(ai>=p) fl=1; 81 if(ai%p==0) return 0; 82 if(ai==1) {fl=0;return 1;} 83 int o=0;ai%=p; 84 if(i==ma) return ai; 85 ull ans=euler(i+1,ma,phi[p],o); 86 if(o==1) fl=1; 87 /*if(ans==0){ 88 if(o){ 89 if(ai%p==0) {return 0;} 90 else return qpow(ai,ans+phi[p],p,fl); 91 }else{ 92 return 1; 93 } 94 }else{ 95 if(o) return qpow(ai,ans+phi[p],p,fl); 96 else return qpow(ai,ans,p,fl); 97 }*/ 98 if(ans==0&&o&&ai==0) return 0; 99 else if(o) return qpow(ai,ans+phi[p],p,fl); 100 else return qpow(ai,ans,p,fl); 101 } 102 103 int main() 104 { 105 //freopen("t2.in","r",stdin); 106 //freopen("testdata.in","r",stdin); 107 //freopen("a.out","w",stdout); 108 int fl,x,y,z; 109 scanf("%d%d",&n,&m); 110 Pre(); 111 for(int i=1;i<=n;i++) 112 x=gint(),s.update(i,i,1,n,1,x); 113 for(int i=1;i<=m;i++) 114 { 115 fl=gint(),x=gint(),y=gint(),z=gint(); 116 if(fl==1){ 117 s.update(x,y,1,n,1,z); 118 }else{ 119 int fl=0; 120 printf("%llu\n",euler(x,y,z,fl)); 121 } 122 } 123 return 0; 124 }

?

轉載于:https://www.cnblogs.com/guapisolo/p/10001422.html

總結

以上是生活随笔為你收集整理的BZOJ 5394 [Ynoi2016]炸脖龙 (线段树+拓展欧拉定理)的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。