OFFSET
1,1
COMMENTS
Positions of 4's in A071625.
Is a(n) ~ c * n for some c? - David A. Corneth, Jan 09 2019
The asymptotic density of this sequence is (6/Pi^2) * Sum_{n>=2, n squarefree} r(n)/((n-1)*psi(n)) = 0.00035750... (corresponding to c = 2797.1... in the question above, whose answer is affirmative), where psi is the Dedekind psi function (A001615), and r(n) = Sum_{d_1|n, 1<d_1<n} (1/(d_1-1)) * Sum_{d_2|d_1, 1<d_2<d_1} 1/(d_2-1) (Sanna, 2020). - Amiram Eldar, Oct 18 2020
LINKS
David A. Corneth, Table of n, a(n) for n = 1..10000
Carlo Sanna, On the number of distinct exponents in the prime factorization of an integer, Proceedings - Mathematical Sciences, Indian Academy of Sciences, Vol. 130, No. 1 (2020), Article 27, alternative link.
EXAMPLE
126000 = 2^4 * 3^2 * 5^3 * 7^1 has four distinct exponents {1, 2, 3, 4}, so belongs to the sequence.
831600 = 2^4 * 3^3 * 5^2 * 7^1 * 11^1 has four distinct exponents {1, 2, 3, 4}, so belongs to the sequence.
MATHEMATICA
tom[n_]:=Length[Union[Last/@If[n==1, {}, FactorInteger[n]]]];
Select[Range[100000], tom[#]==4&]
PROG
(PARI) is(n) = #Set(factor(n)[, 2]) == 4 \\ David A. Corneth, Jan 09 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jan 02 2019
STATUS
approved