# difficult.py

import math

def isPerfect(x):
  sx = math.sqrt(x)
  return sx == int(sx)

print("  a    b   q")
for a in range(1, 151):
  for b in range(1, 1001):
    q = (a*a + b*b)/(a*b+1)
    if (q == int(q)) and isPerfect(q):
      print(f"{a:3d} {b:4d} {int(q):3d}")