OFFSET
1,1
COMMENTS
The asymptotic density of this sequence is 1 - 2/zeta(5) + Product_{p prime} (1 - 2/p^5) = 0.000284512101137896862... - Amiram Eldar, Feb 16 2021
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Harvey P. Dale)
MATHEMATICA
Select[ Range[2, 250000], Max[ Transpose[ FactorInteger[ # ]] [[2]]] > 4 && Max[ Transpose[ FactorInteger[ # + 1]] [[2]]] > 4 &]
SequencePosition[Table[If[Max[FactorInteger[n][[All, 2]]]>4, 1, 0], {n, 120000}], {1, 1}][[All, 1]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Dec 25 2018 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Mar 04 2002
STATUS
approved