結果

問題 No.187 中華風 (Hard)
ユーザー kimiyukikimiyuki
提出日時 2016-12-29 23:38:25
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 830 ms / 3,000 ms
コード長 752 bytes
コンパイル時間 263 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 11,264 KB
最終ジャッジ日時 2024-12-15 16:29:30
合計ジャッジ時間 11,569 ms
ジャッジサーバーID
(参考情報)
judge2 / judge5
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 404 ms
11,136 KB
testcase_01 AC 395 ms
11,264 KB
testcase_02 AC 78 ms
11,008 KB
testcase_03 AC 78 ms
11,136 KB
testcase_04 AC 763 ms
11,264 KB
testcase_05 AC 761 ms
11,136 KB
testcase_06 AC 763 ms
11,264 KB
testcase_07 AC 761 ms
11,136 KB
testcase_08 AC 830 ms
11,264 KB
testcase_09 AC 825 ms
11,264 KB
testcase_10 AC 830 ms
11,264 KB
testcase_11 AC 764 ms
11,136 KB
testcase_12 AC 763 ms
11,264 KB
testcase_13 AC 37 ms
11,008 KB
testcase_14 AC 37 ms
11,008 KB
testcase_15 AC 91 ms
11,136 KB
testcase_16 AC 93 ms
11,008 KB
testcase_17 AC 31 ms
10,752 KB
testcase_18 AC 269 ms
11,136 KB
testcase_19 AC 31 ms
10,880 KB
testcase_20 AC 531 ms
11,136 KB
testcase_21 AC 33 ms
10,880 KB
testcase_22 AC 762 ms
11,136 KB
testcase_23 AC 33 ms
10,752 KB
testcase_24 AC 33 ms
10,752 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#!/usr/bin/env python3
import functools
import math
lcm = lambda a, b: a * b // math.gcd(a, b)
def extgcd(a, b):
    if b == 0:
        return 1, 0
    else:
        na, nb = extgcd(b, a % b)
        return nb, na - a//b * nb
def modinv(a, n):
    return extgcd(a, n)[0] % n
def crt(eqn1, eqn2):
    x1, m1 = eqn1
    x2, m2 = eqn2
    d = math.gcd(m1, m2)
    x = x1 + (m1 // d) * (x2 - x1) * modinv(m1 // d, m2 // d)
    m = lcm(m1, m2)
    return x % m, m
n = int(input())
eqn = [ tuple(map(int, input().split())) for _ in range(n) ]
ans, _ = functools.reduce(crt, eqn)
if ans == 0:
    ans += functools.reduce(lcm, map(lambda it: it[1], eqn))
for x, y in eqn:
    if ans % y != x:
        ans = -1
        break
else:
    ans %= 10**9+7
print(ans)
0