結果
| 問題 |
No.181 A↑↑N mod M
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2016-12-28 16:44:53 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 2 ms / 5,000 ms |
| コード長 | 925 bytes |
| コンパイル時間 | 1,520 ms |
| コンパイル使用メモリ | 169,476 KB |
| 実行使用メモリ | 6,824 KB |
| 最終ジャッジ日時 | 2024-12-15 07:51:38 |
| 合計ジャッジ時間 | 2,726 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 6 |
| other | AC * 37 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
int phi[2001];
void sievePhi(int n) {
for (int i = 0; i < n; i++) {
phi[i] = i;
}
for (int i = 2; i < n; i++) {
if (phi[i] == i) {
for (int j = i; j < n; j += i) {
phi[j] = phi[j] / i * (i - 1);
}
}
}
}
int modpow(int a, long long b, int mod) {
if (b == 0) return 1;
return 1LL * modpow(1LL * a * a % mod, b / 2, mod) * (b & 1 ? a : 1) % mod;
}
int tetration(int a, int n) {
if (a == 1 || n == 0) return 1;
if (n == 1) return a;
if (n == 2) return a * a;
if (n == 3 && a == 2) return 16;
return 123456;
}
int modtet(int a, int n, int m) {
if (m == 1) return 0;
if (n == 0) return 1;
if (tetration(a, n - 1) < phi[m]) {
return modpow(a, modtet(a, n - 1, phi[m]), m);
} else {
return modpow(a, modtet(a, n - 1, phi[m]) + phi[m], m);
}
}
int main() {
int a, n, m;
cin >> a >> n >> m;
sievePhi(m + 1);
cout << modtet(a, n, m) << endl;
}