Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #45 Sep 08 2022 08:44:56
%S 11,13,17,19,23,29,37,43,47,53,59,67,73,79,89,103,107,109,113,137,157,
%T 163,173,179,197,199,229,233,239,257,263,269,277,283,313,317,337,353,
%U 359,373,379,389,439,449,463,467,509,547,563,569,577,593,599,607,613
%N Primes p such that p-6 is also prime. (Upper of a pair of sexy primes.)
%H T. D. Noe, <a href="/A046117/b046117.txt">Table of n, a(n) for n = 1..1000</a>
%H Ernest G. Hibbs, <a href="https://www.proquest.com/openview/4012f0286b785cd732c78eb0fc6fce80">Component Interactions of the Prime Numbers</a>, Ph. D. Thesis, Capitol Technology Univ. (2022), see p. 33.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/SexyPrimes.html">Sexy Primes</a>. [The definition in this webpage is unsatisfactory, because it defines a "sexy prime" as a pair of primes.- _N. J. A. Sloane_, Mar 07 2021].
%F a(n) = A087695(n) + 3.
%F a(n) = A023201(n) + 6. - _M. F. Hasler_, Jan 02 2020
%t q=6; a={}; Do[p1=Prime[n]; p2=p1+q; If[PrimeQ[p2], AppendTo[a, p2]], {n, 7^2}]; a "and/or" Select[Prime[Range[3, 7^2]], PrimeQ[ # ]&&PrimeQ[ #-6]&] (* _Vladimir Joseph Stephan Orlovsky_, Aug 07 2008 *)
%t Select[Prime[Range[4,200]],PrimeQ[#-6]&] (* _Harvey P. Dale_, Mar 31 2018 *)
%o (PARI) forprime(p=2,1e4,if(isprime(p-6),print1(p", "))) \\ _Charles R Greathouse IV_, Jul 15 2011
%o (Magma) [p:p in PrimesInInterval(7,650)| IsPrime(p-6)]; // _Marius A. Burtea_, Jan 03 2020
%Y Cf. A023201, A098428, A087695.
%K nonn
%O 1,1
%A _Eric W. Weisstein_
%E Name edited by _M. F. Hasler_, Jan 02 2020