%I #6 Jul 19 2024 14:28:11
%S 1,4,6,12,10,24,14,32,27,8,22,72,26,56,60,80,34,54,38,120,4,88,46,192,
%T 75,104,54,56,58,48,62,64,132,136,28,162,74,152,52,320,82,16,86,264,
%U 135,184,94,160,49,30,204,104,106,216,20,64,76,232,118,720,122,248
%N Denominator of the mean abundancy index of the divisors of n.
%H Amiram Eldar, <a href="/A374778/b374778.txt">Table of n, a(n) for n = 1..10000</a>
%e For n = 2, n has 2 divisors, 1 and 2. Their abundancy indices are sigma(1)/1 = 1 and sigma(2)/2 = 3/2, and their mean abundancy index is (1 + 3/2)/2 = 5/4. Therefore a(2) = denominator(5/4) = 4.
%t f[p_, e_] := ((e+1)*p^2 - (e+2)*p + p^(-e))/((e+1)*(p-1)^2); a[1] = 1; a[n_] := Denominator[Times @@ f @@@ FactorInteger[n]]; Array[a, 100]
%o (PARI) a(n) = {my(f = factor(n), p, e); denominator(prod(i = 1, #f~, p = f[i, 1]; e = f[i, 2]; ((e+1)*p^2 - (e+2)*p + p^(-e))/((e+1)*(p-1)^2)));}
%Y Cf. A374777 (numerators).
%K nonn,easy,frac
%O 1,2
%A _Amiram Eldar_, Jul 19 2024