結果

問題 No.1936 Rational Approximation
ユーザー ripityripity
提出日時 2022-05-13 23:19:00
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
RE  
実行時間 -
コード長 503 bytes
コンパイル時間 92 ms
コンパイル使用メモリ 10,988 KB
実行使用メモリ 8,096 KB
最終ジャッジ日時 2023-09-29 08:58:49
合計ジャッジ時間 1,152 ms
ジャッジサーバーID
(参考情報)
judge15 / judge12
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
権限があれば一括ダウンロードができます

ソースコード

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