結果

問題 No.109 N! mod M
コンテスト
ユーザー Hachimori
提出日時 2014-12-23 19:52:09
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 1,434 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 586 ms
コンパイル使用メモリ 62,460 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-06-12 04:12:46
合計ジャッジ時間 1,572 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 2 WA * 6 RE * 1
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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


int n, m;

void read() {
    cin >> n >> m;
}


vector<int> getFactor(int M) {
    vector<int> factor;
    for (int i = 2; i * i <= M; ++i)
        if (M % i == 0) {
            while (M % i == 0) {
                M /= i;
            }
            factor.push_back(i);
        }
    if (M > 1)
        factor.push_back(M);
    return factor;
}


int modpow(int p, int nMul, int mod) {
    int ret = 1;
    while (nMul) {
        if (nMul & 1) ret = 1LL * ret * p % mod;
        ret = 1LL * ret * ret % mod;
        nMul >>= 1;
    }
    return ret;
}


// mod: 素数
int inv(int num, int mod) {
    return modpow(num, mod - 2, mod);
}


void work() {
    if (n >= m) {
        cout << 0 << endl;
        return;
    }

    vector<int> factor = getFactor(m);
    
    if (factor.size() == 1) {
        int toInv = 1;
        for (int i = m - 1; i > n; --i)
            toInv = 1LL * toInv * i % m;
        cout << 1LL * (m - 1) * inv(toInv, m) % m << endl;
    }
    else if (n > factor.back()) {
        cout << 0 << endl;
    }
    else {
        int ans = 1;
        for (int i = 1; i <= n; ++i) {
            ans = 1LL * ans * i % m;
            if (ans == 0)
                break;
        }
        cout << ans << endl;
    }
}


int main() {
    int cases;
    cin >> cases;
    for (int i = 0; i < cases; ++i) {
        read();
        work();
    }
    return 0;
}
0