A366329 a(n) = Product_{k=1..n} Sum_{j=1..k} j^n.

1, 5, 324, 589764, 52393770000, 347773153451938500, 244632735619259069507040000, 24547871392966749661547369532868031040, 455140097017244017295446005144727669016636127744000, 1960564895414510364772369567330640938816177001699555385515625000000

Offset: 1

Author

Vaclav Kotesovec, Oct 07 2023

Keywords

• nonn

Crossrefs

Cf. A031971, A036740, A366305, A366342.

Programs

Mathematica

Table[Product[Sum[j^n, {j, 1, k}], {k, 1, n}], {n, 1, 12}]
Table[Product[HarmonicNumber[k, -n], {k, 1, n}], {n, 1, 12}] // FunctionExpand

Formula

a(n) = A036740(n) * Product_{k=1..n} Sum_{j=1..k} (j/k)^n.

a(n) ~ n!^n * c * d^n, where d = exp(-Integral_{x=0..1} log(1 - exp(-1/x)) dx) = 1.187538543919977798892363400109897833660222697152558038684860736484... and c = exp(1 - 1/(exp(1) - 1)) / (exp(1) - 1) = 0.88399704290317414073109479991305699114875723090346..., updated Apr 19 2024

a(n) ~ c * d^n * (2*Pi)^(n/2) * n^(n*(2*n+1)/2) / exp(n^2 - 1/12).