login
A018928
Define {b(n)} by b(1)=3, b(n) (n >= 2) is the smallest number such that b(1)^2 + ... + b(n)^2 = m^2 for some m and all b(i) are distinct. Sequence gives values of m.
6
3, 5, 13, 85, 157, 12325, 12461, 106285, 276341, 339709, 10363909, 17238541, 1936511509, 51335823965, 133473142309, 872709007405, 1574530008629, 667511933218429, 698925273030725, 707670964169285, 1839944506840141
OFFSET
1,1
COMMENTS
Also: Begin with the least length of a Pythagorean Triangle (PT), a(1)=3. Then a(n) is the least hypotenuse of a PT which has a(n-1) as one of its legs. - Robert G. Wilson v, Mar 17 2014
LINKS
Robert G. Wilson v, Table of n, a(n) for n = 1..165 (first 100 terms from Lei Zhou)
MATHEMATICA
NextA018928[n_] := Block[{a = n^2, b, l, i, c, d, f}, b = Divisors[a]; l = Length[b]; i = l; While[i--; c = b[[i]]; d = a/c - (c - 1); (d <= 1) || EvenQ[d]]; f = (a/c + (c - 1) + 1)/2]; Table[If[i == 1, a = 3, a = NextA018928[a]]; a, {i, 1, 21}](* Lei Zhou, Feb 20 2014 *)
f[s_List] := Block[{x = s[[-1]]}, Append[s, Transpose[ Solve[ x^2 + y^2 == z^2 && y > 0 && z > 0, {y, z}, Integers]][[-1, 1, 2]]]]; lst = Nest[f, lst, 15] (* Robert G. Wilson v, Mar 17 2014 *)
CROSSREFS
Sequence in context: A159293 A051901 A268021 * A239381 A180313 A053630
KEYWORD
nonn
AUTHOR
Charles Reed (charles.reed(AT)bbs.ewgateway.org)
EXTENSIONS
More terms from David W. Wilson
STATUS
approved