p142

p142.py

def is_sq(x):
    sqrt = (x ** (1 / 2))
    if round(sqrt) ** 2 == x:
        return True
    return False

def compute():
    for b in range(1, 1000):
        for a in range(b + 2, 1000, 2):
            x = (a*a + b*b)//2
            y = a*a - x
            if x < y:
                break
            else:
                for c in range(int(math.sqrt(x)), a):
                    z = c*c - x
                    if z > y:
                        break
                    else:
                        if all([is_sq(x - z), is_sq(y + z), is_sq(y - z)]):
                            return x + y + z