login
A091640
Number of primes less than 10^n which do not contain the digit 6.
10
4, 23, 136, 897, 6367, 46706, 355148, 2770239, 21984207, 176966593, 1440765209, 11838096715, 98014747908, 816769206831
OFFSET
1,1
FORMULA
Number of primes less than 10^n after removing any primes with at least one digit 6.
a(n) = A006880(n) - A091707(n).
EXAMPLE
a(2) = 23 because of the 25 primes less than 10^2, 2 have at least one digit 6; 25-2 = 23.
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], 6] == {}, c++ ]]; Print[c]; p--, {n, 1, 8}] (* Robert G. Wilson v, Feb 02 2004 *)
PROG
(Python)
from sympy import sieve # use slower primerange for larger terms
def a(n): return sum('6' not in str(p) for p in sieve.primerange(2, 10**n))
print([a(n) for n in range(1, 8)]) # Michael S. Branicky, Apr 23 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