結果

問題 No.109 N! mod M
コンテスト
ユーザー Mi_Sawa
提出日時 2015-04-08 22:55:32
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,838 bytes
コンパイル時間 1,277 ms
コンパイル使用メモリ 159,508 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-04 13:20:17
合計ジャッジ時間 3,148 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 8 WA * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define all(x) begin(x),end(x)
#define rall(x) (x).rbegin(),(x).rend()
#define REP(i,b,n) for(int i=(int)(b);i<(int)(n);++i)
#define rep(i,n) REP(i,0,n)
#define repsz(i,v) rep(i,(v).size())
#define aur auto&
#define bit(n) (1LL<<(n))
#define eb emplace_back
#define mt make_tuple
#define fst first
#define snd second
using namespace std;
typedef long long ll;
//#define int long long
template<class C>int size(const C &c){ return c.size(); }
template<class T>bool chmin(T&a,const T&b){if(a<=b)return false;a=b;return true;}
template<class T>bool chmax(T&a,const T&b){if(a>=b)return false;a=b;return true;}

bool is_prime(ll n){
    for(int p = 2; p * p <= n; ++p) if(n % p == 0) return false;
    return true;
}
ll exgcd(ll a, ll b, ll &x, ll &y){//{{{
    ll u[] = {a, 1, 0}, v[] = {b, 0, 1};
    while(*v){
        ll t = *u / *v;
        rep(i, 3) u[i] -= t * v[i];
        swap(u, v);
    }
    if(*u < 0) rep(i, 3) u[i] = -u[i];
    x = u[1]; y = u[2];
    return u[0];
}//}}}
ll invMod(ll a, ll m){//{{{
    ll p, q;
    exgcd(a, m, p, q);
    return (p%m+m)%m;
}//}}}
bool solve(){
    int T;
    cin >> T;
    rep(_, T){
        ll n, m;
        cin >> n >> m;
        if(n >= m){ cout << 0 << endl; continue; }
        if(n < 1000 * 1000){
            ll res = 1;
            rep(i, n+1) if(i) (res *= i) %= m;
            cout << res << endl;
            continue;
        }
        if(!is_prime(m)){ cout << 0 << endl; continue; }
        ll res = -1;
        for(int i = m-1; i > n; --i) (res *= invMod(i, m)) %= m;
        if(res < 0) res += m;
        cout << res << endl;
    }
    return true;
}
signed main(){
    cin.tie(nullptr);
    ios_base::sync_with_stdio(false);
    cout << std::fixed << std::setprecision(10);
    solve();
    return 0;
}
// vim:set foldmethod=marker commentstring=//%s:
0