#define _USE_MATH_DEFINES
#include <bits/stdc++.h>
using namespace std;
#define FOR(i,m,n) for(int i=(m);i<(n);++i)
#define REP(i,n) FOR(i,0,n)
#define ALL(v) (v).begin(),(v).end()
using ll = long long;
constexpr int INF = 0x3f3f3f3f;
constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL;
constexpr double EPS = 1e-8;
constexpr int MOD = 998244353;
// constexpr int MOD = 1000000007;
constexpr int DY4[]{1, 0, -1, 0}, DX4[]{0, -1, 0, 1};
constexpr int DY8[]{1, 1, 0, -1, -1, -1, 0, 1};
constexpr int DX8[]{0, -1, -1, -1, 0, 1, 1, 1};
template <typename T, typename U>
inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; }
template <typename T, typename U>
inline bool chmin(T& a, U b) { return a > b ? (a = b, true) : false; }
struct IOSetup {
  IOSetup() {
    std::cin.tie(nullptr);
    std::ios_base::sync_with_stdio(false);
    std::cout << fixed << setprecision(20);
  }
} iosetup;

template <typename T>
std::vector<T> divisor(const T n) {
  std::vector<T> res;
  T i = 1;
  for (; i * i < n; ++i) {
    if (n % i == 0) {
      res.emplace_back(i);
      res.emplace_back(n / i);
    }
  }
  if (i * i == n && n % i == 0) res.emplace_back(i);
  std::sort(res.begin(), res.end());
  return res;
}

int main() {
  int p, q; cin >> p >> q;
  const int g = gcd(p, q);
  p /= g; q /= g;
  const vector<int> d = divisor(q);
  set<ll> divisors;
  for (const int d1 : d) {
    for (const int d2 : d) {
      divisors.emplace(1LL * d1 * d2);
    }
  }
  set<pair<ll, ll>> ans;
  for (const ll a : divisors) {
    const ll b = 1LL * q * q / a;
    if ((a + q) % p == 0 && (b + q) % p == 0) ans.emplace((a + q) / p, (b + q) / p);
  }
  cout << ans.size() << '\n';
  for (const auto& [n, m] : ans) cout << n << ' ' << m << '\n';
  return 0;
}