from math import sqrt rac = 1000 carres = [n*n - 1 for n in range(rac)] for a in range(1,rac): for b in range(a,rac): if a*b in carres: for c in range(b,rac): if a*c in carres and b*c in carres: d = a + b + c + 2*a*b*c + 2*sqrt(a*b+1)*sqrt(a*c+1)*sqrt(b*c+1) print(a,b,c,int(d))