cp's OEIS Frontend

A034030 Imprimitively represented by x^2+2y^2.

%I A034030 #15 Jul 08 2025 20:03:08
%S A034030 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 A034030 96,98,99,100,108,121,128,132,136,144,147,150,152,153,162,164,169,171,
%U A034030 172,176,192,196,198,200,204,216,225,228
%N A034030 Imprimitively represented by x^2+2y^2.
%p A034030 # Maple code for A002479, A057127, A034030-A034034 from _N. J. A. Sloane_, Apr 30 2015
%p A034030 lis:={}; lisP:={}; lisI:={};
%p A034030 M:=50; M2:=M^2;
%p A034030 for x from 0 to M do
%p A034030    x2:=x^2;
%p A034030 for y from 0 to M do
%p A034030    N:=x2+2*y^2;
%p A034030 if N <= M2 then
%p A034030    if gcd(x,y) = 1 then lisP:={op(lisP),N}; else lisI:={op(lisI),N} fi;
%p A034030    lis:={op(lis),N};
%p A034030 fi;
%p A034030 od: od:
%p A034030 lprint("lis");
%p A034030 Lis:=sort(convert(lis,list));
%p A034030 lprint("lisP");
%p A034030 LisP:=sort(convert(lisP,list));
%p A034030 lprint("lisI");
%p A034030 LisI:=sort(convert(lisI,list));
%p A034030 lprint("lisPnotI");
%p A034030 LisPnotI:=sort(convert(lisP minus lisI, list));
%p A034030 lprint("lisInotP");
%p A034030 LisInotP:=sort(convert(lisI minus lisP,list));
%p A034030 lprint("lisIandP");
%p A034030 LisIandP:=sort(convert(lisI intersect lisP,list));
%p A034030 lprint("liseither");
%p A034030 Liseither:=sort(convert(lis minus (lisI intersect lisP),list));
%Y A034030 Cf. A002479, A097700, A057127, A034031-A034034.
%K A034030 nonn
%O A034030 1,2
%A A034030 _N. J. A. Sloane_
%E A034030 Corrected by _N. J. A. Sloane_, Apr 30 2015