A158476 Lesser of twin pairs in A086006.

107, 71549, 97379, 103067, 136709, 342449, 403829, 686969, 695879, 758267, 795797, 805499, 887399, 941489, 945881, 1023227, 1025747, 1081709, 1081979, 1169027, 1179989, 1261259, 1456919, 1554101, 2110877, 2148659, 2167469, 2433059, 2494439, 2554397, 2623571, 2803121

Offset: 1

Author

Zak Seidov, Mar 20 2009

Keywords

• nonn

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.

Crossrefs

Cf. A086005, A086006.

Programs

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 *)

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

Extensions

More terms from Amiram Eldar, Apr 06 2023