%I #13 Apr 30 2015 13:59:20
%S 0,4,8,9,12,16,18,24,25,27,32,36,44,48,49,50,54,64,68,72,75,76,81,88,
%T 96,98,99,100,108,121,128,132,136,144,147,150,152,153,162,164,169,171,
%U 172,176,192,196,198,200,204,216,225,228
%N Imprimitively represented by x^2+2y^2.
%p # Maple code for A002479, A057127, A034030-A034034 from _N. J. A. Sloane_, Apr 30 2015
%p lis:={}; lisP:={}; lisI:={};
%p M:=50; M2:=M^2;
%p for x from 0 to M do
%p x2:=x^2;
%p for y from 0 to M do
%p N:=x2+2*y^2;
%p if N <= M2 then
%p if gcd(x,y) = 1 then lisP:={op(lisP),N}; else lisI:={op(lisI),N} fi;
%p lis:={op(lis),N};
%p fi;
%p od: od:
%p lprint("lis");
%p Lis:=sort(convert(lis,list));
%p lprint("lisP");
%p LisP:=sort(convert(lisP,list));
%p lprint("lisI");
%p LisI:=sort(convert(lisI,list));
%p lprint("lisPnotI");
%p LisPnotI:=sort(convert(lisP minus lisI, list));
%p lprint("lisInotP");
%p LisInotP:=sort(convert(lisI minus lisP,list));
%p lprint("lisIandP");
%p LisIandP:=sort(convert(lisI intersect lisP,list));
%p lprint("liseither");
%p Liseither:=sort(convert(lis minus (lisI intersect lisP),list));
%Y Cf. A002479, A097700, A057127, A034031-A034034.
%K nonn
%O 1,2
%A _N. J. A. Sloane_.
%E Corrected by _N. J. A. Sloane_, Apr 30 2015