This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A224320 #11 Aug 20 2014 12:32:22 %S A224320 2,5,7,11,17,41,47,71,107,167,179,197,449,859,1019,1061,1499,2089, %T A224320 16901,47717,56269,86269,11917049 %N A224320 Primes without "3" as a digit that remain prime when any single digit is replaced with "3". %C A224320 No more terms < 10^13. %H A224320 Carlos Rivera, <a href="http://www.primepuzzles.net/puzzles/puzz_591.htm">Puzzle 591</a> %t A224320 lst = {}; n = 3; 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, 86269}]; lst %t A224320 p3Q[n_]:=Module[{idn=IntegerDigits[n]},FreeQ[idn,3] && AllTrue[ FromDigits/@ Table[ReplacePart[idn,i->3],{i,IntegerLength[n]}],PrimeQ]]; Select[Prime[Range[10^6]],p3Q] (* The program uses the function AllTrue from Mathematica version 10 *) (* _Harvey P. Dale_, Aug 20 2014 *) %Y A224320 Cf. A224319, A224321-A224322. Subsequence of A038611. %K A224320 base,more,nonn %O A224320 1,1 %A A224320 _Arkadiusz Wesolowski_, Apr 03 2013