p9.py

from math import sqrt
for a in range(1000):
    for b in range(a, 1000, 1):
        c = sqrt(a**2 + b**2)
        if a + b + c == 1000:
            print(a * b * c)