A224322 Primes without "9" as a digit that remain prime when any single digit is replaced with "9".

17, 107, 137, 347, 1013, 1433, 1613, 3767, 4337, 8623, 25633, 31114073

Offset

1,1

Comments

No more terms < 10^13.

Links

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

Carlos Rivera, Puzzle 591

Mathematica

lst = {}; n = 9; Do[If[PrimeQ[p], i = IntegerDigits[p]; If[FreeQ[i, n], t = 0; s = IntegerLength[p]; Do[If[PrimeQ@FromDigits@Insert[Drop[i, {d}], n, d], t++, Break[]], {d, s}]; If[t == s, AppendTo[lst, p]]]], {p, 25633}]; lst
pr9Q[n_]:=Module[{idn=IntegerDigits[n]}, FreeQ[idn, 9]&&AllTrue[Table[ FromDigits[ ReplacePart[ idn, i->9]], {i, Length[idn]}], PrimeQ]]; Select[ Prime[Range[2*10^6]], pr9Q] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Mar 30 2015 *)

Crossrefs

Cf. A224319-A224321. Subsequence of A038617.

Sequence in context: A268263 A343706 A282323 * A142321 A007682 A125327

Adjacent sequences: A224319 A224320 A224321 * A224323 A224324 A224325

Keyword

base,more,nonn

Author

Arkadiusz Wesolowski, Apr 03 2013

Status

approved