#974930 (Python3) No.186 中華風 (Easy)

• 問題: No.186 中華風 (Easy)
• ユーザー: shinnshinn
• 提出日時: 2024-04-20 10:22:05
• 言語: Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0)
• 結果: AC
• 実行時間: 35 ms / 2,000 ms
• コード長: 812 bytes
• コンパイル時間: 96 ms
• コンパイル使用メモリ: 12,672 KB
• 実行使用メモリ: 10,880 KB
• 最終ジャッジ日時: 2024-10-12 05:53:02
• 合計ジャッジ時間: 1,903 ms
• ジャッジサーバーID （参考情報）: judge3 / judge2

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	35 ms; 10,752 KB
testcase_01	AC	32 ms; 10,752 KB
testcase_02	AC	30 ms; 10,752 KB
testcase_03	AC	32 ms; 10,752 KB
testcase_04	AC	29 ms; 10,752 KB
testcase_05	AC	30 ms; 10,752 KB
testcase_06	AC	29 ms; 10,880 KB
testcase_07	AC	29 ms; 10,752 KB
testcase_08	AC	30 ms; 10,880 KB
testcase_09	AC	30 ms; 10,752 KB
testcase_10	AC	30 ms; 10,752 KB
testcase_11	AC	30 ms; 10,752 KB
testcase_12	AC	30 ms; 10,752 KB
testcase_13	AC	29 ms; 10,752 KB
testcase_14	AC	29 ms; 10,752 KB
testcase_15	AC	28 ms; 10,752 KB
testcase_16	AC	28 ms; 10,752 KB
testcase_17	AC	28 ms; 10,880 KB
testcase_18	AC	31 ms; 10,752 KB
testcase_19	AC	28 ms; 10,752 KB
testcase_20	AC	29 ms; 10,880 KB
testcase_21	AC	31 ms; 10,752 KB
testcase_22	AC	29 ms; 10,752 KB

ソースコード

def gcd(a, b):
    if b == 0:
        return a
    else:
        return gcd(b, a%b)

def extgcd(a, b):
    #拡張ユークリッド互除法
    #ax + by = gcd(a, b)となる(x, y)を求めます
    if b == 0:
        return [1, 0]

    q = a//b
    r = a%b
    
    s, t = extgcd(b, r)
    y = s - q*t
    return [t, y]

def CRT(rs, ms):
    R = 0
    M = 1
    for r, m in zip(rs, ms):
        g = gcd(M, m)
        
        # assert r%g == R%g
        if r%g != R%g:
            return [-1, -1]
        s, _ = extgcd(M, m)

        lcm = m*M // g
        R += (r - R)//g * M * s
        R %= lcm
        M = lcm
    return [R, M]

x = []
y = []
for i in range(3):
    X, Y = map(int, input().split())
    x.append(X)
    y.append(Y)

ans = CRT(x, y)
if ans[0] == 0:
    print(ans[1])
else:
    print(ans[0])