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, 101): for b in range(a + 1, 101, 2): if gcd(a, b) == 1: is_ps, c = is_perfect_square((a * a) + (b * b)) if is_ps: print a, b, c