|
|
|
|
1, 2, 5, 6, 10, 13, 14, 21, 22, 29, 30, 33, 34, 37, 38, 41, 42, 44, 46, 49, 57, 58, 61, 65, 66, 69, 70, 73, 75, 77, 78, 80, 82, 85, 86, 93, 94, 98, 101, 102, 105, 106, 109, 110, 113, 114, 116, 118, 122, 129, 130, 133, 135, 137, 138, 141, 142, 145, 147, 154, 157
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
First differs from A358817 at n = 165.
First differs from A369166 at n = 558. a(558) = 1520 is the least term that is not in A369166. A369166(144273) = 397952 is the least term of A369166 that is not a term of this sequence.
The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 5, 38, 368, 3638, 36337, 363163, 3631569, 36315800, 363156839, 3631559150, ... . Apparently, the asymptotic density of this sequence exists and equals 0.36315... .
|
|
LINKS
|
|
|
MATHEMATICA
|
s[n_] := s[n] = Times @@ FactorInteger[n][[;; , 2]]; Select[Range[300], s[#] == s[# + 1] &]
|
|
PROG
|
(PARI) lista(kmax) = {my(c1 = 1, c2); for(k = 2, kmax, c2 = vecprod(factor(k)[, 2]); if(c1 == c2, print1(k-1, ", ")); c1 = c2); }
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|