OFFSET
1,1
COMMENTS
The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 0, 0, 1, 4, 57, 555, 5492, 55078, 551443, 5512825, ... . Apparently, the asymptotic density of this sequence exists and equals 0.000551... . - Amiram Eldar, Jan 18 2023
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..5000 from Michael De Vlieger)
EXAMPLE
k = 844 is a term since 844 = 2^2*211, k+1 = 845 = 5*13^2, k+2 = 846 = 2*3^2*47, and k+4 = 847 = 7*11^2.
MATHEMATICA
With[{s = Select[Range[10^5], And[MemberQ[#, 2], FreeQ[#, k_ /; k > 2]] &@ FactorInteger[#][[All, -1]] &]}, Function[t, Part[s, #] &@ SequencePosition[t, {1, 1, 1}][[All, 1]]]@ Differences@ s] (* Michael De Vlieger, Jul 30 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, May 29 2002
STATUS
approved