OFFSET
1,1
COMMENTS
a(n) = p+q*(r+s) where p = A343448(n) and q,r,s are the next three primes after p.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
For n = 3, A343448(3) = 11 and the next three primes are 13, 17, 19, so a(3) = 11+13*(17+19) = 479.
MAPLE
B:= NULL: q:= 2: r:= 3: s:= 5: count:= 0:
while count < 100 do
p:= q; q:= r; r:= s; s:= nextprime(s);
v:= p+q*(r+s);
if isprime(v) then B:= B, v; count:= count+1 fi
od:
B;
PROG
(Python)
from sympy import isprime, nextprime
def aupto(limit):
p, q, r, s, alst = 2, 3, 5, 7, []
t = p + q*(r+s)
while t <= limit:
if isprime(t): alst.append(t)
p, q, r, s = q, r, s, nextprime(s)
t = p + q*(r+s)
return alst
print(aupto(2067529)) # Michael S. Branicky, Apr 15 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Apr 15 2021
STATUS
approved