結果
| 問題 |
No.181 A↑↑N mod M
|
| コンテスト | |
| ユーザー |
beet
|
| 提出日時 | 2018-11-14 21:24:12 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 2 ms / 5,000 ms |
| コード長 | 1,011 bytes |
| コンパイル時間 | 1,699 ms |
| コンパイル使用メモリ | 192,812 KB |
| 最終ジャッジ日時 | 2025-01-06 16:40:33 |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 6 |
| other | AC * 37 |
ソースコード
#include<bits/stdc++.h>
using namespace std;
using Int = long long;
template<typename T1,typename T2> inline void chmin(T1 &a,T2 b){if(a>b) a=b;}
template<typename T1,typename T2> inline void chmax(T1 &a,T2 b){if(a<b) a=b;}
//BEGIN CUT HERE
Int phi(Int n){
Int res=n;
for(Int i=2;i*i<=n;i++){
if(n%i==0){
res=res/i*(i-1);
for(;n%i==0;n/=i);
}
}
if(n!=1) res=res/n*(n-1);
return res;
}
Int mpow(Int x,Int n,Int m,Int &f){
if(x==1||n==0) return 1;
Int v=1;
f|=x>=m;
x%=m;
while(1){
if(n&1) v*=x;
if(v>=m) f=1,v%=m;
n>>=1;
if(!n) break;
x=x*x;
if(x>=m) f=1,x%=m;
}
return v;
}
Int mtetra(Int a,Int n,Int m,Int &f){
if(a==0) return n==0;
if(m==1) return f=1;
if(n==0) return 1;
if(n==1) return a+f*m;
Int z=mtetra(a,n-1,phi(m),f);
Int r=mpow(a,z,m,f);
return r+f*m;
}
//END CUT HERE
//INSERT ABOVE HERE
signed solve(){
Int a,n,m,f=0;
cin>>a>>n>>m;
cout<<mtetra(a,n,m,f)%m<<endl;
return 0;
}
signed main(){
solve();
}
beet