OFFSET
1,3
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
M. F. Hasler, Numbers avoiding certain digits OEIS wiki, Jan 12 2020.
FORMULA
a(n) = replace digits d > 5 by d + 1 in base-9 representation of n - 1. - Reinhard Zumkeller, Oct 07 2014
Sum_{k>1} 1/a(k) = A082835 = 22.205598... (Kempner series). - Bernard Schott, Jan 12 2020, edited by M. F. Hasler, Jan 13 2020
MAPLE
a:= proc(n) local l, m; l, m:= 0, n-1;
while m>0 do l:= (d->
`if`(d<6, d, d+1))(irem(m, 9, 'm')), l
od; parse(cat(l))/10
end:
seq(a(n), n=1..100); # Alois P. Heinz, Aug 01 2016
MATHEMATICA
Select[Range[0, 100], DigitCount[#, 10, 6]==0&] (* Harvey P. Dale, Jun 20 2013 *)
PROG
(Magma) [ n: n in [0..89] | not 6 in Intseq(n) ]; // Bruno Berselli, May 28 2011
(sh) seq 0 1000 | grep -v 6; # Joerg Arndt, May 29 2011
(Haskell)
a052414 = f . subtract 1 where
f 0 = 0
f v = 10 * f w + if r > 5 then r + 1 else r where (w, r) = divMod v 9
-- Reinhard Zumkeller, Oct 07 2014
(PARI) lista(nn)=for (n=0, nn, if (!vecsearch(vecsort(digits(n), , 8), 6), print1(n, ", ")); ); \\ Michel Marcus, Feb 22 2015
(PARI) /* See OEIS wiki page (cf. LINKS) for more programs */
apply( {A052414(n)=fromdigits(apply(d->d+(d>5), digits(n-1, 9)))}, [1..99]) \\ a(n)
next_A052414(n, d=digits(n+=1))={for(i=1, #d, d[i]==6&&return((1+n\d=10^(#d-i))*d)); n} \\ least a(k) > n, used in A038614. - M. F. Hasler, Jan 11 2020
CROSSREFS
KEYWORD
base,easy,nonn
AUTHOR
Henry Bottomley, Mar 13 2000
EXTENSIONS
Offset changed by Reinhard Zumkeller, Oct 07 2014
STATUS
approved