OFFSET
1,1
COMMENTS
If you graph a(n) versus n, a clear pattern emerges.
As you go farther along the n-axis, greater are the number of consecutive sexy primes, on average, within each interval obtained.
If one could prove that there is at least one consecutive sexy prime within each interval, this would imply that consecutive sexy primes are infinite.
I suspect all numbers in the sequence are > 0.
LINKS
J. S. Cheema, Table of n, a(n) for n = 1..1762
Rick Aster, Prime number sieve SAS prime sieve program
Eric Weisstein's World of Mathematics, Sexy Primes. [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].
Wikipedia, Sexy Primes
EXAMPLE
The first sexy prime pair with consecutive primes is (23,29) = A031924(1) and A031925(1). Square the first term, you get 529, then take the product of the two primes, you get 667.
Between these two numbers, namely (529,667), there are ten consecutive sexy primes: (541,547), (557,563), (563,569),
(571,577), (587,593), (593,599), (601,607), (607,613), (647,653), and (653 659).
Hence the very first term of the sequence is 10.
MAPLE
isA031924 := proc(p) return (isprime(p) and (nextprime(p)-p) = 6 ); end proc:
A031924 := proc(n) local p; if n = 1 then 23; else p := nextprime(procname(n-1)) ; while not isA031924(p) do p := nextprime(p) ; end do ; return p ; end if ; end proc:
A173198 := proc(n) local ulim, llim, a, i ; llim := A031924(n)^2 ; ulim := A031924(n)*A031925(n) ; a := 0 ; for i from llim to ulim-6 do if isA031924(i) then a := a+1 ; end if; end do ; a ; end proc:
seq(A173198(n), n=1..80) ; # R. J. Mathar, Feb 15 2010
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaspal Singh Cheema, Feb 12 2010
EXTENSIONS
Comments condensed by R. J. Mathar, Feb 15 2010
STATUS
approved