A346167 Primes p such that p*p! + 1 is also prime.

2, 3, 5, 7, 19, 23, 149, 151, 197, 37691

Offset

1,1

Links

Table of n, a(n) for n=1..10.

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 *)

Prog

(PARI) a = List(); for(p=2, 200, if(isprime(p) && isprime(p*p!+1), listput(a, p))); a
(SageMath) [p for p in range(2, 200) if is_prime(p) and is_prime(p*factorial(p) + 1)]

Crossrefs

Prime terms of A090703.

Cf. A346168.

Sequence in context: A345335 A025019 A140327 * A163074 A230041 A068803

Adjacent sequences: A346164 A346165 A346166 * A346168 A346169 A346170

Keyword

nonn,more

Author

Reza K Ghazi, Jul 08 2021

Extensions

a(10) from Georg Grasegger, Apr 07 2025

Status

approved