N = 10
for a in range(1, N):
for sa in [1, -1]:
for b in range(a + 1, N):
for sb in [1, -1]:
for c in range(b + 1, N):
for sc in [1, -1]:
for d in range(c + 1, N):
for sd in [1, -1]:
x = sa * a + sb * b + sc * c + sd * d
if x <= 0:
continue
a2 = a ** 2
b2 = b ** 2
c2 = c ** 2
d2 = d ** 2
x2 = sa * a2 + sb * b2 + sc * c2 + sd * d2
def pm(x):
return "+" if x > 0 else "-"
if x * x == x2:
print(
f"{pm(sa)}sqrt({a2}) {pm(sb)} sqrt({b2}) {pm(sc)} sqrt({c2}) {pm(sd)} sqrt({d2})"
f" = sqrt({pm(sa)}{a2} {pm(sb)} {b2} {pm(sc)} {c2} {pm(sd)} {d2})"
)