OFFSET
1,1
COMMENTS
Primes in A335233.
MATHEMATICA
Select[Prime@ Range@ 100, (f = #!; NoneTrue[f*Range[#] + 1, PrimeQ ]) &] (* Robert Price, Sep 14 2020 *)
PROG
(Python)
from sympy import isprime, nextprime, factorial
A335361_list, p = [], 2
while p < 500:
f, g = factorial(p), 1
for i in range(1, p+1):
g += f
if isprime(g):
break
else:
A335361_list.append(p)
p = nextprime(p)
(PARI) is(p) = if(isprime(p), for(i=1, p, if(ispseudoprime(i*p!+1), return(0))); 1, 0); \\ Jinyuan Wang, Jun 21 2020
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Chai Wah Wu, Jun 10 2020
EXTENSIONS
a(34) from Jinyuan Wang, Jun 21 2020
STATUS
approved