OFFSET
1,1
COMMENTS
a(17) > 5830. - Jinyuan Wang, Feb 05 2020
a(19) > 12000. - Michael S. Branicky, Jul 04 2024
EXAMPLE
k = 3 is here because 7*3! + 1 = 43 is prime.
PROG
(PARI) is(k) = ispseudoprime(7*k!+1); \\ Jinyuan Wang, Feb 05 2020
(Python)
from sympy import isprime
from math import factorial
def aupto(m): return [k for k in range(m+1) if isprime(7*factorial(k)+1)]
print(aupto(300)) # Michael S. Branicky, Mar 07 2021
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Phillip L. Poplin (plpoplin(AT)bellsouth.net), Oct 25 2002
EXTENSIONS
a(17)-a(18) from Michael S. Branicky, Jul 04 2024
STATUS
approved