結果
問題 |
No.2125 Inverse Sum
|
ユーザー |
![]() |
提出日時 | 2022-11-18 21:32:51 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 74 ms / 2,000 ms |
コード長 | 1,253 bytes |
コンパイル時間 | 1,885 ms |
コンパイル使用メモリ | 181,332 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-20 01:57:01 |
合計ジャッジ時間 | 3,181 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#include <bits/stdc++.h> using namespace std; int main(){ int P, Q; cin >> P >> Q; int Q2 = Q; long long X = (long long) Q * Q; vector<pair<long long, int>> F1; for (int i = 2; i * i <= Q2; i++){ if (Q2 % i == 0){ int cnt = 0; while (Q2 % i == 0){ cnt++; Q2 /= i; } F1.push_back(make_pair(i, cnt * 2)); } } if (Q2 > 1){ F1.push_back(make_pair(Q2, 2)); } int cnt = F1.size(); vector<long long> F = {1}; for (int i = 0; i < cnt; i++){ int sz = F.size(); vector<long long> F2; for (int j = 0; j < sz; j++){ long long x = F[j]; for (int k = 0; k <= F1[i].second; k++){ F2.push_back(x); if (k < F1[i].second){ x *= F1[i].first; } } } F = F2; } sort(F.begin(), F.end()); int cnt2 = F.size(); vector<pair<long long, long long>> ans; for (int i = 0; i < cnt2; i++){ long long x = F[i], y = X / F[i]; if ((x + Q) % P == 0 && (y + Q) % P == 0){ long long N = (x + Q) / P; long long M = (y + Q) / P; ans.push_back(make_pair(N, M)); } } int T = ans.size(); cout << T << endl; for (int i = 0; i < T; i++){ cout << ans[i].first << ' ' << ans[i].second << endl; } }