OFFSET
0,2
COMMENTS
Clearly each term is even as ax^3 + bx^2 + cx + d = 0 and -ax^3 - bx^2 - cx - d = 0 have the same roots.
The variable D in the PARI program below is the discriminant of the reduced form y^3 + py + q = 0.
REFERENCES
Jan Gullberg, Mathematics, From the Birth of Numbers, W. W. Norton & Co., NY, pages 318-9.
PROG
(PARI) {for(n=0, 32, c=0; forvec(xx=[[ -n, n], [ -n, n], [ -n, n], [ -n, n]],
if(xx[1]==0, next, z=Pol(xx); x=y-xx[2]/(3*xx[1]);
zz=eval(z); if(polcoeff(zz, 3)<>1, zz=zz/polcoeff(zz, 3));
p=polcoeff(zz, 1); q=polcoeff(zz, 0); D=(q/2)^2+(p/3)^3;
if(D>0, c++))); print1(c, ", "))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Rick L. Shepherd, Jan 21 2009
STATUS
approved