OFFSET
1,1
COMMENTS
All terms == 1 (mod 6). The sequence is apparently infinite. There are 16486 terms up to 10^9. - Zak Seidov, Jan 17 2014
LINKS
Zak Seidov, Table of n, a(n) for n = 1..10000
EXAMPLE
53 is a term because 1531, 3533, 7537 and 9539 are primes.
PROG
(Python)
from sympy import isprime, nextprime
from itertools import islice
def agen(): # generator of terms
p = 2
while True:
sp = str(p)
if all(isprime(int(d+sp+d)) for d in "1379"):
yield p
p = nextprime(p)
print(list(islice(agen(), 40))) # Michael S. Branicky, Feb 23 2023
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Patrick De Geest, Feb 07 2001
STATUS
approved