OFFSET
1,1
LINKS
John Cerkan, Table of n, a(n) for n = 1..10000
EXAMPLE
4^8 - 4^7 - 4^6 - 4^5 - 4^4 - 4^3 - 4^2 - 4^1 - 1 = 43691 is prime. Thus 4 is a term.
MATHEMATICA
Rest@ Select[Range@ 450, Function[n, PrimeQ[Fold[#1 - n^#2 &, n^8, Range@ 7] - 1]]] (* Michael De Vlieger, Apr 03 2017 *)
PROG
(Python)
import sympy
from sympy import isprime
{print(n, end=', ') for n in range(10**3) if isprime(n**8-n**7-n**6-n**5-n**4-n**3-n**2-n-1)}
(PARI) for(n=1, 10^3, if(ispseudoprime(n^8-sum(i=0, 7, n^i)), print1(n, ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Jun 02 2014
STATUS
approved