OFFSET
1,1
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000
Patrick De Geest, World of Palindromic Primes
PROG
(Python)
from sympy import isprime
from itertools import count, islice, combinations_with_replacement as mc
def agen():
yield from (2, 3, 5, 7, 11)
for d in count(2):
ndni = (int("".join(m+m[:-1][::-1])) for m in mc("123456789", d))
yield from filter(isprime, ndni)
print(list(islice(agen(), 40))) # Michael S. Branicky, Jun 26 2022
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Patrick De Geest, Jun 07 2003
EXTENSIONS
a(39) and beyond from Michael S. Branicky, Jun 26 2022
STATUS
approved