OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 1000: # to get all terms <= N
X:= floor(sqrt(N/3)):
V:= Vector(N):
for x from 2 to X do
if x^3 > N then
y0:= iroot(x^3-N, 3);
if x^3 - y0^3 > N then y0:= y0+1 fi;
else y0:= 1 fi;
for y from y0 to x-1 do
V[x^3 - y^3] := V[x^3 - y^3]+1
od
od: select(t -> V[t] = 1, [$1..N]); # Robert Israel, Dec 11 2015
MATHEMATICA
r = 992; p = 3; Sort@Drop[Flatten@Select[Tally@Reap[Do[n = i^p - j^p; If[n <= r, Sow[n]], {i, Ceiling[(r/p)^(1/(p - 1))]}, {j, i}]][[2, 1]], #[[2]] == 1 &], {2, -1, 2}] (* Arkadiusz Wesolowski, Dec 10 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Glen Burch (gburch(AT)erols.com)
EXTENSIONS
Corrected and extended by Don Reble, Nov 19 2006
STATUS
approved