結果

問題 No.109 N! mod M
ユーザー なおなお
提出日時 2014-12-22 02:00:19
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 1,274 ms / 5,000 ms
コード長 1,741 bytes
コンパイル時間 1,209 ms
コンパイル使用メモリ 146,512 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-09-04 05:21:59
合計ジャッジ時間 3,179 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge12
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,376 KB
testcase_01 AC 77 ms
4,380 KB
testcase_02 AC 98 ms
4,376 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 18 ms
4,380 KB
testcase_05 AC 1,274 ms
4,380 KB
testcase_06 AC 5 ms
4,376 KB
testcase_07 AC 9 ms
4,380 KB
testcase_08 AC 2 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define REP(i, n)           for(int(i)=0;(i)<(n);++(i))
const int MOD = 1000000007;

const int max_r = 32000;
bool prime[max_r];

bool isprime(int n){
    if(n < max_r) return prime[n];
    if(n % 2 == 0) return false;
    for(int i = 3; i < max_r; i += 2){
        if(!prime[i]) continue;
        if(n % i == 0) return false;
    }
    return true;
}

ll extgcd(ll a,ll b,ll &m,ll &n){ll g=a;m=1;n=0;if(b)g=extgcd(b,a%b,n,m),n-=(a/b)*m;return g;}
ll divmod(ll n, ll m, ll mod){
    ll a,b;
    extgcd(m,mod,a,b);
    return (n * a) % mod;
}

ll solve(int N, int M){
    if(N >= M) return 0;    // N >= M の場合、N!の因数にMが含まれるため必ずMの倍数
    if(N == 0) return 1 % M;

    if(M < 200000){
        ll res = 1;
        for(int i = 2; i <= N; i++){
            res = (res * i) % M;
        }
        return res;
    }
    if(isprime(M)){
        // ウィルソンの定理より (M-1)! ≡ M-1 (mod M) なので、
        // N! * (N+1)(N+2)..(M-1) = (M-1)!
        // N! = M! / (N+1)(N+2)..(M-1)
        ll m = -1;
        for(int i = N+1; i < M; i++){
            m = divmod(m, i, M);
        }
        return (m+M)%M;
    }

    // M は 2e5以上の合成数
    // なんやかんやで 0 になる(´･_･`)

    return 0;
}


int T,N,M;
set<int> s;
int main(){
    memset(prime, 1, sizeof(prime));
    prime[0] = prime[1] = false;
    for(int i = 4; i < max_r; i += 2) prime[i] = false;
    for(int i = 3; i < max_r; i += 2){
        if(!prime[i]) continue;
        for(int j = i+i; j < max_r; j += i) prime[j] = false;
    }

    cin >> T;
    REP(i,T){
        cin >> N >> M;
        cout << solve(N,M) << endl;
    }
    return 0;
}
0