結果

問題 No.2104 Multiply-Add
ユーザー 👑 KazunKazun
提出日時 2022-10-21 22:47:56
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 40 ms / 2,000 ms
コード長 1,468 bytes
コンパイル時間 150 ms
コンパイル使用メモリ 82,200 KB
実行使用メモリ 54,316 KB
最終ジャッジ日時 2024-07-01 07:13:05
合計ジャッジ時間 2,914 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 39 ms
52,652 KB
testcase_01 AC 38 ms
53,072 KB
testcase_02 AC 39 ms
53,852 KB
testcase_03 AC 39 ms
52,948 KB
testcase_04 AC 39 ms
53,048 KB
testcase_05 AC 38 ms
52,868 KB
testcase_06 AC 38 ms
53,532 KB
testcase_07 AC 38 ms
53,332 KB
testcase_08 AC 40 ms
54,316 KB
testcase_09 AC 39 ms
53,336 KB
testcase_10 AC 39 ms
52,936 KB
testcase_11 AC 39 ms
53,680 KB
testcase_12 AC 39 ms
53,368 KB
testcase_13 AC 40 ms
52,996 KB
testcase_14 AC 39 ms
53,436 KB
testcase_15 AC 39 ms
53,072 KB
testcase_16 AC 39 ms
52,576 KB
testcase_17 AC 39 ms
53,244 KB
testcase_18 AC 39 ms
53,020 KB
testcase_19 AC 40 ms
52,932 KB
testcase_20 AC 38 ms
53,436 KB
testcase_21 AC 38 ms
53,564 KB
testcase_22 AC 40 ms
52,872 KB
testcase_23 AC 39 ms
53,380 KB
testcase_24 AC 38 ms
53,264 KB
testcase_25 AC 39 ms
53,708 KB
testcase_26 AC 40 ms
52,948 KB
testcase_27 AC 40 ms
53,260 KB
testcase_28 AC 40 ms
53,264 KB
testcase_29 AC 39 ms
53,376 KB
testcase_30 AC 39 ms
53,624 KB
testcase_31 AC 39 ms
53,536 KB
testcase_32 AC 38 ms
53,744 KB
testcase_33 AC 39 ms
53,136 KB
testcase_34 AC 38 ms
53,800 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def operate(x,y,T):
        for t,k in T:
            if t==1:
                x+=k*y
            else:
                y+=k*x
        return x,y

def solve():
    from math import gcd

    a,b,c,d=map(int,input().split())
    if gcd(a,b)!=gcd(c,d):
        print(-1)
        return

    if gcd(a,b)==0:
        print(0)
        return

    def euclid(a,b):
        res=[]
        while a and b:
            if abs(a)>abs(b):
                k=abs(a)//abs(b)
                x=a+k*b; y=a-k*b
                if abs(x)<abs(y):
                    res.append((1,k))
                    a+=k*b
                else:
                    res.append((1,-k))
                    a-=k*b
            else:
                k=abs(b)//abs(a)
                x=b+k*a; y=b-k*a
                if abs(x)<abs(y):
                    res.append((2,k))
                    b+=k*a
                else:
                    res.append((2,-k))
                    b-=k*a
        return res

    T1=euclid(a,b); a2,b2=operate(a,b,T1)
    T2=euclid(c,d); c2,d2=operate(c,d,T2)

    T3=[]
    if a2==0:
        a2,b2=b2,a2
        T3=[(1,1),(2,-1)]

    T4=[]
    if c2==0:
        c2,d2=d2,c2
        T4=[(1,1),(2,-1)]

    T5=[]
    if a2*c2<0:
        T5=[(2,-2),(1,1),(2,-2)]

    T2=[(t,-k) for t,k in T2][::-1]
    T4=[(t,-k) for t,k in T4][::-1]

    T=T1+T3+T5+T4+T2

    print(len(T))
    for t,k in T:
        print(t,k)

#==================================================
solve()
0