OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..5000
EXAMPLE
a(1)=3 because 3^4-(3+1)^3=81-64=17 is prime.
MATHEMATICA
lst={}; Do[If[PrimeQ[n^4-(n+1)^3], AppendTo[lst, n]], {n, 3, 1000}]; lst (* Vincenzo Librandi, Jul 16 2012 *)
PROG
(Magma) [n: n in [3..300]| IsPrime(n^4-(n+1)^3)]; // Vincenzo Librandi, Jul 16 2012
(PARI) is(n)=isprime(n^4-(n+1)^3) \\ Charles R Greathouse IV, Jun 05 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Hugo Pfoertner, Aug 23 2003
STATUS
approved