day1/exercise/pytriads.py
author Puneeth Chaganti <punchagan@fossee.in>
Wed, 28 Jul 2010 19:50:00 +0530
branchscipy2010
changeset 438 8af5dfa5432b
parent 428 ee689f7c3122
permissions -rw-r--r--
Added FFT stuff to day1/cheatsheet6.

def is_perfect_square(n):
    i = 1
    while i * i < n:
        i += 1
    return i * i == n, i

def 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