OFFSET
0,3
LINKS
Robert Israel, Table of n, a(n) for n = 0..10000
MAPLE
f:= proc(n) local S, Sa, k, m, x, a;
m:= 1 + ilog10(n);
S:= map(t -> subs(t, x), [msolve(x^3=n, 10^m)]);
if S = [] then return 0 fi;
for k from 0 do
a:= ceil((n*10^k)^(1/3));
Sa:= map(s -> (s-a) mod 10^m, S);
x:= a + min(Sa);
if x^3 < (n+1)*10^k then return x fi;
od
end proc:
f(0):= 0:
map(f, [$0..100]); # Robert Israel, Nov 19 2019
CROSSREFS
KEYWORD
AUTHOR
Amarnath Murthy, Nov 20 2002
EXTENSIONS
More terms from Sascha Kurz, Jan 07 2003
STATUS
approved