OFFSET
1,1
LINKS
Pierre CAMI, Table of n, a(n) for n = 1..500
EXAMPLE
30030*j-1 or 30030*j+1 are not both primes for j=1,2,3,4,5. But for j=6 {180179,180181} are twin primes. So a(6)=6.
MATHEMATICA
smp[n_]:=Module[{k=1}, While[!PrimeQ[k*n+1]||!PrimeQ[k*n-1], k++]; k]; Table[ smp[n], {n, FoldList[Times, Prime[Range[70]]]}] (* Harvey P. Dale, Oct 27 2016 *)
PROG
(PARI) a(n)=p=vecprod(primes(n)); for(k=1, +oo, ispseudoprime(k*p+1)&&ispseudoprime(k*p-1)&&return(k)) \\ Jeppe Stig Nielsen, Nov 09 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Mar 22 2001
EXTENSIONS
Corrected and extended by Ray Chandler, Apr 03 2009
STATUS
approved