OFFSET
1,1
LINKS
K. D. Bajpai, Table of n, a(n) for n = 1..3379
EXAMPLE
12 appears in the sequence because:
12^2 + 2 = 146 = 2*73
12^2 - 2 = 142 = 2*71
12 + 2 = 14 = 2*7
12 - 2 = 10 = 2*5 are all semiprimes.
MAPLE
with(numtheory): select(n -> (bigomega(n^2 + 2)=2 and bigomega(n^2 - 2)=2 and bigomega(n + 2)=2 and bigomega(n - 2)=2), [seq(n, n=1..10000)]);
MATHEMATICA
Select[Range[10000], PrimeOmega[#^2 + 2] == PrimeOmega[#^2 - 2] == PrimeOmega[# + 2] == PrimeOmega[# - 2] == 2 &]
PROG
(PARI) for(n = 1, 10000, if(bigomega(n^2 + 2) == 2 && bigomega(n^2 - 2) == 2 && bigomega(n + 2) == 2 && bigomega(n - 2) == 2, print1(n, ", ")))
(Magma) IsSemiprime:=func<i | &+[d[2]: d in Factorization(i)] eq 2>; [ n : n in [2..10000] | IsSemiprime(n^2 + 2) and IsSemiprime(n^2 - 2) and IsSemiprime(n + 2) and IsSemiprime(n - 2)];
CROSSREFS
KEYWORD
nonn,less
AUTHOR
K. D. Bajpai, Jan 28 2016
STATUS
approved