def div(n): S = set() i = 1 while i * i <= n: if n % i == 0: S.add(i) S.add(n // i) i += 1 S.add(1) return list(S) P, Q = map(int, input().split()) S = set() L = div(Q) for a in L: for b in L: S.add(a * b) S.add(-a * b) ans = [] Q2 = Q * Q for d1 in S: d2 = Q2 // d1 if (d1 + Q) <= 0 or (d2 + Q) <= 0: continue if (d1 + Q) % P or (d2 + Q) % P: continue ans.append(((d1 + Q)//P, (d2 + Q)//P)) print(len(ans)) for a, b in ans: print(a, b)