login
A344098
a(n) = [x^n] Product_{k>=1} (1 + x^k)^binomial(k+n-1,n-1).
3
1, 1, 4, 29, 221, 2027, 21022, 242209, 3060262, 41936745, 618154670, 9735013136, 162892047930, 2882449728121, 53727527279464, 1051276401060921, 21529017626095851, 460231878244308738, 10246160509840187387, 237067632496414877363, 5689786581042000827057, 141415234722601777758232
OFFSET
0,3
MATHEMATICA
Table[SeriesCoefficient[Product[(1 + x^k)^Binomial[k + n - 1, n - 1], {k, 1, n}], {x, 0, n}], {n, 0, 21}]
A[n_, k_] := A[n, k] = If[n == 0, 1, (1/n) Sum[Sum[(-1)^(j/d + 1) d Binomial[d + k - 1, k - 1], {d, Divisors[j]}] A[n - j, k], {j, 1, n}]]; a[n_] := A[n, n]; Table[a[n], {n, 0, 21}]
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, May 09 2021
STATUS
approved