結果

問題 No.2125 Inverse Sum
コンテスト
ユーザー LyricalMaestro
提出日時 2024-12-29 13:32:57
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,470 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 217 ms
コンパイル使用メモリ 82,604 KB
実行使用メモリ 62,020 KB
最終ジャッジ日時 2024-12-29 13:33:00
合計ジャッジ時間 3,060 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2 WA * 1
other AC * 17 WA * 13
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

## https://yukicoder.me/problems/no/2125

import math

def calc_gcd(A, B):
    """
    正の整数A, Bの最大公約数を計算する
    """
    a = max(A, B)
    b = min(A, B)
    while a % b > 0:
        c = a % b
        a = b
        b = c
    return b

def main():
    P, Q = map(int, input().split())

    # P, Qを互いに素な整数として解く
    gcd = calc_gcd(P, Q)
    P //= gcd
    Q //= gcd

    # Qの約数分解
    q_divisors = []
    sqrt_q = int(math.sqrt(Q))
    for p in range(1, sqrt_q + 1):
        if Q % p == 0:
            q = Q // p
            q_divisors.append(p)
            if q != p:
                q_divisors.append(q)
    
    answers = []
    # nが取りうる値で全探索
    q_divisors_set = set(q_divisors)
    for q in q_divisors:
        l = Q - q * P
        r = -Q * q
        if abs(l) > 0 and abs(r) % abs(l) == 0:
            if r * l > 0:
                m = abs(r) // abs(l)
                if m >= 1:
                    answers.append((q, m))
        
        # 2倍したものも考慮
        q *= 2
        if q not in q_divisors_set:
            l = Q - q * P
            r = -Q * q
            if abs(l) > 0 and abs(r) % abs(l) == 0:
                if r * l > 0:
                    m = abs(r) // abs(l)
                    if m >= 1:
                        answers.append((q, m))

    print(len(answers))
    for n, m in answers:
        print(n, m)









if __name__ == "__main__":
    main()
0