OFFSET
1,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..300
FORMULA
Completely multiplicative with a(p) = p!.
Sum_{n>=1} 1/a(n) = 1/Product_{p prime} (1 - 1/p!) = 2.42064841562992124691... . - Amiram Eldar, Dec 24 2022
MAPLE
a:= n-> mul(i[1]!^i[2], i=ifactors(n)[2]):
seq(a(n), n=1..30); # Alois P. Heinz, Mar 30 2014
MATHEMATICA
f[p_, e_] := p!^e; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 30] (* Amiram Eldar, Aug 04 2022 *)
PROG
(PARI) a(n)=local(fm=factor(n)); prod(k=1, matsize(fm)[1], fm[k, 1]!^fm[k, 2])
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Franklin T. Adams-Watters, Mar 30 2014
STATUS
approved