Pythonで基本ピタゴラス数の生成 - yattのブログ

Pythonで基本ピタゴラス数の生成 - yattのブログ

import itertools def gcd(x, y): while y: x,y = y,x%y return x def pytanum(): # 基本ピタゴラス数の...概要を表示 import itertools def gcd(x, y): while y: x,y = y,x%y return x def pytanum(): # 基本ピタゴラス数の生成 # [ピタゴラス数] # x^2 + y^2 = z^2 となる (x,y,z) の組 # [基本ピタゴラス数] # x,y,zが互いに素なピタゴラス数 # [基本ピタゴラス数の生成] # gcd a b = 1, a,bの一方が偶数かつ他方が奇数, a>b>0なら # 基本ピタゴラス数(a^2-b^2, 2ab, a^2+b^2) for a in itertools.count(2): a2 = a*a for b in xrange([1, 2][a%2], a, 2): if gcd(a, b) == 1: b2 = b*b yield a2-b2, 2*a*b, a2+b2 確認 >>> g = py