結果

問題 No.2104 Multiply-Add
ユーザー 👑 KazunKazun
提出日時 2022-10-21 22:47:56
言語 PyPy3
(7.3.13)
結果
AC  
実行時間 70 ms / 2,000 ms
コード長 1,468 bytes
コンパイル時間 813 ms
コンパイル使用メモリ 86,656 KB
実行使用メモリ 71,724 KB
最終ジャッジ日時 2023-09-13 23:01:37
合計ジャッジ時間 4,952 ms
ジャッジサーバーID
(参考情報)
judge11 / judge15
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 68 ms
71,596 KB
testcase_01 AC 67 ms
71,188 KB
testcase_02 AC 69 ms
71,276 KB
testcase_03 AC 68 ms
71,312 KB
testcase_04 AC 69 ms
71,212 KB
testcase_05 AC 69 ms
71,196 KB
testcase_06 AC 66 ms
71,224 KB
testcase_07 AC 69 ms
71,596 KB
testcase_08 AC 65 ms
71,308 KB
testcase_09 AC 68 ms
71,460 KB
testcase_10 AC 67 ms
71,188 KB
testcase_11 AC 65 ms
71,324 KB
testcase_12 AC 66 ms
71,332 KB
testcase_13 AC 67 ms
71,312 KB
testcase_14 AC 64 ms
71,724 KB
testcase_15 AC 64 ms
71,452 KB
testcase_16 AC 65 ms
71,408 KB
testcase_17 AC 70 ms
71,156 KB
testcase_18 AC 66 ms
71,460 KB
testcase_19 AC 64 ms
71,376 KB
testcase_20 AC 63 ms
71,524 KB
testcase_21 AC 64 ms
71,552 KB
testcase_22 AC 63 ms
71,588 KB
testcase_23 AC 65 ms
71,196 KB
testcase_24 AC 63 ms
71,184 KB
testcase_25 AC 65 ms
71,672 KB
testcase_26 AC 63 ms
71,524 KB
testcase_27 AC 64 ms
71,552 KB
testcase_28 AC 63 ms
71,508 KB
testcase_29 AC 63 ms
71,512 KB
testcase_30 AC 63 ms
71,296 KB
testcase_31 AC 63 ms
71,580 KB
testcase_32 AC 65 ms
71,408 KB
testcase_33 AC 64 ms
71,508 KB
testcase_34 AC 63 ms
71,404 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