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

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

Offset

1,1

Comments

a(10) > 10^4.

Links

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

Maple

select(p -> isprime(p) and isprime(p*factorial(p) + 1), [$2 .. 200])

Mathematica

Select[Range[2, 200], PrimeQ[#] && PrimeQ[#*#! + 1] &]

Prog

(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)]

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

Status

approved