OFFSET
1,1
COMMENTS
All terms are odd.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 77 = 7*11 is a term because p = 7, q = 11, (7*11+1+11-7)/2 = 41 and (7*11+1+7-11)/2 = 37 are prime.
MAPLE
filter:= proc(n) local p, q, t, k;
if issqr(n) or numtheory:-bigomega(n) <> 2 then return false fi;
p, q:= op(numtheory:-factorset(n));
isprime((n+1+p-q)/2) and isprime((n+1+q-p)/2)
end proc:
select(filter, [seq(i, i=3..10000, 2)]);
MATHEMATICA
Select[Range[1, 1400, 2], (f = FactorInteger[#])[[;; , 2]] == {1, 1} && PrimeQ[((p = f[[1, 1]])*(q = f[[2, 1]]) + 1 + p - q)/2] && PrimeQ[(p*q + 1 + q - p)/2] &] (* Amiram Eldar, Aug 16 2022 *)
sfspQ[n_]:=Module[{p, q}, {p, q}=FactorInteger[n][[;; , 1]]; AllTrue[{(n+1+p-q)/2, (n+1+q-p)/2}, PrimeQ]]; Select[Range[1500], SquareFreeQ[#]&&PrimeOmega[#]==2&&sfspQ[#]&] (* Harvey P. Dale, Mar 01 2024 *)
CROSSREFS
KEYWORD
nonn,less
AUTHOR
J. M. Bergot and Robert Israel, Aug 15 2022
STATUS
approved