OFFSET
1,2
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
FORMULA
Sum_{n>=1} 1/a(n) = (7*10)/((7-1)*(10-1)) = 35/27. - Amiram Eldar, Sep 25 2020
a(n) ~ exp(sqrt(2*log(7)*log(10)*n)) / sqrt(70). - Vaclav Kotesovec, Sep 25 2020
MATHEMATICA
n = 10^6; Flatten[Table[7^i*10^j, {i, 0, Log[7, n]}, {j, 0, Log10[n/7^i]}]] // Sort (* Amiram Eldar, Sep 25 2020 *)
PROG
(Haskell)
import Data.Set (singleton, deleteFindMin, insert)
a025632 n = a025632_list !! (n-1)
a025632_list = f $ singleton (1, 0, 0) where
f s = y : f (insert (7 * y, i + 1, j) $ insert (10 * y, i, j + 1) s')
where ((y, i, j), s') = deleteFindMin s
-- Reinhard Zumkeller, May 15 2015
(PARI) list(lim)=my(v=List(), N); for(n=0, logint(lim\=1, 10), N=10^n; while(N<=lim, listput(v, N); N*=7)); Set(v) \\ Charles R Greathouse IV, Jan 10 2018
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
STATUS
approved