結果

問題 No.1936 Rational Approximation
ユーザー ripity
提出日時 2022-05-13 23:19:46
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
WA  
実行時間 -
コード長 503 bytes
コンパイル時間 270 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 11,008 KB
最終ジャッジ日時 2024-07-22 03:36:25
合計ジャッジ時間 1,391 ms
ジャッジサーバーID
(参考情報)
judge3 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample WA * 1
other WA * 14
権限があれば一括ダウンロードができます

ソースコード

diff #

from math import gcd

def modinv(a, b):
    p = b
    x, y, u, v = 1, 0, 0, 1
    while b:
        k = a // b
        x -= k * u
        y -= k * v
        x, u = u, x
        y, v = v, y
        a, b = b, a % b
    x %= p
    if x < 0:
        x += p
    return x

def calcR(p, q):
    inv = modinv(-p, q)
    G = gcd(inv-q, p)
    return (inv-q+p)/G

def calcL(p, q):
    inv = modinv(q, -p)
    G = gcd(inv-p, q)
    return (inv-p+q)/G

p, q = map(int, input().split())
print(calcL(p, q)+calcR(p, q))
0