結果
問題 | No.2324 Two Countries within UEC |
ユーザー |
|
提出日時 | 2023-05-28 14:52:49 |
言語 | D (dmd 2.109.1) |
結果 |
AC
|
実行時間 | 222 ms / 2,000 ms |
コード長 | 1,316 bytes |
コンパイル時間 | 3,764 ms |
コンパイル使用メモリ | 173,096 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-26 23:23:52 |
合計ジャッジ時間 | 9,731 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 41 |
ソースコード
import std; void main () { long N, M, P, Q; { auto buf = readln.split.to!(long[]); N = buf[0], M = buf[1], P = buf[2], Q = buf[3]; } solve(N, M, P, Q); } void solve (long N, long M, long P, long Q) { foreach (_; 0..Q) { auto buf = readln.split.to!(long[]); long x = buf[0], f = buf[1]; if (x % P != 0) { // x^{-1}が存在 long inv_x = modinv(x, P); long modular = (f * inv_x) % P; // y == modular(MOD P) if (M < modular) { writeln(0); } else if (modular != 0) { writeln(1 + (M - modular)/P); } else { writeln((M - modular)/P); } } else { // x^{-1}が存在しない(任意のyについて親交度が0である) if (f == 0) { writeln(M); } else { writeln(0); } } } } long modpow (long x, long a, const long MOD) { if (a == 0) { return 1; } if (a == 1) { return x; } if (a % 2 == 1) { return (x * modpow(x, a-1, MOD)) % MOD; } long t = modpow(x, a/2, MOD); return (t * t) % MOD; } long modinv (long x, const long MOD) { return modpow(x, MOD - 2, MOD); }