OFFSET
1,1
COMMENTS
Also primes of the form 4*k^3 + 4*k^2 - 1.
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..3240
EXAMPLE
k=15: (15^3 - 15^2 - 15 - 1)/2 = 1567 (is prime).
MATHEMATICA
Select[Table[(n^3 - n^2 - n - 1) / 2, {n, 200}], PrimeQ] (* Vincenzo Librandi, Jan 26 2016 *)
PROG
(Sage) [(k^3-k^2-k-1)/2 for k in [2*i+1 for i in [1..100]] if is_prime(Integer((k^3-k^2-k-1)/2))] # Tom Edgar, Jan 25 2016
(Magma) [a: n in [0..200] | IsPrime(a) where a is (n^3-n^2-n-1) div 2 ]; // Vincenzo Librandi, Jan 26 2016
(PARI) lista(nn) = for(n=1, nn, if(ispseudoprime(p=4*n^3+4*n^2-1), print1(p, ", "))); \\ Altug Alkan, Mar 14 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Emre APARI, Jan 25 2016
EXTENSIONS
More terms from Tom Edgar, Jan 25 2016
STATUS
approved