結果

問題 No.2125 Inverse Sum
ユーザー FromBooskaFromBooska
提出日時 2023-03-06 17:20:57
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 884 bytes
コンパイル時間 179 ms
コンパイル使用メモリ 81,856 KB
実行使用メモリ 53,816 KB
最終ジャッジ日時 2023-10-18 05:12:17
合計ジャッジ時間 4,226 ms
ジャッジサーバーID
(参考情報)
judge15 / judge12
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 39 ms
53,816 KB
testcase_01 AC 39 ms
53,696 KB
testcase_02 AC 39 ms
53,696 KB
testcase_03 AC 39 ms
53,696 KB
testcase_04 AC 39 ms
53,696 KB
testcase_05 TLE -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
権限があれば一括ダウンロードができます

ソースコード

diff #

# 方針変える
# 1/(3x) + 1/(3y) = 1/2の解き方と同じだった
# 変形してQM+QN=PMN
# PPMN - PQM - PQN = 0
# (PM-Q)*(PN-Q) - QQ=0
# (PM-Q)*(PN-Q) = Q**2

from math import gcd
P, Q = map(int, input().split())
g = gcd(P, Q)
P = P//g
Q = Q//g

def divisors(n):
    lower_divisors , upper_divisors = [], []
    i = 1
    while i*i <= n:
        if n % i == 0:
            lower_divisors.append(i)
            if i != n // i:
                upper_divisors.append(n//i)
        i += 1
    return lower_divisors + upper_divisors[::-1]
    
divs = divisors(Q**2)
ans_set = set()
#print('P', P, 'Q', Q, 'divs', divs)
for d in divs:
    d2 = (Q**2)//d
    if (Q+d)%P == 0 and (Q+d2)%P == 0:
        m = (Q+d)//P
        n = (Q+d2)//P
        ans_set.add((n, m))
        #print(d, d2, n, m)
ans_list = sorted(list(ans_set))
print(len(ans_list))
for n, m in ans_list:
    print(n, m)
0