login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A091639
Number of primes less than 10^n which do not contain the digit 5.
10
3, 22, 136, 905, 6310, 46549, 354910, 2765749, 21955845, 176781643, 1439380189, 11827571824, 97933795005, 816144146010
OFFSET
1,1
COMMENTS
Number of primes less than 10^n after removing any primes with at least one digit 5.
FORMULA
a(n) = A006880(n) - A091706(n).
EXAMPLE
a(2) = 22 because of the 25 primes less than 10^2, 3 have at least one digit 5; 25-3 = 22.
MATHEMATICA
NextPrim[n_] := Block[{k = n + 1}, While[ !PrimeQ[k], k++ ]; k]; c = 0; p = 1; Do[ While[ p = NextPrim[p]; p < 10^n, If[ Position[ IntegerDigits[p], 5] == {}, c++ ]]; Print[c]; p--, {n, 1, 8}] (* Robert G. Wilson v, Feb 02 2004 *)
PROG
(Python)
from sympy import primerange
def a(n): return sum('5' not in str(p) for p in primerange(2, 10**n))
print([a(n) for n in range(1, 7)]) # Michael S. Branicky, Mar 22 2021
KEYWORD
more,nonn,base
AUTHOR
Enoch Haga, Jan 30 2004
EXTENSIONS
Edited and extended by Robert G. Wilson v, Feb 02 2004
a(9)-a(12) from Donovan Johnson, Feb 14 2008
a(13) from Robert Price, Nov 08 2013
a(14) from Giovanni Resta, Mar 20 2017
STATUS
approved