OFFSET
1,1
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000
FORMULA
a(n) ~ prime(n)^18.
EXAMPLE
115925123 is a term since 41^5 + 41^3 + 1 = 115925123 is prime.
MATHEMATICA
Select[ #^5 + #^3 + 1 & /@ Prime@ Range@ 90, PrimeQ] (* Robert G. Wilson v, Sep 10 2010 *)
PROG
(Python)
from itertools import islice
from sympy import isprime, nextprime
def agen(): # generator of terms
p = 2
while True:
t = p**5 + p**3 + 1
if isprime(t):
yield t
p = nextprime(p)
print(list(islice(agen(), 21))) # Michael S. Branicky, Mar 12 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Carmine Suriano, Sep 07 2010
EXTENSIONS
a(19) and beyond from Michael S. Branicky, Mar 12 2022
STATUS
approved