結果

問題 No.186 中華風 (Easy)
ユーザー keitel339keitel339
提出日時 2021-03-02 22:56:22
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 42 ms / 2,000 ms
コード長 1,122 bytes
コンパイル時間 155 ms
コンパイル使用メモリ 82,428 KB
実行使用メモリ 52,608 KB
最終ジャッジ日時 2024-10-03 02:12:03
合計ジャッジ時間 1,743 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 39 ms
52,352 KB
testcase_01 AC 38 ms
51,584 KB
testcase_02 AC 38 ms
51,968 KB
testcase_03 AC 37 ms
52,480 KB
testcase_04 AC 38 ms
51,712 KB
testcase_05 AC 37 ms
52,096 KB
testcase_06 AC 37 ms
52,096 KB
testcase_07 AC 36 ms
52,352 KB
testcase_08 AC 35 ms
52,224 KB
testcase_09 AC 36 ms
51,840 KB
testcase_10 AC 37 ms
52,352 KB
testcase_11 AC 42 ms
52,096 KB
testcase_12 AC 35 ms
51,712 KB
testcase_13 AC 37 ms
52,224 KB
testcase_14 AC 38 ms
52,096 KB
testcase_15 AC 42 ms
51,968 KB
testcase_16 AC 38 ms
52,224 KB
testcase_17 AC 36 ms
52,608 KB
testcase_18 AC 39 ms
51,968 KB
testcase_19 AC 37 ms
52,008 KB
testcase_20 AC 37 ms
52,224 KB
testcase_21 AC 37 ms
51,840 KB
testcase_22 AC 38 ms
51,968 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def gcd_ext(a, b):   
    '''
    1次不定方程式 ax+by=gcd(a, b) の解と最大公約数を返す。負の値には未対応
    x, y, gcd(a, b)
    '''
    if a < b:
        x, y, s = gcd_ext(b, a)
        return y, x, s
    s, xs, ys = a, 1, 0
    t, xt, yt = b, 0, 1

    while t > 0:
        n = s // t
        s, t = t, s % t
        xs, xt = xt, xs - n * xt
        ys, yt = yt, ys - n * yt
    return xs, ys, s


def crt(b, m):
    '''
    すべてのiについて x≡b[i](mod, m[i]) の論理積を満たす x と lcm(m[:]) を求める。
    解がない場合 -1, 0 を返す
    '''
    n = len(b)
    assert n == len(m)
    m0, b0 = 1, 0
    for i in range(n):
        m1, b1 = m[i], b[i] % m[i]
        p, _, d = gcd_ext(m0, m1) # m0 * p + m1 * q = gcd(m0, m1)
        if (b1 - b0) % d: # b0 ≡ b1 (mod gcd(m0, m1)) の不成立
            return -1, 0
        lcm = m0 * m1 // d
        b0 = (b0 + (b1 - b0) // d * p * m0) % lcm
        m0 = lcm
    return b0, lcm

import sys
xy = list(map(int, sys.stdin.read().split()))
x = xy[0::2]
y = xy[1::2]
ans, lcm = crt(x, y)
print(ans if ans else lcm)
0