A344063 Expansion of Product_{k>=1} (1 + 4^(k-1)*x^k).

1, 1, 4, 20, 80, 384, 1600, 7424, 30720, 143360, 593920, 2703360, 11403264, 51118080, 214958080, 965738496, 4047503360, 17951621120, 76168560640, 334202142720, 1411970498560, 6211596451840, 26203595472896, 114246130073600, 484815908372480, 2101441598586880, 8896148580335616

Offset

0,3

Links

Table of n, a(n) for n=0..26.

Formula

a(n) = Sum_{k=0..A003056(n)} q(n,k) * 4^(n-k), where q(n,k) is the number of partitions of n into k distinct parts.

a(n) ~ (-polylog(2, -1/4))^(1/4) * 4^n * exp(2*sqrt(-polylog(2, -1/4)*n)) / (2*sqrt(5*Pi/4)*n^(3/4)). - Vaclav Kotesovec, May 09 2021

Mathematica

nmax = 26; CoefficientList[Series[Product[(1 + 4^(k - 1) x^k), {k, 1, nmax}], {x, 0, nmax}], x]
Table[Sum[Length[Select[IntegerPartitions[n, {k}], UnsameQ @@ # &]] 4^(n - k), {k, 0, Floor[(Sqrt[8 n + 1] - 1)/2]}], {n, 0, 26}]

Prog

(PARI) seq(n)={Vec(prod(k=1, n, 1 + 4^(k-1)*x^k + O(x*x^n)))} \\ Andrew Howroyd, May 08 2021

Crossrefs

Cf. A003056, A008289, A261568, A304961, A338673, A340103, A344062, A344064, A344065, A344066, A344067, A344068.

Sequence in context: A371408 A292540 A320934 * A055296 A140532 A217482

Adjacent sequences: A344060 A344061 A344062 * A344064 A344065 A344066

Keyword

nonn

Author

Ilya Gutkovskiy, May 08 2021

Status

approved