def fact(n):
  #nの素因数分解
  arr=[]
  temp=n
  for i in range(2,int(-(-n**0.5//1))+1):
    if temp%i==0:
      cnt=0
      while temp%i==0:
        cnt+=1
        temp//=i
      arr.append([i, cnt])
  if temp!=1:
    arr.append([temp, 1])
  if arr==[]:
    arr.append([n, 1])
  return arr

def div(n):
  #n^2の約数列挙
  arr=[1]
  for i,j in fact(n):
    tmp=arr.copy()
    arr=[]
    for a in range(2*j+1):
      for k in tmp:
        arr.append(k*pow(i,a))
  return sorted(arr)


def solve(p,q):
  #(pa-q)*(pb-q)=q^2 を解く
  arr=[]
  check=set()
  for i in div(q):
    j=q**2//i
    if (i+q)%p==0 and (j+q)%p==0:
      if (i+q)//p not in check:
        print((i+q)//p,(j+q)//p)
        check.add((i+q)//p)

p,q=map(int,input().split())
solve(p,q)