OFFSET
1,1
COMMENTS
This sequence was originally described as the list of "congrua". But that name more properly refers to A256418.
Numbers of the form 4*(x^3*y-x*y^3) (where x,y are integers and x>=y). Squares of these numbers are of the form N^4-K^2 (where N belongs to A135786 and K to A135789 or A135790). Proof uses identity: (4(x^3y-xy^3))^2=(x^2+y^2)^4-(x^4 - 6x^2 y^2 + y^4)^2. - Artur Jasinski, Nov 29 2007, Nov 14 2008
MAPLE
N:= 10^5: # to get all terms <= N
select(`<=`, {seq(seq(4*(x^3*y-x*y^3), y=1..x-1), x=1..floor(sqrt(N/4+1)))}, N);
# If using Maple 11 or earlier, uncomment the following line
# sort(convert(%, list)); # Robert Israel, Apr 06 2015
MATHEMATICA
a = {}; Do[Do[w = 4x^3y - 4x y^3; If[w > 0 && w < 10000, AppendTo[a, w]], {x, y, 1000}], {y, 1, 1000}]; Union[a] (* Artur Jasinski, Nov 29 2007 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Henry Bottomley, Aug 02 2000
EXTENSIONS
Edited by N. J. A. Sloane, Apr 06 2015 at the suggestion of Robert Israel, Apr 03 2015
STATUS
approved