|
|
A308944
|
|
a(n) = Product_{k=1..n} lcm(n,k) / (k * gcd(n,k)).
|
|
0
|
|
|
1, 1, 3, 4, 125, 9, 16807, 1024, 59049, 15625, 2357947691, 5184, 1792160394037, 282475249, 474609375, 17179869184, 2862423051509815793, 3486784401, 5480386857784802185939, 250000000000, 10382917022245341, 5559917313492231481, 39471584120695485887249589623
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,3
|
|
LINKS
|
|
|
FORMULA
|
a(n) = Product_{d|n} d^(phi(d)-phi(n/d)).
a(n) = n^n / Product_{d|n} d^(2*phi(n/d)).
a(n) = n^(-n) * Product_{d|n} d^(2*phi(d)).
a(n) = n^n / Product_{k=1..n} gcd(n,k)^2.
a(n) = n^(-n) * Product_{k=1..n} lcm(n,k)^2/k^2.
a(p) = p^(p-2), where p is a prime.
|
|
MATHEMATICA
|
Table[Product[LCM[n, k]/(k GCD[n, k]), {k, 1, n}], {n, 1, 23}]
Table[Product[d^(EulerPhi[d] - EulerPhi[n/d]), {d, Divisors[n]}], {n, 1, 23}]
|
|
PROG
|
(PARI) a(n) = prod(k=1, n, lcm(n, k)/(k*gcd(n, k))); \\ Michel Marcus, Jul 02 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|