OFFSET
0,1
COMMENTS
LINKS
Robert G. Wilson v, Table of n, a(n) for n = 0..10000
MATHEMATICA
f[s_List] := Block[{n = s[[-1]]}, sol = Solve[ x^2 + y^2 == z^2 && GCD[x, y, z] == 1 && x > 0 && y > 0 && z > 0 && (x == n || z == n), {x, y, z}, Integers]; Append[s, Min[ Complement[ Union[ Extract[ sol, Position[ sol, _Integer]]], s]]]]; lst = Nest[f, {3}, 25]
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Mar 16 2014
STATUS
approved