結果
問題 | No.2183 LCA on Rational Tree |
ユーザー | shobonvip |
提出日時 | 2023-01-08 04:01:03 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 443 ms / 2,000 ms |
コード長 | 2,664 bytes |
コンパイル時間 | 706 ms |
コンパイル使用メモリ | 82,304 KB |
実行使用メモリ | 78,808 KB |
最終ジャッジ日時 | 2024-05-08 21:57:37 |
合計ジャッジ時間 | 3,000 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 59 ms
62,592 KB |
testcase_01 | AC | 443 ms
78,808 KB |
testcase_02 | AC | 407 ms
78,720 KB |
testcase_03 | AC | 420 ms
78,464 KB |
testcase_04 | AC | 202 ms
78,400 KB |
testcase_05 | AC | 416 ms
77,952 KB |
ソースコード
import os import sys from io import BytesIO, IOBase BUFSIZE = 8192 class FastIO(IOBase): newlines = 0 def __init__(self, file): self._fd = file.fileno() self.buffer = BytesIO() self.writable = "x" in file.mode or "r" not in file.mode self.write = self.buffer.write if self.writable else None def read(self): while True: b = os.read(self._fd, max(os.fstat(self._fd).st_size, BUFSIZE)) if not b: break ptr = self.buffer.tell() self.buffer.seek(0, 2), self.buffer.write(b), self.buffer.seek(ptr) self.newlines = 0 return self.buffer.read() def readline(self): while self.newlines == 0: b = os.read(self._fd, max(os.fstat(self._fd).st_size, BUFSIZE)) self.newlines = b.count(b"\n")+(not b) ptr = self.buffer.tell() self.buffer.seek(0, 2), self.buffer.write(b), self.buffer.seek(ptr) self.newlines -= 1 return self.buffer.readline() def flush(self): if self.writable: os.write(self._fd, self.buffer.getvalue()) self.buffer.truncate(0), self.buffer.seek(0) class IOWrapper(IOBase): def __init__(self, file): self.buffer = FastIO(file) self.flush = self.buffer.flush self.writable = self.buffer.writable self.write = lambda s:self.buffer.write(s.encode("ascii")) self.read = lambda:self.buffer.read().decode("ascii") self.readline = lambda:self.buffer.readline().decode("ascii") sys.stdin, sys.stdout = IOWrapper(sys.stdin), IOWrapper(sys.stdout) input = lambda:sys.stdin.readline().rstrip("\r\n") from math import gcd # plist を素数リストを使って高速化 def plist(m:int, ps:list) -> list: pf = [] for i in ps: flag = False while m%i == 0: flag = True m //= i if flag: pf.append(i) if m>1 : pf.append(m) return pf def gogo(p:int, q:int, lst:list): pl = plist(q-p, ps) while q-p > 1: x = 10 ** 18 for i in pl: if (q-p)%i == 0: x = min(x, i-(p%i)) g = gcd(p+x, q+x) lst.append((q-p, p, p+x)) p, q = (p+x)//g, (q+x)//g lst.append((1, p, 9 * 10 ** 18)) # 素数列挙 ps = [] er = [0] * 32001 for i in range(2, 32001): if er[i] == 0: ps.append(i) for j in range(2*i, 32001, i): er[j] = 1 # ここから本編 Q = int(input()) for _ in range(Q): lst1 = [] lst2 = [] p1,q1,p2,q2 = map(int,input().split()) assert p1 < q1 assert p2 < q2 assert gcd(p1, q1) == 1 assert gcd(p2, q2) == 1 gogo(p1,q1,lst1) gogo(p2,q2,lst2) g_begin = dict() g_end = dict() for i, st, ed in lst1: g_begin[i] = st g_end[i] = ed ans1 = 0 ans2 = 0 for i, st, ed in lst2: if i in g_begin.keys(): if max(0, min(ed, g_end[i]) - max(st, g_begin[i])) > 0: ans1 = max(st, g_begin[i]) ans2 = i + ans1 break print(ans1, ans2)