結果
問題 | No.2125 Inverse Sum |
ユーザー |
|
提出日時 | 2022-11-18 22:03:47 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 17 ms / 2,000 ms |
コード長 | 1,314 bytes |
コンパイル時間 | 2,222 ms |
コンパイル使用メモリ | 207,372 KB |
最終ジャッジ日時 | 2025-02-08 21:44:46 |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#include<bits/stdc++.h>using namespace std;using ll=long long;#define rep(i,n) for(int i=0;i<n;i++)#define rrep(i,n) for(int i=(n)-1;i>=0;i--)#define all(v) v.begin(),v.end()#define rall(v) v.rbegin(),v.rend()template<class T> bool chmax(T &a, T b){if (a < b){a = b;return true;} else return false;}template<class T> bool chmin(T &a, T b){if (a > b){a = b;return true;} else return false;}ll p,q;vector<ll>f(int n,vector<pair<ll,int>>&pr){if(n==pr.size())return{1};vector<ll>V=f(n+1,pr),re;auto[x,c]=pr[n];ll z=1;rep(j,c+1){for(auto y:V){re.push_back(y*z);}z*=x;}return re;}int main(){ios::sync_with_stdio(false);cin.tie(nullptr);cin>>p>>q;ll r=q;vector<pair<ll,ll>>ans;ll x=2;vector<pair<ll,int>>pr;while(x*x<=q){int cnt=0;while(q%x==0){cnt++;q/=x;}if(cnt>0)pr.push_back(make_pair(x,2*cnt));x++;}if(q>1)pr.push_back(make_pair(q,2));vector<ll>V=f(0,pr);for(auto y:V){if((y+r)%p==0&&(r*r/y+r)%p==0){ans.push_back(make_pair((y+r)/p,(r*r/y+r)/p));}}cout<<ans.size()<<"\n";sort(all(ans));for(auto[n,m]:ans){cout<<n<<" "<<m<<"\n";}}