login
Imprimitively represented by x^2+2y^2.
5

%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