OFFSET
1,2
COMMENTS
From G. C. Greubel, Dec 31 2015: (Start)
for n>=1: 10^3|a(10*n), 10^6|a(20*n), 10^9|a(30*n).
for n>=0: 10^6|a(60*n+50), 10^9|a(60*n+70). (End)
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..1200
FORMULA
a(p) = p for p = prime.
a(n) = Product_{d|n} d^tau(n/d). - Ridouane Oudra, Apr 09 2023
EXAMPLE
For n = 6; with b(n) = A007955(n); a(6) = b(1)*b(2)*b(3)*b(6) = 1*2*3*36 = 216.
MAPLE
A266265 := proc(n)
mul( A007955(d), d=numtheory[divisors](n)) ;
end proc:
seq(A266265(n), n=1..10) ; # R. J. Mathar, Feb 13 2019
MATHEMATICA
Table[Product[Times @@ Divisors@ d, {d, Divisors@ n}], {n, 38}] (* Michael De Vlieger, Dec 31 2015 *)
PROG
(Magma) [&*[&*[b: b in Divisors(d)]: d in Divisors(n)]: n in [1..100]]
(PARI) a(n) = {my(d = divisors(n)); prod(k=1, #d, dd = divisors(d[k]); prod(kk=1, #dd, dd[kk])); } \\ Michel Marcus, Dec 27 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Dec 25 2015
STATUS
approved