結果

問題 No.186 中華風 (Easy)
ユーザー shinnshinnshinnshinn
提出日時 2024-04-20 10:20:02
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
WA  
実行時間 -
コード長 758 bytes
コンパイル時間 89 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 10,880 KB
最終ジャッジ日時 2024-10-12 05:50:32
合計ジャッジ時間 1,808 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
10,752 KB
testcase_01 AC 32 ms
10,880 KB
testcase_02 AC 32 ms
10,752 KB
testcase_03 AC 32 ms
10,752 KB
testcase_04 AC 30 ms
10,752 KB
testcase_05 AC 30 ms
10,752 KB
testcase_06 AC 30 ms
10,752 KB
testcase_07 AC 32 ms
10,752 KB
testcase_08 AC 30 ms
10,752 KB
testcase_09 AC 29 ms
10,752 KB
testcase_10 AC 29 ms
10,752 KB
testcase_11 AC 31 ms
10,752 KB
testcase_12 AC 30 ms
10,880 KB
testcase_13 AC 31 ms
10,752 KB
testcase_14 AC 30 ms
10,752 KB
testcase_15 AC 30 ms
10,752 KB
testcase_16 WA -
testcase_17 WA -
testcase_18 AC 29 ms
10,880 KB
testcase_19 AC 28 ms
10,752 KB
testcase_20 AC 30 ms
10,752 KB
testcase_21 AC 30 ms
10,880 KB
testcase_22 AC 30 ms
10,752 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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)

print(CRT(x, y)[0])
0