結果

問題 No.187 中華風 (Hard)
ユーザー kimiyukikimiyuki
提出日時 2016-12-29 23:38:25
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 820 ms / 3,000 ms
コード長 752 bytes
コンパイル時間 144 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 11,264 KB
最終ジャッジ日時 2024-05-09 06:13:12
合計ジャッジ時間 11,111 ms
ジャッジサーバーID
(参考情報)
judge4 / judge3
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 393 ms
11,136 KB
testcase_01 AC 391 ms
11,136 KB
testcase_02 AC 74 ms
11,008 KB
testcase_03 AC 75 ms
11,136 KB
testcase_04 AC 743 ms
11,136 KB
testcase_05 AC 735 ms
11,264 KB
testcase_06 AC 734 ms
11,264 KB
testcase_07 AC 732 ms
11,264 KB
testcase_08 AC 811 ms
11,264 KB
testcase_09 AC 820 ms
11,264 KB
testcase_10 AC 814 ms
11,136 KB
testcase_11 AC 757 ms
11,264 KB
testcase_12 AC 741 ms
11,136 KB
testcase_13 AC 35 ms
10,880 KB
testcase_14 AC 35 ms
10,880 KB
testcase_15 AC 88 ms
11,136 KB
testcase_16 AC 89 ms
11,008 KB
testcase_17 AC 30 ms
10,880 KB
testcase_18 AC 264 ms
11,008 KB
testcase_19 AC 30 ms
10,752 KB
testcase_20 AC 524 ms
11,008 KB
testcase_21 AC 29 ms
10,752 KB
testcase_22 AC 750 ms
11,264 KB
testcase_23 AC 31 ms
10,880 KB
testcase_24 AC 30 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