結果

問題 No.181 A↑↑N mod M
ユーザー glreto
提出日時 2021-03-20 17:20:09
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 1,828 bytes
コンパイル時間 809 ms
コンパイル使用メモリ 97,220 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-21 05:51:28
合計ジャッジ時間 2,056 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 6
other AC * 37
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include<vector>
#include<string>
#include<algorithm>
#include<iomanip>
#include<map>
#include<complex>
#include<functional>
#include<stack>
#include<queue>
#include<set>
using namespace std;
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define req(i,n) for(int i = 1;i <=  n; i++)
#define rrep(i,n) for(ll i = n-1;i >= 0;i--)
#define ALL(obj) begin(obj), end(obj)
#define RALL(a) rbegin(a),rend(a)
typedef long long ll;
typedef long double ld;
template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; }
template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; }
const ll INF = 1e18;
const int inf = 0x3fffffff;
typedef vector<ll> vec;
typedef vector<vec> mat;
const ll mod = 1e9 + 7;
template<typename T>
T totient(T n) {
    T res = n;
    for (T 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;
}
template<typename T>
T mpow(T a, T n, T m, T& f) {
    using ull = unsigned long long;
    if (a == 1 || n == 0) return 1;
    ull v = 1, x = a, z = m;
    f |= x >= z;x %= m;
    while (1) {
        if (n & 1) v *= x;
        if (v >= z) f = 1, v %= m;
        n >>= 1;
        if (!n) break;
        x = x * x;
        if (x >= z) f = 1, x %= m;
    }return v;
}template<typename T>
T tetration(T a, T n, T m, T& f) {
    if (a == 0) return ~n & 1;
    if (m == 1) return f = 1;
    if (a == 1 || n == 0) return 1;
    if (n == 1) {
        f |= a >= m;
        return a % m + f * m;
    }T z = tetration(a, n - 1, totient(m), f);
    T r = mpow(a, z, m, f);return r + f * m;
}
int main() {
    int a, n, m,f=0; cin >> a >> n >> m;
    cout << tetration(a, n, m,f)%m << endl;
}
0