import math

p, q = map(int, input().split())
g = math.gcd(p, q)
p //= g
q //= g

divs = [q]
for x in range(1, q):
  if (q * q) % x == 0:
    divs.append(x)
    if x != (q * q) // x:
      divs.append((q * q) // x)

ans = []
for x in divs:
  y = (q * q) // x
  if (x + q) % p != 0 or (y + q) % p != 0:
    continue
  m = (x + q) // p
  n = (y + q) // p
  if m > 0 and n > 0:
    ans.append((m, n))

ans.sort()
print(len(ans))
for m, n in ans:
  print(f"{m} {n}")