OFFSET
1,1
COMMENTS
Sophie Germain primes of the form p*q*r - 1, where p, q and r are primes.
Except for 11, all terms == 5 (mod 12).
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 41 is a term because it is prime, 41 + 1 = 42 = 2 * 3 * 7 is a triprime, and 41 * 2 + 1 = 83 is prime.
MAPLE
select(p -> isprime(p) and isprime(2*p+1) and numtheory:-bigomega(p+1) = 3, [seq(i, i=3..30000, 2)]);
MATHEMATICA
s= {}; Do[p = Prime[k]; If[3 == PrimeOmega[p + 1] && PrimeQ[2*p +1], AppendTo[s, p]], {k, 2000}]; s
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov and Robert Israel, Jul 17 2025
STATUS
approved
