結果

問題 No.2183 LCA on Rational Tree
ユーザー tassei903
提出日時 2023-01-14 16:26:54
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 489 ms / 2,000 ms
コード長 1,516 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 605 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 77,952 KB
最終ジャッジ日時 2024-12-26 01:32:58
合計ジャッジ時間 3,734 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 6
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import sys
input = lambda :sys.stdin.readline()[:-1]
ni = lambda :int(input())
na = lambda :list(map(int,input().split()))
yes = lambda :print("yes");Yes = lambda :print("Yes");YES = lambda : print("YES")
no = lambda :print("no");No = lambda :print("No");NO = lambda : print("NO")
#######################################################################

pmax = 31625

p = [-1] * pmax
primes = []
for i in range(2, pmax):
    if p[i] == -1:
        primes.append(i)
        for j in range(i, pmax, i):
            p[j] = i

def fac(x):
    ans = []
    for p in primes:
        if x % p == 0:
            ans.append(p)
            while x % p == 0:
                x //= p
    if x != 1:
        return ans + [x]
    return ans

def nxt(p, q, f):
    x = 10**18
    for i in f:
        if (q-p) % i == 0:
            x = min(x, (p + i-1)//i*i-p)
    g = gcd(p+x,q+x)
    return (p+x)//g, (q+x)//g, x
from math import gcd
for _ in range(ni()):
    p1, q1, p2, q2 = na()
    f1 = fac(q1-p1)
    f2 = fac(q2-p2)
    while True:
        #print(p1,q1,p2,q2)
        if q1-p1 == q2-p2:
            if q1 > q2:
                q1,q2 = q2,q1
                p1,p2 = p2,p1
                f1,f2 = f2,f1
            p3,q3,x = nxt(p1, q1, f1)
            if q1 + x >= q2:
                print(p2, q2)
                break
            p1, q1 = p3, q3
        elif q1 - p1 > q2 - p2:
            q1,q2 = q2,q1
            p1,p2 = p2,p1
            f1,f2 = f2,f1
        p3, q3, x = nxt(p2, q2, f2)
        p2, q2 = p3, q3
0