login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A018928 Define {b(n)} by b(1)=3, b(n) (n >= 2) is 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 (list; graph; refs; listen; history; text; internal format)
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

Cf. A018929, A018930.

Sequence in context: A159293 A051901 A268021 * A239381 A180313 A053630

Adjacent sequences:  A018925 A018926 A018927 * A018929 A018930 A018931

KEYWORD

nonn

AUTHOR

Charles Reed (charles.reed(AT)bbs.ewgateway.org)

EXTENSIONS

More terms from David W. Wilson

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 5 00:43 EDT 2020. Contains 333238 sequences. (Running on oeis4.)