|
%I #18 Sep 08 2022 08:45:12
%S 3,11,13,19,31,41,59,71,73,101,139,193,239,269,271,313,349,433,521,
%T 643,823,829,881,1051,1061,1093,1621,1669,1723,1951,2549,2999,3359,
%U 3373,3463,3469,3583,4019,4219,4481,4483,4519,5233,5639,5881,6089,6131,6133
%N Primes p such that 6*p-7 and 6*p-5 are twin primes and p is also a twin prime.
%H Robert Israel, <a href="/A089152/b089152.txt">Table of n, a(n) for n = 1..5462</a>
%e 6*19-7 = 107, 6*19-5 = 109, 107 and 109 are twin primes and 19 has 17 as twin.
%p P:= {seq(ithprime(i),i=1..20000)}:
%p A:= P intersect (map(`-`,P,2) union map(`+`,P,2)) intersect map(t -> (t+5)/6, P) intersect map(t -> (t+7)/6, P):
%p sort(convert(A,list)); # _Robert Israel_, Jan 16 2019
%o (Magma) [p: p in PrimesUpTo(6500) |(IsPrime(p-2) or IsPrime(p+2)) and IsPrime(6*p-5) and IsPrime(6*p-7)]; // _Vincenzo Librandi_, Jan 17 2019
%Y Subsequence of A089151.
%K easy,nonn
%O 1,1
%A _Pierre CAMI_, Dec 06 2003
%E 2 and 173 removed by _R. J. Mathar_, Nov 13 2009
|
