D. Power Tower (欧拉降幂&幂塔函数)

最新推荐文章于 2022-01-05 18:31:30 发布

酷酷的Herio

最新推荐文章于 2022-01-05 18:31:30 发布

阅读量1k

点赞数 2

分类专栏：数论

本文链接：https://blog.csdn.net/weixin_45750972/article/details/113548525

版权

数论专栏收录该内容

105 篇文章 2 订阅

订阅专栏

D. Power Tower (欧拉降幂&幂塔函数)

题意

求给定区间的数的幂塔值。

思路

扩展欧拉定理降幂的裸题，主要是扩欧+重新定义取模，欧拉函数可以 $\sqrt n$ +记忆化，主要是递归的思想，当递归到 $p h i (x) = 1$ 直接返回，因为任何数模 $1$ 为0，

则 $(sum)\% phi(x)+phi(x)=1$ 。

因此时间复杂度： $O(qlog^2n)$

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+5,M=2e4+5,inf=0x3f3f3f3f,mod=1e9+7;
#define mst(a,b) memset(a,b,sizeof a)
#define PII pair<int,int>
#define fi first
#define se second
#define pb push_back
//重新定义取模 
#define mod(a,b) a<b?a:a%b+b
int n,m;
int a[N];
unordered_map<int,int>mp;
ll ksm(ll a,ll n,ll m){
   ll ans=1;
   while(n){
   	if(n&1) ans=mod(ans*a,m);
   	a=mod(a*a,m);
   	n>>=1;
   }return ans;
}
int phi(int n){		//求phi(n)
   if(mp[n]) return mp[n];
   int s=n,nn=n;
   for(int i=2;i*i<=n;i++){
   	if(n%i==0){
   		s-=s/i;
   		while(n%i==0) n/=i;
   	}
   }
   if(n>1) s-=s/n;return mp[nn]=s;
}
ll solve(int l,int r,int m){
   if(l==r||m==1) return mod(a[l],m);
   return ksm(a[l],solve(l+1,r,phi(m)),m);
} 
int main(){
   scanf("%d%d",&n,&m);for(int i=1;i<=n;i++) scanf("%d",&a[i]);
   int q;scanf("%d",&q);while(q--){
   	int l,r;scanf("%d%d",&l,&r);
   	printf("%lld\n",solve(l,r,m)%m);
   } 
   return 0;
}