結果
| 問題 | No.109 N! mod M |
| コンテスト | |
| ユーザー |
古寺いろは
|
| 提出日時 | 2015-04-08 22:38:44 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0 + boost 1.89.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 854 bytes |
| 記録 | |
| コンパイル時間 | 1,369 ms |
| コンパイル使用メモリ | 158,308 KB |
| 実行使用メモリ | 6,948 KB |
| 最終ジャッジ日時 | 2024-07-04 13:19:58 |
| 合計ジャッジ時間 | 1,918 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 2 WA * 7 |
ソースコード
#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 <= 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 - 2; i > N; i--)
{
div *= i;
}
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;
}
}
古寺いろは