結果

問題 No.2125 Inverse Sum
ユーザー 👑 potato167potato167
提出日時 2022-11-17 19:47:11
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 47 ms / 2,000 ms
コード長 888 bytes
コンパイル時間 2,112 ms
コンパイル使用メモリ 210,396 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-20 01:29:55
合計ジャッジ時間 3,302 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll=long long;
/*
N=an
M=am
1/N + 1/M = (n+m)/anm 
P/Q=(n+m)/anm
a=(n+m)Q/nmP
*/

//Nの正の約数を列挙する
vector<long long> Divisors(long long N){
    vector<long long> p,q;
    long long i=1,K=0;
    while(i*i<N){
        if(N%i==0){
            p.push_back(i);
            q.push_back(N/i);
            K++;
        }
        i++;
    }
    if(i*i==N) p.push_back(i);
    for(int i=K-1;i>=0;i--){
        p.push_back(q[i]);
    }
    return p;
}
int main() {
	ll P,Q;
	cin>>P>>Q;
	auto div=Divisors(Q);
	vector<pair<ll,ll>> ans;
	for(auto n:div) for(auto m:div){
		if(Q%(n*m)!=0) continue;
		if(__gcd(n,m)!=1) continue;
		ll X=(n+m)*Q;
		ll Y=n*m*P;
		if(X%Y!=0) continue;
		ll a=X/Y;
		ans.push_back({a*n,a*m});
	}
	sort(ans.begin(),ans.end());
	cout<<ans.size()<<"\n";
	for(auto x:ans) cout<<x.first<<" "<<x.second<<"\n";
}
0