OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..10000
Ken Ono and Scott Ahlgren, Weierstrass points on X0(p) and supersingular j-invariants, Mathematische Annalen 325, 2003, pp. 355-368.
MATHEMATICA
Select[Table[n^3-n-1, {n, 0, 800}], PrimeQ] (* Vincenzo Librandi, Dec 07 2011 *)
PROG
(Magma) [ a: n in [1..200] | IsPrime(a) where a is n^3-n-1 ]; // Vincenzo Librandi, Dec 07 2011
(Python)
from sympy import isprime
def aupton(terms):
k, alst = 2, []
while len(alst) < terms:
if isprime(k**3-k-1): alst.append(k**3-k-1)
k += 1
return alst
print(aupton(37)) # Michael S. Branicky, May 23 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Roger L. Bagula, Mar 22 2006
EXTENSIONS
Edited by N. J. A. Sloane, Jan 01 2007
More terms from Artur Jasinski, Jan 01 2007
STATUS
approved