OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 10^4: # for terms <= N
S:= {}:
for x from 1 to floor(N/2) do
S:= S union {seq(3*x*y^2 - x^3, y = ceil(x/sqrt(3)) .. floor(sqrt((x^2+N/x)/3)))}
od:
sort(convert(S, list)); # Robert Israel, Aug 14 2023
MATHEMATICA
a = {}; Do[Do[w = 3x y^2 - x^3; If[(w > 0) && w < 3000, AppendTo[a, w]], {x, 1, 1000}], {y, 1, 1000}]; Union[a]
With[{nn=300}, Select[Flatten[Table[3*x*y^2-x^3, {x, nn}, {y, nn}]], 0<#<600&]//Union] (* Harvey P. Dale, Jul 03 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Nov 29 2007
STATUS
approved