import math

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

divs = set()
for x in range(1, int(math.sqrt(q)) + 1):
  if q % x == 0:
    y = q // x
    divs.add(x)
    divs.add((q * q) // x)
    divs.add(y)
    divs.add((q * q) // y)

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}")