%I #16 Oct 18 2020 03:13:12
%S 12,18,20,24,28,40,44,45,48,50,52,54,56,60,63,68,72,75,76,80,84,88,90,
%T 92,96,98,99,104,108,112,116,117,120,124,126,132,135,136,140,144,147,
%U 148,150,152,153,156,160,162,164,168,171,172,175,176,180,184,188
%N Numbers with exactly two distinct exponents in their prime factorization, or two distinct parts in their prime signature.
%C The first term is A006939(2) = 12.
%C First differs from A059404 in lacking 360, whose prime signature has three distinct parts.
%C Positions of 2's in A071625.
%C Numbers k such that A001221(A181819(k)) = 2.
%C The asymptotic density of this sequence is (6/Pi^2) * Sum_{n>=2, n squarefree} 1/((n-1)*psi(n)) = 0.3611398..., where psi is the Dedekind psi function (A001615) (Sanna, 2020). - _Amiram Eldar_, Oct 18 2020
%H Amiram Eldar, <a href="/A323055/b323055.txt">Table of n, a(n) for n = 1..10000</a>
%H Carlo Sanna, <a href="https://doi.org/10.1007/s12044-020-0556-y">On the number of distinct exponents in the prime factorization of an integer</a>, Proceedings - Mathematical Sciences, Indian Academy of Sciences, Vol. 130, No. 1 (2020), Article 27, <a href="https://www.ias.ac.in/describe/article/pmsc/130/0027">alternative link</a>.
%e 3000 = 2^3 * 3^1 * 5^3 has two distinct exponents {1, 3}, so belongs to the sequence.
%p isA323055 := proc(n)
%p local eset;
%p eset := {};
%p for pf in ifactors(n)[2] do
%p eset := eset union {pf[2]} ;
%p end do:
%p simplify(nops(eset) = 2 ) ;
%p end proc:
%p for n from 12 to 1000 do
%p if isA323055(n) then
%p printf("%d,",n) ;
%p end if;
%p end do: # _R. J. Mathar_, Jan 09 2019
%t Select[Range[100],Length[Union[Last/@FactorInteger[#]]]==2&]
%Y One distinct exponent: A062770 or A072774.
%Y Two distinct exponents: A323055.
%Y Three distinct exponents: A323024.
%Y Four distinct exponents: A323025.
%Y Five distinct exponents: A323056.
%Y Cf. A001221, A001222, A001615, A006939, A007774, A059404, A071625, A118914, A181819, A323014, A323022.
%K nonn
%O 1,1
%A _Gus Wiseman_, Jan 03 2019
|