OFFSET
1,1
COMMENTS
This is a subsequence of A019546.
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000
MATHEMATICA
Select[Table[Prime[n], {n, 7200}], ContainsExactly[IntegerDigits[#], {2, 3, 5, 7}]&] (* James C. McMahon, Mar 05 2024 *)
PROG
(Python)
from sympy import isprime
from itertools import count, islice, product
def agen(): yield from (t for d in count(4) for b in product("2357", repeat=d-1) for e in "37" if len(set(b+(e, )))==4 and isprime(t:=int("".join(b)+e)))
print(list(islice(agen(), 40))) # Michael S. Branicky, Mar 05 2024
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Rick L. Shepherd, Jun 02 2005
EXTENSIONS
Offset corrected by Arkadiusz Wesolowski, Oct 18 2011
STATUS
approved