結果

問題 No.1232 2^x = x
コンテスト
ユーザー t33f
提出日時 2019-08-16 13:35:04
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 125 ms / 2,000 ms
コード長 1,065 bytes
コンパイル時間 587 ms
コンパイル使用メモリ 64,528 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-11-21 03:26:47
合計ジャッジ時間 1,357 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 3
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cassert>
#include <iostream>
using namespace std;

int power(long long e, int m) {
    assert (m % 2 == 1);
    long long ans = 1, p = 2;
    while (e > 0) {
        if (e & 1) ans = ans * p % m;
        p = p * p % m;
        e /= 2;
    }
    return ans;
}

int totient(int n) {
    int ans = 1;
    for (int i = 2; i * i <= n; i++) {
        int q = 1;
        while (n % i == 0) { n /= i; q *= i; }
        if (q > 1) ans *= q / i * (i-1);
    }
    if (n > 1) ans *= n-1;
    return ans;
}

int tower(int n, int m) {
    int r = 1, e = 0;
    while (m % 2 == 0) { r *= 2; e++; m /= 2; }
    if (m == 1) return 0;
    int t = totient(m);
    return power(tower(n-1, t) + t - e, m) * r;
}

long long solve(int p) {
    if (p == 2) return 2;
    long long a = tower(50, p), b = tower(50, p-1);
    if (a > b) b += p-1;
    return a + (b-a) * p;
}

int main() {
    int n; cin >> n;
    while (n--) {
        int p; cin >> p;
        long long ans = solve(p);
        assert (p == 2 || power(ans, p) % p == ans % p);
        cout << ans << endl;
    }
}
0