login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A083121
Bases of the squares arising in A083119.
2
1, 3, 4, 5, 6, 7, 9, 10, 12, 13, 15, 16, 17, 18, 19, 20, 22, 23, 24, 25, 27, 28, 30, 31, 32, 33, 35, 36, 38, 39, 41, 42, 43, 44, 45, 46, 48, 50, 51, 52, 53, 55, 57, 58, 59, 60, 61, 62, 63, 65, 67, 68, 70, 71, 72, 73, 74, 75, 76, 77
OFFSET
1,2
MAPLE
A083121 := proc(nmax) local a119, anew, asum, n, a; a119 := [1] ; a := [1] ; while nops(a119) < nmax do n := nops(a119)+1 ; if n mod 2 = 0 then anew := 1 ; while anew in a119 do anew := anew+1 ; od ; a119 := [op(a119), anew] ; else asum := add(op(i, a119), i=ceil(n/2)..n-1) ; anew := 1 ; while anew in a119 or not issqr(asum+anew) do anew := anew+1 ; od ; a119 := [op(a119), anew] ; a := [op(a), sqrt(asum+anew)] ; fi ; od; a ; end: A083121(120) ; # R. J. Mathar, Jul 17 2007
CROSSREFS
Cf. A083119.
Sequence in context: A253200 A227938 A039096 * A298006 A026438 A026442
KEYWORD
nonn
AUTHOR
Amarnath Murthy and Meenakshi Srikanth (menakan_s(AT)yahoo.com), Apr 23 2003
EXTENSIONS
More terms from R. J. Mathar, Jul 17 2007
STATUS
approved