|
%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_
|
