OFFSET
1,3
COMMENTS
For n>=0 and n = x*y*z*(x+y-z) with integers x>=y>=z then we can even find nonnegative solutions (x,y,z). However, if we restrict to z>=0 then there are no solutions (x,y,z) in case n<0.
The negative integers of the form x*y*z*(x+y-z) with integers x>=y>=z are the negatives of A213158 and in that case z<0.
Nonnegative integers of the form (a^2-c^2)*(b^2-c^2) with integers a>=b>=c.
Note that we must allow c<0 to represent n=12, 24, 40, ....
The negative integers of the form (a^2-c^2)*(b^2-c^2) with integers a>=b>=c are the negatives of A213158.
EXAMPLE
12 = (1)*(-2)*(-3)*((1)+(-2)-(-3)) with (x,y,z) = (1,-2,-3).
12 = 2*2*1*(2+2-1) with (x,y,z) = (2,2,1).
12 = ((0)^2-(-2)^2)*((-1)^2-(-2)^2) with (a,b,c) = (0,-1,-2).
12 = ((1)^2-(-2)^2)*((0)^2-(-2)^2) with (a,b,c) = (1,0,-2).
PROG
(PARI) {isa(n) = forvec( v = vector(3, i, [0, ceil(n^(1/2))]), if( n == v[1] * v[2] * v[3] * (v[3] + v[2] - v[1]), return(1)), 1)}
CROSSREFS
KEYWORD
nonn
AUTHOR
_Michael Somos_, May 18 2013
STATUS
approved