#include<bits/stdc++.h>
using namespace std;
long long P,Q;
vector<array<long long,2>>ans;
 void dfs(long long t,vector<array<long long,2>>& v,int c){
  if(c==v.size()){
    if((t+Q)%P==0&&(Q*Q/t+Q)%P==0){
      ans.push_back({(t+Q)/P,(Q*Q/t+Q)/P});
      ans.push_back({(Q*Q/t+Q)/P,(t+Q)/P});
    }
    if((-t+Q)%P==0&&(-Q*Q/t+Q)%P==0&&-t+Q!=0&&-Q*Q/t+Q!=0&&-t+Q>0&&-Q*Q/t+Q>0){
      ans.push_back({(-t+Q)/P,(-Q*Q/t+Q)/P});
      ans.push_back({(-Q*Q/t+Q)/P,(-t+Q)/P});
    }
    return;
  }
  dfs(t,v,c+1);
  for(int i=0;i<v[c][1];i++){
    t*=v[c][0];
    dfs(t,v,c+1);
  }
  for(int i=0;i<v[c][1];i++){
    t/=v[c][0];
  }
}

int main(){
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  cin>>P>>Q;
  long long g=gcd(P,Q);
  P/=g,Q/=g;
  vector<array<long long,2>>d;
  long long Q2=Q;
  for(long long i=2;i*i<=Q;i++){
    if(Q%i)continue;
    long long cnt=0;
    while(Q2%i==0){
      Q2/=i;
      ++cnt;
    }
    d.push_back({i,cnt});
  }
  if(Q2!=1)d.push_back({Q2,1});
  for(auto &i:d)i[1]*=2;
  dfs(1,d,0);
  sort(ans.begin(),ans.end());
  ans.erase(unique(ans.begin(),ans.end()),ans.end());
  cout<<ans.size()<<'\n';
  for(auto i:ans)cout<<i[0]<<' '<<i[1]<<'\n';
}