%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