OFFSET
1,1
COMMENTS
Primes of the form (10^(d+1)+1) * x + 10^d for some x where 10^(d-1) <= x < 10^d.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 17117 is a term because it is prime and is the concatenation of 17, 1 and 17.
MAPLE
R:= NULL:
for d from 1 to 3 do
for x from 10^(d-1) to 10^d-1 do
y:= x * (1 + 10^(d+1))+10^d;
if isprime(y) then
R:= R, y;
fi
od od:
R;
PROG
(Python)
from sympy import isprime
from itertools import count, islice
def b(n): return 10**len(str(n))*(10*n+1) + n # A392225
def agen(): # generator of terms
yield from (filter(isprime, (b(n) for n in count(1) if n%10 in {1, 3, 7, 9})))
print(list(islice(agen(), 37))) # Michael S. Branicky, Jan 04 2026
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Robert Israel, Jan 03 2026
STATUS
approved
