結果

問題 No.2125 Inverse Sum
ユーザー eve__fuyuki
提出日時 2024-05-07 14:48:51
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 304 ms / 2,000 ms
コード長 1,482 bytes
コンパイル時間 2,510 ms
コンパイル使用メモリ 210,456 KB
最終ジャッジ日時 2025-02-21 11:35:54
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;
void fast_io() {
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);
}

int main() {
    fast_io();
    long long p, q;
    cin >> p >> q;
    long long g = gcd(p, q);
    p /= g;
    q /= g;
    map<long long, int> fac;
    {
        long long t = q;
        for (long long i = 2; i * i <= t; i++) {
            while (t % i == 0) {
                fac[i] += 2;
                t /= i;
            }
        }
        if (t > 1) {
            fac[t] += 2;
        }
    }
    vector<long long> divs;
    divs.push_back(1);
    for (auto &i : fac) {
        int sz = divs.size();
        for (int j = 0; j < sz; j++) {
            long long cur = 1;
            for (int k = 0; k < i.second; k++) {
                cur *= i.first;
                divs.push_back(divs[j] * cur);
                divs.push_back(-divs[j] * cur);
            }
        }
    }
    vector<pair<long long, long long>> ans;
    for (long long d : divs) {
        if (d + q <= 0 || (d + q) % p != 0) {
            continue;
        }
        long long n = (d + q) / p;
        long long e = q * q / d;
        if (e + q <= 0 || (e + q) % p != 0) {
            continue;
        }
        long long m = (e + q) / p;
        ans.push_back({n, m});
    }
    sort(ans.begin(), ans.end());
    ans.erase(unique(ans.begin(), ans.end()), ans.end());
    cout << ans.size() << "\n";
    for (auto [n, m] : ans) {
        cout << n << " " << m << "\n";
    }
}
0