OFFSET
1,1
LINKS
Daniel Starodubtsev, Table of n, a(n) for n = 1..10000
MATHEMATICA
a = {}; Do[If[PrimeQ[1 + n^2 + n^4 + n^6 + n^8 + n^10 + n^12 + n^13], AppendTo[a, n]], {n, 1, 1400}]; a
PROG
(PARI) is(n)=isprime(1+n^2+n^4+n^6+n^8+n^10+n^12+n^13) \\ Charles R Greathouse IV, Jun 13 2017
(Python)
from sympy import isprime
def ok(k): return isprime(1+sum(k**i for i in [2, 4, 6, 8, 10, 12, 13]))
print([k for k in range(2000) if ok(k)]) # Michael S. Branicky, Oct 24 2021
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Artur Jasinski, Dec 31 2006
STATUS
approved