OFFSET
1,1
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..1000
MATHEMATICA
w = {3};
Do[k = 1;
q = Monitor[
Parallelize[
While[True,
If[PrimeQ[FromDigits[
Join @@ IntegerDigits /@ Append[w, Prime[k]]]], Break[]]; k++];
Prime[k]], k];
w = Append[w, q], {i, 2, 57}];
w
PROG
(Python)
from itertools import count, islice
from gmpy2 import digits, is_prime, mpz, next_prime
def agen(): # generator of terms
s, an = "", 3
while True:
yield int(an)
s += digits(an)
p = 3
while not is_prime(mpz(s+digits(p))): p = next_prime(p)
an = p
print(list(islice(agen(), 57))) # Michael S. Branicky, Dec 21 2024
CROSSREFS
KEYWORD
base,nonn
AUTHOR
J.W.L. (Jan) Eerland, Dec 21 2024
STATUS
approved