OFFSET
0,3
LINKS
Robert Israel, Table of n, a(n) for n = 0..10000
MAPLE
f:= proc(n) local m, d, x, R, Rn, n1, n2, n0;
m:= 1+ilog10(n);
R:= map(t -> rhs(op(t)), {msolve(x^3=n, 10^m)});
if R = {} then return 0 fi;
for d from 0 do
n1:= ceil(surd(n*10^d, 3));
n2:= ceil(surd((n+1)*10^d, 3))-1;
Rn:= (R -~ n1) mod 10^m;
n0:= min(Rn) + n1;
if n0 <= n2 then return n0^3 fi
od
end proc:
f(0):= 0:
map(f, [$0..50]); # Robert Israel, Apr 16 2026
CROSSREFS
KEYWORD
AUTHOR
Amarnath Murthy, Nov 20 2002
EXTENSIONS
More terms from Sascha Kurz, Jan 07 2003
Corrected and extended by Franklin T. Adams-Watters, Jul 19 2006
STATUS
approved
