cp's OEIS Frontend

A324981 a(n) = Product_{d|n} (d*pod(d)) where pod(k) = the product of the divisors of k (A007955).

%I A324981 #12 Sep 08 2022 08:46:24
%S A324981 1,4,9,128,25,7776,49,65536,2187,100000,121,5159780352,169,537824,
%T A324981 759375,1073741824,289,198359290368,361,512000000000,4084101,5153632,
%U A324981 529,21035720123168587776,78125,11881376,43046721,10578455953408,841,15943230000000000000,961
%N A324981 a(n) = Product_{d|n} (d*pod(d)) where pod(k) = the product of the divisors of k (A007955).
%F A324981 a(n) = (Product_{d|n} d) * (Product_{d|n} pod(d)) = A007955(n) * A266265(n).
%F A324981 a(p) = p^2 for p = primes (A000040).
%e A324981 a(6) = 1*pod(1) * 2*pod(2) * 3*pod(3) * 6*pod(6) = (1*1) * (2*2) * (3*3) * (6*36) = 7776.
%t A324981 Array[Times @@ Map[# Apply[Times, Divisors@ #] &, Divisors@ #] &, 31] (* _Michael De Vlieger_, Mar 24 2019 *)
%o A324981 (Magma) [&*[d * &*[c: c in Divisors(d)]: d in Divisors(n)]: n in [1..100]]
%o A324981 (PARI) a(n) = my(d=divisors(n), p=1); fordiv(n, d, p*=d*vecprod(divisors(d))); p;  \\ _Michel Marcus_, Mar 22 2019
%Y A324981 Cf. A007955, A174933 (Sum_{d|n} (d*pod(d))), A266265.
%K A324981 nonn
%O A324981 1,2
%A A324981 _Jaroslav Krizek_, Mar 22 2019