結果
問題 |
No.2558 中国剰余定理
|
ユーザー |
|
提出日時 | 2024-04-09 22:17:57 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 41 ms / 2,000 ms |
コード長 | 690 bytes |
コンパイル時間 | 262 ms |
コンパイル使用メモリ | 82,588 KB |
実行使用メモリ | 52,480 KB |
最終ジャッジ日時 | 2024-10-01 22:52:27 |
合計ジャッジ時間 | 2,510 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 29 |
ソースコード
import math def CRT(rem_pair, mod_pair): """中国剰余定理 - rem_pair = (r1, r2) - mod_pair = (m1, m2) - m1 != 0 and m2 != 0 x = m1 (mod r1) x = m2 (mod r2) となる整数xが存在する場合0以上の最小の自然数を返す(存在しない場合は-1) """ r1, r2 = rem_pair m1, m2 = mod_pair assert (m1 != 0 and m2 != 0) assert len(rem_pair) == len(mod_pair) == 2 g = math.gcd(m1, m2) if (r2 - r1) % g != 0: return -1 M1, M2, R = m1//g, (-m2)//g, (r2 - r1)//g inv = pow(M1, -1, M2) * R % M2 return (m1*inv + r1) % abs(m1*m2//g) A, B, a, b = map(int, input().split()) print(CRT((a, b), (A, B)))