結果

問題 No.181 A↑↑N mod M
コンテスト
ユーザー ゴリポン先生
提出日時 2026-01-09 10:43:02
言語 D
(dmd 2.111.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 1,019 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 4,952 ms
コンパイル使用メモリ 195,068 KB
実行使用メモリ 7,852 KB
最終ジャッジ日時 2026-01-09 10:43:10
合計ジャッジ時間 7,315 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 6
other AC * 37
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

module main;
// // https://ei1333.github.io/library/math/number-theory/ より
import std;

// オイラーのファイ関数
long eulerPhi(long n)
{
	long ret = n;
	for (long i = 2; i * i <= n; i++) {
		if (n % i) continue;
		ret -= ret / i;
		while (n % i == 0) n /= i;
	}
	if (n > 1) ret -= ret / n;
	return ret;
}
long modPow(long x, long n, long MOD)
{
	long ret = 1;
	while (n > 0) {
		if (n & 1) ret = ret * x % MOD;	 // n の最下位bitが 1 ならば x^(2^i) をかける
		x = x * x % MOD;
		n >>= 1;	// n を1bit右にずらす
	}
	return ret;
}
long modTetration(long a, long b, long m)
{
	if (m == 1) return 0;
	if (a == 0) return !(b & 1);
	if (b == 0) return 1;
	if (b == 1) return a % m;
	if (b == 2) return modPow(a, a, m);
	auto phi = eulerPhi(m);
	auto tmp = modTetration(a, b - 1, phi);
	if (tmp == 0) tmp += phi;
	return modPow(a, tmp, m);
}

void main()
{
	// 入力
	long A, N, M;
	readln.chomp.formattedRead("%d %d %d", A, N, M);
	// 答えの計算と出力
	writeln(modTetration(A, N, M));
}
0