結果
問題 |
No.187 中華風 (Hard)
|
ユーザー |
![]() |
提出日時 | 2024-10-24 21:13:52 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 173 ms / 3,000 ms |
コード長 | 920 bytes |
コンパイル時間 | 305 ms |
コンパイル使用メモリ | 82,704 KB |
実行使用メモリ | 76,504 KB |
最終ジャッジ日時 | 2024-10-24 21:13:57 |
合計ジャッジ時間 | 4,356 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 25 |
ソースコード
#yukicoder 187 中華風(Hard) 多倍長解法 #exec gcd #Reference: https://x.com/_miz_tom/status/1842086118990479854 exec('def gcd(x, y):\n' + ' x, y = abs(x), abs(y)\n' + ' if x == 0: return y\n' + ' if not (y := y % x): return x\n if not (x := x % y): return y\n' * 300 + ' return gcd(x, y)') def CRT(n1, M1, n2, M2): #smallest n ≡ n1 mod M1 ≡ n2 mod M2 if (n2 - n1) % ( G := gcd(M1, M2) ) != 0: return -1 return (n2 - n1) * pow(M1 // G, -1, M2 // G) % M2 * M1 // G + n1 #入力受取 N = int(input()) P = [tuple(map(int, input().split())) for _ in range(N)] MOD = 10 ** 9 + 7 #CRTを実行 rem, mod = 0, 1 for Xi, Yi in P: rem = CRT(rem, mod, Xi, Yi) if rem == -1: break mod = mod * Yi // gcd(mod, Yi) if rem == -1: #解なし print(-1) elif rem == 0: #正整数を探す → lcm(Yi) print(mod % MOD) else: #remを出力 print(rem % MOD)