A346167 Primes p such that p*p! + 1 is also prime.
2, 3, 5, 7, 19, 23, 149, 151, 197, 37691
Offset: 1
Programs
-
Maple
select(p -> isprime(p) and isprime(p*factorial(p) + 1), [$2 .. 200]);
-
Mathematica
Select[Range[2, 200], PrimeQ[#] && PrimeQ[#*#! + 1] &] Select[Prime[Range[100]],PrimeQ[#*#!+1]&] (* Harvey P. Dale, Mar 21 2025 *)
-
PARI
a = List(); for(p=2, 200, if(isprime(p) && isprime(p*p!+1), listput(a, p))); a
-
Sage
[p for p in range(2, 200) if is_prime(p) and is_prime(p*factorial(p) + 1)]
Extensions
a(10) from Georg Grasegger, Apr 07 2025