OFFSET
0,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 1..1000 from T. D. Noe)
FORMULA
a(n) >= n^(1/3), for all n > 0, with equality when n is a cube. - Derek Orr, Jul 26 2015
MAPLE
f:= proc(n) local d, m;
for d from 0 do
m:= ceil((10^d*n)^(1/3));
if m^3 < 10^d*(n+1) then return m fi
od
end proc:
map(f, [$1..100]); # Robert Israel, Jul 26 2015
PROG
(Python)
for n in range(1, 10**3):
for k in range(10**3):
if str(k**3).startswith(str(n)):
print(k, end=', ')
break
n += 1 # Derek Orr, Aug 03 2014
(PARI) a(n)=k=1; while(k, d=digits(k^3); D=digits(n); if(#D<=#d, c=1; for(i=1, #D, if(D[i]!=d[i], c=0; break)); if(c, return(k))); k++)
vector(100, n, a(n)) \\ Derek Orr, Jul 26 2015
CROSSREFS
KEYWORD
AUTHOR
EXTENSIONS
0 prepended by Seiichi Manyama, Jan 30 2017
STATUS
approved