結果
問題 | 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);}