cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A270248 Even Löschian numbers.

This page as a plain text file.
%I A270248 #13 Jul 05 2017 13:32:32
%S A270248 0,4,12,16,28,36,48,52,64,76,84,100,108,112,124,144,148,156,172,192,
%T A270248 196,208,228,244,252,256,268,292,300,304,316,324,336,364,372,388,400,
%U A270248 412,432,436,444,448,468,484,496,508,516,532,556,576,588,592,604,624,628,652,676
%N A270248 Even Löschian numbers.
%C A270248 Even numbers of the form x^2 - xy + y^2.
%H A270248 Charles R Greathouse IV, <a href="/A270248/b270248.txt">Table of n, a(n) for n = 1..10000</a>
%F A270248 a(n) = 2 * A270050(n) = 4 * A003136(n).
%e A270248 Even number 12 is a term because 12 = 2^2 + 2*2 + 2^2.
%t A270248 Select[Range[0, 680, 2], Resolve@ Exists[{x, y}, Reduce[# == (x^2 - x y + y^2), {x, y}, Integers]] &] (* _Michael De Vlieger_, Mar 15 2016 *)
%o A270248 (PARI) x='x+O('x^800); p=eta(x)^3/eta(x^3); for(n=0, 799, if(polcoeff(p, n) != 0 && n % 2 == 0, print1(n, ", ")));
%o A270248 (PARI) list(lim)=my(v=List(), y, t); forstep(x=0, sqrtint(lim\3), 2, my(y=x, t); while((t=x^2+x*y+y^2)<=lim, listput(v, t); y+=2)); Set(v) \\ _Charles R Greathouse IV_, Jul 05 2017
%Y A270248 Cf. Loeschian numbers: A003136 (all), A266836 (2*k+1), A202822 (3*k+1), A260682 (6*k+1).
%K A270248 nonn
%O A270248 1,2
%A A270248 _Altug Alkan_, Mar 14 2016