OFFSET
1,3
COMMENTS
For k > 0, write k = s * 10^t, 1 <= s < 10, then k is a term if and only if s is in [1, 2^(1/3)) U (20^(1/3), 3) U (30^(1/3), 40^(1/3)) U (30^(2/3), 10). - Jianing Song, Dec 26 2022
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
EXAMPLE
a(13)=99 and 970299=99^3 start both with 9.
MATHEMATICA
Select[Range[0, 300], First[IntegerDigits[#]]== First[IntegerDigits[#^3]]&] (* Harvey P. Dale, Mar 14 2011 *)
PROG
(PARI) isok(n) = (n == 0) || (digits(n)[1] == digits(n^3)[1]); \\ Michel Marcus, Mar 18 2015
(Haskell)
a144582 n = a144582_list !! (n-1)
a144582_list = [x | x <- [0..], a000030 x == a000030 (x ^ 3)]
-- Reinhard Zumkeller, Apr 01 2015
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Reinhard Zumkeller, Aug 17 2008
STATUS
approved