def is_perfect_square(n): i = 1 while i * i < n: i += 1 return i * i == n, idef gcd(a, b): if a % b == 0: return b else: return gcd(b, a%b)for a in range(3, 501): for b in range( a+1, 501, 2): if gcd( a, b ) == 1: is_ps, c = is_perfect_square((a * a) + (b * b)) if is_ps: print a, b, c