OFFSET
1,2
COMMENTS
1 is a term if and only if number of terms in polynomial is prime.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..250
MATHEMATICA
a[n_]: = Module[{k = 1}, While[!PrimeQ[1 + k^(2*n+1) + Sum[k^(2*j), {j, 1, n}]], k++]; k]; Array[a, 30] (* Amiram Eldar, Mar 13 2020 *)
PROG
(PARI) a(n) = my(k = 1); while(! isprime(1 + k^(2*n+1) + sum(j=1, n, k^(2*j))), k++); k; \\ Michel Marcus, Mar 13 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Dec 31 2006
EXTENSIONS
More terms from Amiram Eldar, Mar 13 2020
STATUS
approved