reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
Vincenzo Librandi, <a href="/A128106/b128106.txt">Table of n, a(n) for n = 1..10000</a>
approved
editing
editing
approved
(Sage)
def A128106_list(max):
R = []; s = 1; sq = 1
for n in (0..max//2):
if n == s:
sq += 1;
s = sq*sq;
for k in range(sq):
if is_square(n-k*k):
R.append(2*n)
break
R[0] = 1
return R
A128106_list(274) # Peter Luschny, Jun 20 2014
approved
editing
q=12; imax=2*q^2; lst=Select[Union[Flatten[Table[2*x^2+2*y^2, {x, 0, q}, {y, 0, x}]]], #<=imax&]; Join[{1}, Drop[lst, 1]] (* From _Vladimir Joseph Stephan Orlovsky, _, Apr 20 2011 *)
_T. D. Noe (noe(AT)sspectra.com), _, Feb 15 2007
proposed
approved
q=1712; imax=2*q^2; lst=Select[Union[Flatten[Table[2*x^2+2*y^2, {y, x, 0, q}, {x, y, 0, qx}]]], #<=imax&]; Join[{1}, Drop[lst, 1]] (* From Vladimir Joseph Stephan Orlovsky, Apr 20 2011 *)