結果

問題 No.2125 Inverse Sum
ユーザー FromBooskaFromBooska
提出日時 2023-03-06 16:47:48
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,228 bytes
コンパイル時間 284 ms
コンパイル使用メモリ 82,312 KB
実行使用メモリ 130,448 KB
最終ジャッジ日時 2024-09-18 01:52:09
合計ジャッジ時間 2,788 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 36 ms
52,908 KB
testcase_01 AC 43 ms
61,360 KB
testcase_02 AC 43 ms
60,512 KB
testcase_03 AC 43 ms
59,940 KB
testcase_04 AC 46 ms
61,172 KB
testcase_05 WA -
testcase_06 AC 59 ms
69,760 KB
testcase_07 WA -
testcase_08 AC 54 ms
65,720 KB
testcase_09 AC 39 ms
52,428 KB
testcase_10 AC 62 ms
72,976 KB
testcase_11 AC 43 ms
60,068 KB
testcase_12 AC 42 ms
60,460 KB
testcase_13 AC 51 ms
64,028 KB
testcase_14 AC 52 ms
70,356 KB
testcase_15 AC 44 ms
63,248 KB
testcase_16 AC 58 ms
74,060 KB
testcase_17 WA -
testcase_18 AC 45 ms
64,068 KB
testcase_19 AC 45 ms
64,028 KB
testcase_20 AC 63 ms
74,756 KB
testcase_21 WA -
testcase_22 AC 44 ms
63,944 KB
testcase_23 AC 48 ms
64,328 KB
testcase_24 WA -
testcase_25 AC 43 ms
61,804 KB
testcase_26 AC 35 ms
53,268 KB
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
testcase_30 AC 36 ms
53,204 KB
testcase_31 WA -
testcase_32 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

# 変形してQ*(M+N) = N*M*P
# gcd(P, Q)で両辺を割る
# もしP=Q=1の場合、M=N=2
# それ以外の場合、N,MはQの素因数を持つ、N+MはPの倍数
# Qの約数でNを全探索すればいい、約数だけではなく約数の倍数もありうるか

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]

from math import gcd
P, Q = map(int, input().split())
g = gcd(P, Q)
P = P//g
Q = Q//g
if P == 1:
    print(1)
    print(Q*2, Q*2)
else:
    ans_list = []
    divs = divisors(Q)
    visited = set()
    for d in divs:
        for k in range(1, 1000):
            n = d*k
            if n in visited:
                continue
            visited.add(n)
            
            if (P*n-Q) != 0 and (n*Q)%(P*n-Q) == 0 and (n*Q)//(P*n-Q) > 0:
                if Q*(((n*Q)//(P*n-Q))+n) == P*((n*Q)//(P*n-Q))*n:
                    ans_list.append((n, (n*Q)//(P*n-Q)))
                    #print(d, k, n)

    print(len(ans_list))
    for n, m in ans_list:
        print(n, m)
0