%I #34 Mar 18 2021 05:53:25
%S 7,11,17,23,41,47,71,83,101,107,113,131,167,197,227,233,281,311,317,
%T 353,383,401,443,461,467,491,503,617,647,677,743,761,773,827,857,863,
%U 881,887,911,941,971,1013,1091,1097,1217,1283,1301,1307,1427,1433
%N Larger member p+4 of cousin primes (p, p+4).
%C A pair of cousin primes are primes of the form p and p+4 (where p+2 may or may not be a prime). - _N. J. A. Sloane_, Mar 18 2021
%H Reinhard Zumkeller, <a href="/A046132/b046132.txt">Table of n, a(n) for n = 1..10000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/CousinPrimes.html">Cousin Primes</a>
%F a(n) = A023200(n) + 4 = A087679(n) + 2.
%F a(n) = 3*A157834(n-1) + 2 = A029710(n-1) + 4 = 6*A056956(n-1) + 5 (thus a(n) mod 6 == 5), for all n>1. - _M. F. Hasler_, Jan 15 2013
%t lst={};Do[p=Prime[n];If[PrimeQ[p4=p+4], (*Print[p4];*)AppendTo[lst, p4]], {n, 10^2}];lst (* _Vladimir Joseph Stephan Orlovsky_, Aug 21 2008 *)
%t Select[Prime[Range[300]],PrimeQ[#+4]&]+4 (* _Harvey P. Dale_, Dec 15 2017 *)
%o (PARI) forprime(p=2,1e5,if(isprime(p-4),print1(p", "))) \\ _Charles R Greathouse IV_, Jul 15 2011
%o (Haskell)
%o a046132 n = a046132_list !! (n-1)
%o a046132_list = filter ((== 1) . a010051') $ map (+ 4) a000040_list
%o -- _Reinhard Zumkeller_, Aug 01 2014
%Y Essentially the same as A031505. Cf. A023200, A029710, A098429.
%Y Cf. A000040, A010051.
%K nonn
%O 1,1
%A _Eric W. Weisstein_
|