結果
| 問題 |
No.187 中華風 (Hard)
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-04-02 21:25:59 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,832 bytes |
| コンパイル時間 | 288 ms |
| コンパイル使用メモリ | 82,456 KB |
| 実行使用メモリ | 76,288 KB |
| 最終ジャッジ日時 | 2024-09-30 23:19:53 |
| 合計ジャッジ時間 | 2,974 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 23 WA * 2 |
ソースコード
import sys
import math
import bisect
from heapq import heapify, heappop, heappush
from collections import deque, defaultdict, Counter
from functools import lru_cache
from itertools import accumulate, combinations, permutations, product
sys.setrecursionlimit(1000000)
MOD = 10 ** 9 + 7
MOD99 = 998244353
input = lambda: sys.stdin.readline().strip()
NI = lambda: int(input())
NMI = lambda: map(int, input().split())
NLI = lambda: list(NMI())
SI = lambda: input()
SMI = lambda: input().split()
SLI = lambda: list(SMI())
EI = lambda m: [NLI() for _ in range(m)]
def inv_gcd(a, b):
a = a % b
if a == 0:
return b, 0
s = b
t = a
m0 = 0
m1 = 1
while t:
u = s // t
s -= t * u
m0 -= m1 * u
s, t = t, s
m0, m1 = m1, m0
if m0 < 0:
m0 += b // s
return s, m0
def inv_mod(x, m):
assert 1 <= m
z = inv_gcd(x, m)
assert z[0] == 1
return z[1]
def crt(r,m):
# r: 余りのlist
# m: modのlist
assert len(r) == len(m)
n = len(r)
r0 = 0
m0 = 1
for i in range(n):
assert 1 <= m[i]
r1 = r[i] % m[i]
m1 = m[i]
if m0 < m1:
r0, r1 = r1, r0
m0, m1 = m1, m0
if m0 % m1 == 0:
if r0 % m1 != r1:
return 0, 0
continue
g, im = inv_gcd(m0, m1)
u1 = m1 // g
if (r1 - r0) % g:
return 0, 0
x = (r1-r0) // g % u1 * im % u1
r0 += x * m0
m0 *= u1
if r0 < 0:
r0 += m0
return r0,m0
def main():
N = NI()
R = []
M = []
for _ in range(N):
x, y = NMI()
R.append(x)
M.append(y)
r, m = crt(R, M)
if r == 0:
print(-1)
else:
print(r % MOD)
if __name__ == "__main__":
main()