import math for a in range(1, 333): modulos = (1000**2 - 2000 * a) % (2000 - 2 * a) # b = (1000^2 - 2000 * a) / (2000 - 2 * a) if (modulos == 0): b = (1000**2 - 2000 * a) / (2000 - 2 * a) c = math.sqrt(a**2 + b**2) print(a, b, c) print(a * b * c) break # for b in range(333): # c = math.sqrt(a**2 + b**2) # if (a + b + c > 1000): # break # if (a + b + c == 1000): # print('FOUND!') # print(a, b, c)