OFFSET
1,1
COMMENTS
For odd numbers p, q, r, p*q*r - p - q - r == 2 (mod 4).
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 23459 = (31*37*41 - 31 - 37 - 41)/2 where 31, 37 and 41 are consecutive primes, and 23459 is prime.
MAPLE
P:= select(isprime, [seq(i, i=3..10^4, 2)]):
select(isprime, [seq((P[i]*P[i+1]*P[i+2]-P[i]-P[i+1]-P[i+2])/2, i=1..nops(P)-2)]);
MATHEMATICA
s={}; Do[p=Prime[j]; q=NextPrime[p]; r=NextPrime[p, 2]; If[PrimeQ[(p*q*r-p-q-r)/2], AppendTo[s, (p*q*r-p-q-r)/2]], {j, 216}]; s (* James C. McMahon, Nov 25 2025 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Will Gosnell and Robert Israel, Nov 19 2025
STATUS
approved
