[bzoj3884]上帝与集合的正确用法-CSDN博客

题意：求${2^{{2^{{2^{...}}}}}}\bmod p$

解题关键：

因为${a^n} \equiv {a^{n\bmod \varphi (p) + \varphi (p)}}\bmod p,n > \varphi (p)$

所以，

$\begin{array}{l}
f(p) = {2^{{2^{{2^{...}}}}}}\bmod p = {2^{{2^{{2^{...}}}}\bmod \varphi (p) + \varphi (p)}}\bmod p\\
= {2^{f(\varphi (p)) + \varphi (p)}}\bmod p
\end{array}$

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 #include<cstdlib>
 5 #include<iostream>
 6 #include<cmath>
 7 using namespace std;
 8 typedef long long ll;
 9 const int MAXN=10000000+10;
10 int prime[MAXN];//保存素数
11 bool vis[MAXN];//初始化
12 int phi[MAXN];//欧拉函数
13 void Lphisieve(int n){
14     int cnt=0;
15     for(int i=2;i<n;i++){
16         if(!vis[i]){
17             prime[cnt++]=i;
18             phi[i]=i-1;
19         }
20         for(int j=0;j<cnt&&i*prime[j]<n;j++){
21             int k=i*prime[j];
22             vis[k]=true;
23             if(i%prime[j]==0){
24                 phi[k]=phi[i]*prime[j];
25                 break;
26             }
27             else   phi[k]=phi[i]*(prime[j]-1);
28         }
29     }
30 }
31 
32 ll mod_pow(ll x,ll n,ll p){
33     ll res=1;
34     while(n){
35         if(n&1)   res=res*x%p;
36         x=x*x%p;
37         n>>=1;
38     }
39     return res;
40 }
41 
42 ll work(ll n){
43     if(n==1) return 0;
44     return mod_pow(2,work(phi[n])+phi[n],n);
45 }
46 
47 int main(){
48     int T;
49     Lphisieve(10000001);
50     scanf("%d", &T);
51     while(T--){
52         ll n;
53         scanf("%lld",&n);
54         printf("%lld\n",work(n));
55     }
56     return 0;
57 }