OFFSET
1,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
MAPLE
A := {}; for n to 1000 do p := ithprime(n); if isprime(p^3+p+1) then A := `union`(A, {p}) end if end do; A := A
MATHEMATICA
Select[Prime[Range[400]], PrimeQ[#^3 + # + 1] &] (* Bruno Berselli, Sep 01 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(2000) | IsPrime(p^3+p+1)]; // Bruno Berselli, Sep 01 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
César Eliud Lozada, Aug 31 2012
STATUS
approved