OFFSET
1,2
COMMENTS
a(n) = a(n-1)*A062758(n).
a(n) = Product_{k=1..n} k^floor(n/k) * (floor(n/k))!.
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..117
MATHEMATICA
f[n_] := Product[ k^DivisorSigma[0, k], {k, n}]; Array[f, 15] (* Robert G. Wilson v, Jun 11 2010 *)
PROG
(Python)
from sympy import divisor_count
from itertools import count, islice
def agen():
an = 1
for k in count(2):
yield an
an *= k**divisor_count(k)
print(list(islice(agen(), 14))) # Michael S. Branicky, May 03 2022
(PARI) a(n) = prod(k=1, n, k^numdiv(k)); \\ Michel Marcus, May 03 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, May 30 2010
EXTENSIONS
a(6) onwards from Robert G. Wilson v and Jon E. Schoenfield, Jun 11 2010
a(14) and beyond from Michael S. Branicky, May 03 2022
STATUS
approved