Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #37 Oct 26 2020 09:36:21
%S 0,1,4,9,16,25,36,64,81,100,121,144,225,256,289,324,400,484,529,576,
%T 625,729,841,900,1024,1089,1156,1296,1600,1681,1936,2025,2116,2209,
%U 2304,2500,2601,2809,2916,3025,3364,3481,3600,4096,4356,4624,4761,5041,5184,5625,6400,6561,6724,6889,7225,7569
%N Squares that are not of the form x^2 + x*y + y^2, where x and y are positive integers.
%C Or Loeschian numbers (A003136) that are not in A024614.
%C Squares that are not in this sequence are 49, 169, 196, 361, 441, 676, ...
%C This is the list of squares not of the form A050931(k)^2. A number n is in this sequence iff n = m^2 with m having no prime factor == 1 (mod 6). - _M. F. Hasler_, Mar 04 2018
%F a(n) = A230780(n-1)^2 for n > 1.
%e Loeschian number 25 = 5^2 is a term because 25 = x^2 + x*y + y^2 has no solution for positive integers x, y.
%p isA024614:= proc(n) local x,y;
%p for x from 1 to floor(sqrt(n-1)) do
%p if issqr(4*n-3*x^2) then return true fi
%p od:
%p false
%p end proc:
%p isA024614(0):= false:
%p remove(isA024614, [seq(i^2,i=0..200)]); # _Robert Israel_, Mar 02 2018
%t sol[s_] := Solve[0 < x <= y && s == x^2 + x y + y^2, {x, y}, Integers];
%t Select[Range[0, 100]^2, sol[#] == {}&] (* _Jean-François Alcover_, Oct 26 2020 *)
%o (PARI) is(n,m)=issquare(n,m)&&!setsearch(Set(factor(m)[,1]%6),1) \\ second part is equivalent to is_A230780(m), this is sufficient to test (e.g., to produce a list) if we know that n = m^2. - _M. F. Hasler_, Mar 04 2018
%Y Cf. A000548, A003136, A024614, A050931, A230780.
%K nonn,easy
%O 1,3
%A _Altug Alkan_, Mar 02 2018