OFFSET
1,1
COMMENTS
Subsequence of primes of A052419. - Michel Marcus, Feb 22 2015
Maynard proves that this sequence is infinite and in particular contains the expected number of elements up to x, on the order of x^(log 9/log 10)/log x. - Charles R Greathouse IV, Apr 08 2016
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..20000
Marianne Freiberger, Primes without 7s, +plus magazine, August 2016.
M. F. Hasler, Numbers avoiding certain digits, OEIS wiki, Jan 12 2020.
James Maynard, Primes with restricted digits, arXiv:1604.01041 [math.NT], 2016.
James Maynard and Brady Haran, Primes without a 7, Numberphile video (2019).
FORMULA
a(n) ~ n^(log 10/log 9) log n. - Charles R Greathouse IV, Aug 03 2023
MATHEMATICA
Select[Prime[Range[70]], DigitCount[#, 10, 7] == 0 &] (* Vincenzo Librandi, Aug 08 2011 *)
PROG
(Magma) [ p: p in PrimesUpTo(500) | not 7 in Intseq(p) ]; // Bruno Berselli, Aug 08 2011
(PARI) lista(nn)=forprime(p=2, nn, if (!vecsearch(vecsort(digits(p), , 8), 7), print1(p, ", ")); ); \\ Michel Marcus, Feb 22 2015
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
Vasiliy Danilov (danilovv(AT)usa.net), Jul 15 1998
EXTENSIONS
Offset corrected by Arkadiusz Wesolowski, Aug 07 2011
STATUS
approved