|
%I #11 Sep 08 2022 08:46:23
%S 1,1,3,8,5,27,7,128,27,125,11,10368,13,343,3375,131072,17,118098,19,
%T 2000000,9261,1331,23,6879707136,625,2197,19683,15059072,29,
%U 38443359375,31,2147483648,35937,4913,42875,101559956668416,37,6859,59319,10240000000000,41
%N Numerator of Product_{d|n} (pod(d)/tau(d)) where pod(k) = the product of the divisors of k and tau(k) = the number of the divisors of k.
%C Product_{d|n} (pod(d)/tau(d)) > 1 for all n > 2.
%F a(p) = p for primes p > 2.
%e For n=4; Product_{d|4} (pod(d)/tau(d)) = (pod(1)/tau(1))*(pod(2)/tau(2))*(pod(4)/tau(4)) = (1/1)*(2/2)*(8/3) = 8/3; a(4) = 8.
%p A323760 := proc(n)
%p numer(A266265(n)/A211776(n)) ;
%p end proc:
%p seq(A323760(n),n=1..20) ; # _R. J. Mathar_, Feb 13 2019
%o (Magma) [Numerator(&*[&*[c: c in Divisors(d)] / NumberOfDivisors(d): d in Divisors(n)]): n in [1..100]]
%o (PARI) a(n) = my(p=1, vd); fordiv(n, d, vd = divisors(d); p *= vecprod(vd)/#vd); numerator(p); \\ _Michel Marcus_, Jan 27 2019
%Y Cf. A211776, A266265, A323761 (denominator).
%K nonn,frac
%O 1,3
%A _Jaroslav Krizek_, Jan 26 2019
|
