OFFSET
1,2
COMMENTS
k = {1, 2, 4, 67} yields primes p(k) = {2, 3, 7, 331}. There are no more such k up to k=100. Computed in collaboration with Ray Chandler.
Terms a(5) and greater are only probable primes. - Iain Fox, Mar 05 2018
LINKS
J. V. Post, Math Pages.
FORMULA
Numbers k such that (prime(k)-1)! + prime(k)^4 is prime, where prime(k) is the k-th prime.
EXAMPLE
a(3) = 4 because (prime(4)-1)! + prime(4)^4 = (7-1)! + 7^4 = 720 + 2401 = 3121 is the 3rd prime of this form.
MATHEMATICA
lst={}; Do[p=Prime[n]; If[PrimeQ[(p-1)!+p^4], AppendTo[lst, n]], {n, 10^2}]; lst (* Vladimir Joseph Stephan Orlovsky, Sep 08 2008 *)
PROG
(PARI) is(k) = my(p=prime(k)); ispseudoprime((p-1)! + p^4) \\ Iain Fox, Mar 05 2018
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Jonathan Vos Post, Nov 30 2004
STATUS
approved