A046137 - OEIS

%I #24 Apr 18 2022 11:35:14

%S 7,19,67,97,127,229,397,487,739,757,907,1009,1279,1447,1567,1597,1609,

%T 1867,1999,2239,2269,2377,2539,2659,2707,3037,3217,3319,3457,3529,

%U 3697,3877,3907,3919,4639,4789,4999,5167,5437,5569,5647,5689,5737

%N Primes p such that p+4 and p+12 are also prime.

%C All terms == 1 (mod 6). - _Robert Israel_, Jul 24 2015

%H Robert Israel, <a href="/A046137/b046137.txt">Table of n, a(n) for n = 1..10000</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PrimeTriplet.html">Prime Triplet</a>

%F A023200 INTERSECT A046133. - _R. J. Mathar_, Jan 23 2009

%p select(t -> isprime(t) and isprime(t+4) and isprime(t+12), [6*i+1 $ i=1..1000]); # _Robert Israel_, Jul 24 2015

%t Select[Range@ 5760, AllTrue[{#, # + 4, # + 12}, PrimeQ] &] (* _Michael De Vlieger_, Jul 24 2015, Version 10 *)

%t Select[Prime[Range[800]],AllTrue[#+{4,12},PrimeQ]&] (* _Harvey P. Dale_, Apr 18 2022 *)

%o (PARI) lista(nn) = forprime(p=2, nn, if (isprime(p+4) && isprime(p+12), print1(p, ", "))); \\ _Michel Marcus_, Jul 24 2015

%K nonn

%O 1,1

%A _Eric W. Weisstein_