結果

問題 No.186 中華風 (Easy)
ユーザー Kevinrobot34
提出日時 2021-06-26 07:37:25
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
WA  
実行時間 -
コード長 1,180 bytes
コンパイル時間 102 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 10,880 KB
最終ジャッジ日時 2024-06-25 07:51:55
合計ジャッジ時間 1,553 ms
ジャッジサーバーID
(参考情報)
judge4 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 8 WA * 15
権限があれば一括ダウンロードができます

ソースコード

diff #

# return (g, x, y) which satisfies ax + by = g and g = GCD(a,b)
def ext_gcd(a: int, b: int) -> tuple[int, int, int]:
    if b:
        g, y, x = ext_gcd(b, a % b)
        y -= (a // b) * x
        return g, x, y
    return a, 1, 0


# return (r, m) which satisfies x \equiv r (mod m) as solution
# when no solution ( (b1-b2) % gcd(m1, m2) != 0 ), return(0, -1)
def crt(m1: int, m2: int, b1: int, b2: int) -> tuple[int, int]:
    g, p, q = ext_gcd(m1, m2)
    if (b2 - b1) % g != 0:
        return 0, -1
    m = m1 * m2 // g
    s = (b2 - b1) // g
    r = b1 + s * m1 * p
    r %= m
    return r, m



# return (r, m) which satisfies x \equiv r (mod m) as solution
# when no solution ( (b1-b2) % gcd(m1, m2) != 0 ), return(0, -1)
def crt_n(ms: list[int], bs: list[int]) -> tuple[int, int]:
    r, m = 0, 1

    for bi, mi in zip(bs, ms):
        g, p, q = ext_gcd(m, mi)
        if (bi - r) % g != 0:
            return 0, -1
        s = (bi - r) // g * p % (mi // g)
        r += s * m
        m *= mi // g
        r %= m
    return r, m

bs = []
ms = []
for i in range(3):
    bi, mi = map(int, input().split())
    bs.append(bi)
    ms.append(mi)

r, m = crt_n(ms, bs)
print(m)
0