結果

問題 No.2125 Inverse Sum
ユーザー daiotadaiota
提出日時 2022-11-19 10:21:21
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 131 ms / 2,000 ms
コード長 855 bytes
コンパイル時間 1,916 ms
コンパイル使用メモリ 174,552 KB
実行使用メモリ 11,780 KB
最終ジャッジ日時 2024-09-20 11:48:06
合計ジャッジ時間 3,699 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

#define REP(i,n) for(int i=0;i<int(n);i++)
typedef long long ll;
typedef pair<int,int> P;

vector<ll> div(ll n){
	vector<ll> d;
	for(ll i=1;i*i<=n;i++){
		if(n%i==0){
			d.push_back(i);
			if(i!=n/i) d.push_back(n/i);
		}
	}

	return d;
}

int main(void){
	ll i,j;

	cin.tie(0);  ios_base::sync_with_stdio(false);

	ll p,q;
	cin >> p >> q;

	vector<ll> a=div(q);
	vector<ll> b;
	ll s=a.size();
	for(i=0;i<s;i++){
		for(j=i;j<s;j++){
			b.push_back(a[i]*a[j]);
		}
	}

	sort(b.begin(),b.end());
	b.erase(unique(b.begin(),b.end()),b.end());
	ll t=b.size();

	vector<ll> n,m;
	ll T=0;
	for(i=0;i<t;i++){
		if((b[i]+q)%p==0 && (q*q/b[i]+q)%p==0){
			n.push_back((b[i]+q)/p);
			m.push_back((q*q/b[i]+q)/p);
			T++;
		}
	}

	cout << T << endl;
	REP(i,T){
		cout << n[i] << ' ' << m[i] << endl;
	}


	return 0;
}
0