OFFSET
1,2
COMMENTS
Since 8*x^3 - 8*x^2 + 4*x + 1 = (2*x-1)^3 + (2x-1)^2 + (2x-1) + 2, the corresponding primes are of the form k^3 + k^2 + k + 2 and are in A088547.
LINKS
Jens Kruse Andersen, Table of n, a(n) for n = 1..10000
FORMULA
8*a(n)^3 - 8*a(n)^2 + 4*a(n) + 1 = A088547(n+1).
MATHEMATICA
Select[Range[1000], PrimeQ[1 + 4*# - 8*#^2 + 8*#^3] &]
PROG
(PARI) is(n)=isprime(8*n^3-8*n^2+4*n+1)
CROSSREFS
KEYWORD
nonn
AUTHOR
Jean-Christophe Hervé, Sep 17 2014
STATUS
approved