login
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