OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
a(1) = 3 because for k = 1, 1*(1^3) - 1 = 0 is not prime, for k = 2, 2*(1^3) - 1 = 1 is not prime, but for k = 3, 3*(1^3) - 1 = 2 is prime.
MATHEMATICA
sk[n_]:=Module[{k=1, n3=n^3}, While[!PrimeQ[k*n3-1], k++]; k]; Array[sk, 80] (* Harvey P. Dale, Jan 04 2023 *)
PROG
(Python)
import sympy
from sympy import isprime
def f(n):
..for k in range(1, 10**3):
....if isprime(k*(n**3)-1):
......return k
n = 1
while n < 10**3:
..print(f(n))
..n += 1
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Mar 06 2014
STATUS
approved