OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 10^5:
A:= {}:
for v from 1 to floor((N/2)^(1/3)) do
for u from v+1 do
if igcd(u, v) = 1 then
t:= u*v*(u^2-v^2);
if t > N then break fi;
A:= A union {t};
fi
od
od:
A;
# if using Maple 11 or earlier, uncomment the next line
# sort(convert(A, list)); # Robert Israel, Apr 06 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Max Alekseyev, Nov 12 2008
STATUS
approved