A091645 Number of primes less than 10^n having at least one digit 1.

0, 8, 67, 559, 4917, 44073, 402030, 3709989, 34534791, 323587790, 3046685950, 28798331502, 273078628285, 2596399340798

Offset

1,2

Links

Table of n, a(n) for n=1..14.

Example

a(2) = 8 because of the 25 primes less than 10^2, 8 have at least one 1 digit.

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], 1] != {}, c++ ]]; Print[c]; p--, {n, 1, 8}] (* Robert G. Wilson v, Feb 02 2004 *)
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], 2] != {}, c++ ]]; Print[c]; p--, {n, 1, 8}] (* Robert G. Wilson v, Feb 02 2004 *)

Crossrefs

a(n) + A091635(n) = A006880(n).

Cf. A091645, A091646, A091647, A091705, A091706, A091707, A091708, A091709, A091710.

Sequence in context: A037594 A037685 A228419 * A015574 A370028 A152055

Adjacent sequences: A091642 A091643 A091644 * A091646 A091647 A091648

Keyword

more,nonn,base

Author

Enoch Haga, Jan 30 2004

Extensions

Edited, corrected and extended by Robert G. Wilson v, Feb 02 2004

3 more terms from Ryan Propper, Aug 20 2005

a(13) from Robert Price, Nov 11 2013

a(14) from Giovanni Resta, Jul 21 2015

Status

approved