OFFSET
1,1
COMMENTS
Primes in the sequence A236761.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
EXAMPLE
1213 is prime and 1213^4 - 1213 + 1 = 2164926732949 is prime. Thus, 1213 is a member of this sequence.
MATHEMATICA
Select[Prime[Range[10000]], PrimeQ[#^4 - # + 1]&] (* Vincenzo Librandi, Feb 14 2014 *)
PROG
(Python)
import sympy
from sympy import isprime
{print(n) for n in range(10**4) if isprime(n) and isprime(n**4-n+1)}
(PARI)
s=[]; forprime(p=2, 6000, if(isprime(p^4-p+1), s=concat(s, p))); s \\ Colin Barker, Feb 05 2014
(Magma) [p: p in PrimesUpTo(6000) | IsPrime(p^4-p+1)]; // Vincenzo Librandi, Feb 14 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Feb 01 2014
STATUS
approved