OFFSET
1,2
COMMENTS
A001221(a(n)) <= 3 since 10 cannot divide a(n).
It seems that this sequence is finite and contains 12615 terms. - Daniel Mondot, May 03 2022 and Jianing Song, Jan 28 2023
LINKS
Daniel Mondot, Table of n, a(n) for n = 1..12615 (terms 1..10000 from Charles R Greathouse IV)
EXAMPLE
a(12615) = 2^25 * 3^227 * 7^28.
PROG
(Haskell)
import Data.Set (singleton, deleteFindMin, fromList, union)
a238985 n = a238985_list !! (n-1)
a238985_list = filter ((== 1) . a168046) $ f $ singleton 1 where
f s = x : f (s' `union` fromList
(filter ((> 0) . (`mod` 10)) $ map (* x) [2, 3, 5, 7]))
where (x, s') = deleteFindMin s
(PARI) zf(n)=vecmin(digits(n))
list(lim)=my(v=List(), t, t1); for(e=0, log(lim+1)\log(7), t1=7^e; for(f=0, log(lim\t1+1)\log(3), t=t1*3^f; while(t<=lim, if(zf(t), listput(v, t)); t<<=1)); for(f=0, log(lim\t1+1)\log(5), t=t1*5^f; while(t<=lim, if(zf(t), listput(v, t)); t*=3))); Set(v)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Charles R Greathouse IV and Reinhard Zumkeller, Mar 07 2014
EXTENSIONS
Keyword:fini and keyword:full removed by Jianing Song, Jan 28 2023 as finiteness is only conjectured.
STATUS
approved