|
%I #55 Sep 08 2022 08:46:20
%S 1,2,3,8,5,1,7,64,27,25,11,3,13,49,25,1024,17,27,19,1000,441,121,23,
%T 16,125,169,729,7,29,625,31,32768,121,289,1225,5832,37,361,1521,6400,
%U 41,2401,43,1331,3375,529,47,3072,343,31250,289,17576,53,6561,3025,49
%N a(n) = denominator of Product_{d|n} (sigma(d)/d) where sigma(k) = the sum of the divisors of k (A000203).
%F a(n) = 1 for n = 1, 6, ... (no other n <= 5*10^6).
%F a(n) = n for n = primes (A000040).
%F a(n) = n * sqrt(n) for n in A280076 (union of 1 and squares of primes (A001248)).
%e For n=4; Product_{d|4} (sigma(d)/d) = (sigma(1)/1)*(sigma(2)/2)*(sigma(4)/4) = (1/1)*(3/2)*(7/4) = 21/8; a(4) = 8.
%t Array[Denominator@ Product[DivisorSigma[1, d]/d, {d, Divisors@ #}] &, 51]
%o (Magma) [Denominator(&*[&+[c: c in Divisors(d)] / d: d in Divisors(n)]): n in [1..100]]
%o (PARI) a(n) = my(p=1); fordiv(n, d, p *= sigma(d)/d); denominator(p); \\ _Michel Marcus_, Jan 21 2019
%Y Cf. A000203, A318491, A318492, A322673 (numerator).
%K nonn,frac
%O 1,2
%A _Jaroslav Krizek_, Jan 21 2019
|
