A083121 - OEIS

%I #8 Aug 07 2015 03:03:53

%S 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,

%T 33,35,36,38,39,41,42,43,44,45,46,48,50,51,52,53,55,57,58,59,60,61,62,

%U 63,65,67,68,70,71,72,73,74,75,76,77

%N Bases of the squares arising in A083119.

%p 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

%Y Cf. A083119.

%K nonn

%O 1,2

%A _Amarnath Murthy_ and Meenakshi Srikanth (menakan_s(AT)yahoo.com), Apr 23 2003

%E More terms from _R. J. Mathar_, Jul 17 2007