标签：上帝 phi return 10000005 int ll varphi 用法 P4139

求

\[2^{2^{2^{2^{2^{...}}}}}mod\,p \]

\[p\leq 10^7 \]

---

显然硬干是不行的，那么考虑别的思路。设 \(f(p)\) 为原式模 \(p\) 的解，那么 \(f(p)=2^{f(\varphi(p))+\varphi(x)}\) ,递归可以求出上一项的值即可，边界是 \(\varphi(p)=1\) 时 \(f(p)=0\) ，需要预处理出 \(\varphi\) 的值。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
int prime[10000005],phi[10000005],tot;
bool not_prime[10000005];
int t,pp;


void getp(){
	phi[1]=1;
	for(int i=2;i<=10000000;++i){
		if(not_prime[i]==0){
			prime[++tot]=i;
			phi[i]=i-1;
		}
		for(int j=1;(prime[j]*i)<=10000000 && j<=tot;++j){
			not_prime[i*prime[j]]=1;
			if((i%prime[j])==0){
				phi[i*prime[j]]=phi[i]*prime[j];
				break;
			}
			phi[i*prime[j]]=phi[i]*phi[prime[j]];
		}
	}
}

ll Qpow(ll a,ll b,ll p){
	ll ret=1;
	while(b){
		if(b&1)
			ret=(ret*a)%p;
		b>>=1;
		a=(a*a)%p;
	}
	return ret;
}

ll f(ll x){
	if(phi[x]==1)
		return 0;
	ll q=f(phi[x]);
	return Qpow(2,q+phi[x],x);
}



int main(){
	getp();
	scanf("%d",&t);
	while(t--){
		scanf("%d",&pp);
		printf("%lld\n",f(pp));
	}
	return 0;
}

标签：上帝,phi,return,10000005,int,ll,varphi,用法,P4139
来源： https://www.cnblogs.com/zhouzizhe/p/16639845.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。