login
A338697
a(n) = [x^n] Product_{k>=1} 1 / (1 - n^(k-1)*x^k).
4
1, 1, 3, 13, 101, 931, 12391, 178809, 3331721, 66288142, 1589753211, 40104031166, 1183380156013, 36187564837217, 1262524447510383, 45533370885563716, 1834219414937219601, 76016894083755947753, 3479900167920331954531, 162982921698852088968886, 8341707623665223127224821
OFFSET
0,3
LINKS
FORMULA
a(n) = Sum_{k=0..n} p(n,k) * n^(n-k), where p(n,k) is the number of partitions of n into k parts.
a(n) ~ c * n^(n-1), where c = BesselI(1,2) = A096789 = 1.590636854637329... - Vaclav Kotesovec, May 09 2021
MATHEMATICA
Table[SeriesCoefficient[Product[1/(1 - n^(k - 1) x^k), {k, 1, n}], {x, 0, n}], {n, 0, 20}]
Join[{1}, Table[Sum[Length[IntegerPartitions[n, {k}]] n^(n - k), {k, 0, n}], {n, 1, 20}]]
Join[{1}, Table[SeriesCoefficient[x + (n-1)/(n*QPochhammer[1/n, n*x]), {x, 0, n}], {n, 1, 20}]] (* Vaclav Kotesovec, May 09 2021 *)
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Apr 24 2021
STATUS
approved