OFFSET
1,1
COMMENTS
Primes p such that both p and p+2 are terms in A086006.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Zak Seidov, First 100 terms.
MATHEMATICA
Select[Prime[Range[10^6]], PrimeQ[# + 2] && PrimeOmega[2*# - 1] == PrimeOmega[2*# + 1] == 2 && PrimeOmega[2*# + 3] == PrimeOmega[2*# + 5] == 2 &] (* Amiram Eldar, Apr 06 2023 *)
PROG
(PARI) is(p) = (bigomega(2*p+1) == 2) && (bigomega(2*p-1) == 2);
lista(max) = {my(q1 = 0, q2, prev = 0); forprime(p = 2, max, q2 = is(p); if(p == prev + 2 && q1 && q2, print1(prev, ", ")); prev = p; q1 = q2); } \\ Amiram Eldar, Apr 06 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Mar 20 2009
EXTENSIONS
More terms from Amiram Eldar, Apr 06 2023
STATUS
approved