|
%I #39 Mar 28 2019 07:06:35
%S 3,6,7,11,12,21,28,33,35,38,42,44,47,62,70,71,77,83,84,91,102,107,110,
%T 119,132,133,140,147,161,168,182,188,189,203,217,231,238,266,267,278,
%U 280,284,287,294,308,315,318,329,332,335,336,343,357,364,378,385
%N Numbers not of the form x^2+y^2+14*z^2.
%H Robert Israel, <a href="/A306718/b306718.txt">Table of n, a(n) for n = 1..10000</a>
%H Irving Kaplansky, <a href="https://doi.org/10.1090/S0025-5718-1995-1265017-2">The first nontrivial genus of positive definite ternary forms</a>, Mathematics of Computation, Vol. 64 (1995): 341-345.
%p N:= 1000: # for terms <= N
%p V:= Array(0..N):
%p for x from 0 to floor(sqrt(N/2)) do
%p for y from x while x^2 + y^2 <= N do
%p for z from 0 do
%p v:= x^2 + y^2 + 14*z^2;
%p if v > N then break fi;
%p V[v]:= 1;
%p od od od:
%p select(t -> V[t]=0, [$1..N]); # _Robert Israel_, Mar 27 2019
%K nonn
%O 1,1
%A _N. J. A. Sloane_, Mar 26 2019
|
