結果

問題 No.2183 LCA on Rational Tree
ユーザー vwxyzvwxyz
提出日時 2024-06-02 04:46:53
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 1,367 bytes
コンパイル時間 144 ms
コンパイル使用メモリ 82,392 KB
実行使用メモリ 78,172 KB
最終ジャッジ日時 2024-06-02 04:47:04
合計ジャッジ時間 8,499 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 36 ms
58,496 KB
testcase_01 AC 1,805 ms
77,440 KB
testcase_02 TLE -
testcase_03 AC 1,124 ms
77,320 KB
testcase_04 AC 355 ms
77,396 KB
testcase_05 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import math
import pypyjit
pypyjit.set_param('max_unroll_recursion=-1')

def Divisors(N):
    divisors=[]
    for i in range(1,N+1):
        if i**2>=N:
            break
        elif N%i==0:
            divisors.append(i)
    if i**2==N:
        divisors+=[i]+[N//i for i in divisors[::-1]]
    else:
        divisors+=[N//i for i in divisors[::-1]]
    return divisors

Q=int(input())
for q in range(Q):
    a,b,c,d=map(int,input().split())
    M=10**18
    def tree(p,q):
        D=Divisors(q-p)
        def f(p,q):
            if q-p==1:
                return [(1,p,M)]
            lst=[]
            for d in D[1:]:
                if (q-p)%d==0:
                    lst.append((-p)%d)
            mi=min(lst)
            pp=p+mi
            qq=q+mi
            g=math.gcd(pp,qq)
            return [(q-p,p,pp)]+f(pp//g,qq//g)
        return f(p,q)
    tree0=tree(a,b)
    tree1=tree(c,d)
    assert len(tree0)<=30
    assert len(tree1)<=30
    D,A=M,M
    tree0.sort()
    tree1.sort()
    while tree0 and tree1:
        if tree0[-1][0]<tree1[-1][0]:
            tree1.pop()
        elif tree0[-1][0]>tree1[-1][0]:
            tree0.pop()
        else:
            d0,a0,b0=tree0.pop()
            d1,a1,b1=tree1.pop()
            a=max(a0,a1)
            b=min(b0,b1)
            d=d0
            if a<=b:
                D,A=min((D,A),(d,a))
    print(A,A+D)
0