結果

問題 No.181 A↑↑N mod M
コンテスト
ユーザー mamekin
提出日時 2015-04-12 13:49:23
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,690 bytes
コンパイル時間 655 ms
コンパイル使用メモリ 91,596 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-04 13:57:03
合計ジャッジ時間 2,108 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 6
other AC * 22 WA * 13 RE * 2
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <iostream>
#include <sstream>
#include <fstream>
#include <iomanip>
#include <algorithm>
#include <cmath>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <bitset>
#include <numeric>
#include <limits>
#include <climits>
#include <cfloat>
#include <functional>
using namespace std;

int main()
{
    int a, n, m;
    cin >> a >> n >> m;
    if(n == 0){
        cout << 1 << endl;
        return 0;
    }

    vector<int> p(m, -1);
    int p1, p2;
    long long x = 1;
    for(int i=0; ; ++i){
        if(p[x] != -1){
            p1 = p[x];
            p2 = i - p[x];
            break;
        }
        p[x] = i;

        x *= a;
        x %= m;
    }

    vector<int> q(m, -1);
    int q1, q2;
    long long y = 1;
    for(int i=0; ; ++i){
        if(q[y] != -1){
            q1 = q[y];
            q2 = i - q[y];
            break;
        }
        q[y] = i;

        long long y2 = 1;
        for(int j=0; j<y; ++j){
            y2 *= a;
            if(y2 > p1){
                y2 -= p1;
                y2 %= p2;
                y2 += p1;
            }
        }
        y = y2;
    }

    -- n;
    if(n > q1){
        n -= q1;
        n %= q2;
        n += q1;
    }
    long long z = 1;
    for(int i=0; i<n; ++i){
        long long z2 = 1;
        for(int j=0; j<z; ++j){
            z2 *= a;
            if(z2 > p1){
                z2 -= p1;
                z2 %= p2;
                z2 += p1;
            }
        }
        z = z2;
    }

    long long ans = 1;
    for(int j=0; j<z; ++j){
        ans *= a;
        ans %= m;
    }
    cout << ans << endl;

    return 0;
}
0