OFFSET
0,3
COMMENTS
Subset of A007955.
EXAMPLE
a(0) = 1 because 0^0 = 1, whose only divisor is 1, so the product of divisors is 1.
a(1) = 1 because 1^1 = 1, so the product of divisors is 1.
a(3) = 729 because 3^3 = 27, whose divisors are (1, 3, 9, 27), and their product is 729.
PROG
(SageMath)
[ product((1*i^i).divisors()) for i in range(10) ]
(Magma) [&*Divisors(n^n): n in [0..8]]; // Marius A. Burtea, Jul 26 2019
(Python)
from math import isqrt, prod
from sympy import factorint
def A309377(n): return (isqrt(n**n) if (c:=prod(n*e+1 for e in factorint(n).values())) & 1 else 1)*n**(n*(c//2)) # Chai Wah Wu, Jun 25 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Hauke Löffler, Jul 26 2019
STATUS
approved