結果
問題 | No.109 N! mod M |
ユーザー |
![]() |
提出日時 | 2015-04-08 22:46:23 |
言語 | C++11 (gcc 13.3.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 903 bytes |
コンパイル時間 | 1,120 ms |
コンパイル使用メモリ | 158,072 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-07-04 13:20:13 |
合計ジャッジ時間 | 1,791 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 8 WA * 1 |
ソースコード
#include "bits/stdc++.h" using namespace std; bool isPrime(long long N){ for (int i = 2; i * i <= N; i++) { if (N % i == 0) return false; } return true; } long long powmod(long long a, int p, long long M){ if (p == 0) return 1; if (p % 2 == 1){ return a * powmod(a, p - 1, M) % M; } else{ long long c = powmod(a, p / 2, M); return c * c % M; } } long long calc(int N, int M){ if (N >= M) return 0; if (N <= 100005){ long long ret = 1; for (int i = 1; i <= N; i++) { ret *= i; ret %= M; } return ret; } else{ if (!isPrime(M)) return 0; long long ret = M - 1; long long div = 1; for (int i = M - 1; i > N; i--) { div *= i; div %= M; } div = powmod(div, M - 2, M); ret *= div; ret %= M; return ret; } } int main(){ int T; cin >> T; while (T > 0){ T--; long long N, M; cin >> N >> M; cout << calc(N, M) << endl; } cin >> T; }