|
%I #15 Aug 15 2017 20:34:21
%S 1,1,1,2,1,1,1,3,2,1,1,2,1,1,1,4,1,2,1,2,1,1,1,3,2,1,3,2,1,1,1,5,1,1,
%T 1,2,1,1,1,3,1,1,1,2,2,1,1,4,2,2,1,2,1,3,1,3,1,1,1,2,1,1,2,6,1,1,1,2,
%U 1,1,1,6,1,1,2,2,1,1,1,4,4,1,1,2,1,1,1,3,1,2,1,2,1,1,1,5,1,2,2,2,1,1,1,3,1,1,1,6,1,1,1,4,1,1,1,2,2,1,1,3
%N a(1) = 1; for n > 1, a(n) = product of distinct exponents in the prime factorization of n.
%H Antti Karttunen, <a href="/A290107/b290107.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Eu#epf">Index entries for sequences computed from exponents in factorization of n</a>
%F a(n) = A156061(A181819(n)).
%e For n = 36 = 2^2 * 3^2, the only distinct exponent that occurs is 2, thus a(36) = 2.
%e For n = 144 = 2^4 * 3^2, the distinct exponents are 2 and 4, thus a(144) = 2*4 = 8.
%e For n = 4500 = 2^2 * 3^2 * 5^3, the distinct exponents are 2 and 3, thus a(4500) = 2*3 = 6.
%t Table[If[n == 1, 1, Apply[Times, Union[FactorInteger[n][[All, -1]] ]]], {n, 120}] (* _Michael De Vlieger_, Aug 14 2017 *)
%o (PARI) A290107(n) = factorback(vecsort((factor(n)[, 2]), ,8));
%o (Scheme) (define (A290107 n) (A156061 (A181819 n)))
%Y Cf. A156061, A181819.
%Y Differs from A005361 for the first time at n=36.
%Y Differs from A072411 for the first time at n=144, and also from A157754 for the second time (after the initial term).
%K nonn
%O 1,4
%A _Antti Karttunen_, Aug 13 2017
|
