cp's OEIS Frontend

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.

A050666 Primes p such that number of primes produced according to rules stipulated in Honaker's A048853 is 4.

This page as a plain text file.
%I A050666 #8 Oct 14 2019 23:55:22
%S A050666 97,389,659,743,953,2423,2971,3037,3083,3593,3727,4177,4363,4813,4889,
%T A050666 5261,5519,5717,5813,5939,6197,6449,7001,7057,7331,7963,8123,8179,
%U A050666 8353,8501,8707,9293,9587,10247,11261,13147,14813,15263,16547,17573,18077
%N A050666 Primes p such that number of primes produced according to rules stipulated in Honaker's A048853 is 4.
%e A050666 Altering a(1)=97 gives 4 primes: 17, 37, 47 and 67.
%t A050666 A048853[n_] := Module[{idn = IntegerDigits[n], id, np = 0}, Do[id = idn; If[id[[j]] != k, id[[j]] = k; If[id[[1]] != 0 && PrimeQ[FromDigits[id]], np = np + 1]], {j, 1, Length[idn]}, {k, 0, 9}]; np]; Select[Table[{p, A048853[p]}, {p, Prime[Range[3000]]}], #[[2]] == 4&][[All,1]] (* _Jean-François Alcover_, May 09 2017 *)
%Y A050666 Cf. A048853, A050673, A050655.
%K A050666 nonn,base
%O A050666 1,1
%A A050666 _Patrick De Geest_, Jul 15 1999