結果
問題 | No.109 N! mod M |
ユーザー |
|
提出日時 | 2017-11-12 21:27:56 |
言語 | D (dmd 2.109.1) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,107 bytes |
コンパイル時間 | 820 ms |
コンパイル使用メモリ | 117,516 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-12 22:29:51 |
合計ジャッジ時間 | 4,755 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 3 WA * 6 |
ソースコード
import std.stdio, std.array, std.string, std.conv, std.algorithm;import std.typecons, std.range, std.random, std.math, std.container;import std.numeric, std.bigint, core.bitop;bool is_prime(long n) {for (long i = 2; i * i <= n; ++i)if (n % i == 0)return false;return true;}long powmod(long a, long x, long m) {long ret = 1;while (x) {if (x % 2) ret = ret * a % m;a = a * a % m;x /= 2;}return ret;}long solve() {auto s = readln.split.map!(to!long);auto N = s[0];auto M = s[1];if (M == 1)return 0;if (N == 0)return 1;if (N >= M)return 0;if (is_prime(M)) {long ret = M - 1;for (long cnt = M - 1; cnt > N; --cnt)ret = powmod(ret, M-2, M);return ret;} else if (M < 10^^6) {long ret = 1;for (long i = 2; i <= N; ++i)ret = ret * i % M;return ret;} else {return 0;}}void main() {auto T = readln.chomp.to!int;while (T--)solve.writeln;}