A051888 - OEIS

A051888

a(n) is the smallest prime p such that p*n! + 1 is prime.

2, 2, 2, 2, 3, 2, 3, 3, 7, 3, 3, 5, 2, 3, 13, 7, 31, 5, 2, 7, 17, 67, 41, 3, 13, 3, 43, 17, 97, 7, 29, 109, 3, 71, 5, 2, 7, 41, 3, 59, 3, 11, 29, 7, 107, 67, 79, 3, 743, 149, 163, 2, 211, 2, 19, 71, 73, 23, 37, 113, 149, 67, 41, 617, 107, 37, 107, 283, 113, 19, 239, 107, 73, 97, 5

(list; graph; refs; listen; history; text; internal format)

OFFSET

0,1

COMMENTS

Analogous to or subset of A051686; generalization of A005384.

The PFGW program has been used to certify all the primes corresponding to the terms up to a(1000), using a deterministic test which exploits the factorization of a(n) - 1. - Giovanni Resta, May 30 2018

LINKS

Giovanni Resta, Table of n, a(n) for n = 0..1000

FORMULA

a(n) = (A051901(n)-1)/n!. - Amiram Eldar, Feb 25 2025

MATHEMATICA

Do[k = 1; While[ !PrimeQ[ Prime[k]*n! + 1], k++ ]; Print[ Prime[k]], {n, 1, 75} ]

spp[n_]:=Module[{p=2, nf=n!}, While[!PrimeQ[p*nf+1], p=NextPrime[p]]; p]; Array[ spp, 80, 0] (* Harvey P. Dale, May 17 2019 *)

PROG

(PARI) a(n) = {my(p=2); while (!isprime(p*n! + 1), p = nextprime(p+1)); p; } \\ Michel Marcus, May 28 2018

CROSSREFS

Cf. A005384, A051686, A051886, A051887, A051901, A064983.

Sequence in context: A070098 A029233 A147981 * A377014 A327966 A366580