OFFSET
1,1
COMMENTS
Squares of these numbers are of the form N^4 - M^2 (where N belongs to A135786 and M to A057102). Proof is based on the identity (x^4 - 6x^2 * y^2 + y^4)^2 = (x^2 + y^2)^4 - (4(x^3y - xy^3))^2.
Since x^4 - 6x^2 * y^2 + y^4 = d*d' where d = x^2 - y^2 + 2xy and d' = x^2 - y^2 - 2xy, and d - d' = 4xy, the computational technique is to consider the divisors d|n, d'=n/d, to check that the difference is a multiple of 4, and to check x in the range 1..d/3. - R. J. Mathar, Sep 18 2009
Refers to A057102, which had an incorrect description and has been replaced by A256418. As a result the present sequence should be re-checked. - N. J. A. Sloane, Apr 06 2015
MAPLE
isA135789 := proc(n) for d in numtheory[divisors](n) do dprime := n/d ; if abs(d-dprime) mod 4 = 0 then for x from 1 to d/3 do y := (d-dprime)/4/x ; if type(y, 'integer') and y< x and y> 0 then if n = (x^2-y^2+2*x*y)*(x^2-y^2-2*x*y) then RETURN(true); fi; fi; od: fi: od: RETURN(false) ; end: for n from 1 do if isA135789(n) then printf("%d, \n", n) ; fi; od: # R. J. Mathar, Sep 18 2009
MATHEMATICA
a = {}; Do[Do[w = x^4 - 6x^2 y^2 + y^4; If[w > 0&&w<10000, AppendTo[a, w]], {x, y, 2000}], {y, 1, 2000}]; Union[a]
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Nov 29 2007, Nov 14 2008
EXTENSIONS
More terms from R. J. Mathar, Sep 18 2009
STATUS
approved