OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
EXAMPLE
42181 is a term because 42181 and 42187 are sexy primes while 42182-42186 each have 4 prime divisors (counting multiplicity).
MATHEMATICA
Select[Range[2*10^6], AllTrue[{#, #+6}, PrimeQ]&&Length[Union[ PrimeOmega[ Range[ #+1, #+5]]]]==1&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Dec 16 2015 *)
PROG
(Magma) f:=func<n|&+[p[2]: p in Factorization(n)]>; [p:p in PrimesUpTo(2000000)| IsPrime(p+6) and forall{k:k in [2..5]|f(p+k) eq f(p+1)} ]; // Marius A. Burtea, Dec 16 2019
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jason Earls, May 29 2004
STATUS
approved