A320564 - OEIS

The OEIS is supported by the many generous donors to the OEIS Foundation.

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Other ways to Give

A320564

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

1, 1, 3, 10, 30, 87, 249, 705, 1974, 5471, 15032, 40997, 111079, 299151, 801139, 2134251, 5657895, 14930596, 39232009, 102673794, 267692321, 695440442, 1800582809, 4646964755, 11956293758, 30673060344, 78470890246, 200218512582, 509557661691, 1293664233400, 3276659862518

(list; graph; refs; listen; history; text; internal format)

OFFSET

0,3

COMMENTS

First differences of the binomial transform of A026007.

LINKS

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

N. J. A. Sloane, Transforms

FORMULA

G.f.: exp(Sum_{k>=1} (-1)^(k+1)*(1 - x)^k*x^k/(k*((1 - x)^k - x^k)^2)).

a(n) ~ Zeta(3)^(1/6) * 2^(n - 13/12) * exp(3^(4/3) * Zeta(3)^(1/3) * n^(2/3)/4 + (3*Zeta(3))^(2/3) * n^(1/3)/8 - Zeta(3)/16) / (3^(1/3) * sqrt(Pi) * n^(2/3)). - Vaclav Kotesovec, Oct 15 2018

MAPLE

seq(coeff(series(mul((1+x^k/(1-x)^k)^k, k=1..n), x, n+1), x, n), n = 0 .. 30); # Muniru A Asiru, Oct 15 2018

MATHEMATICA

nmax = 30; CoefficientList[Series[Product[(1 + x^k/(1 - x)^k)^k, {k, 1, nmax}], {x, 0, nmax}], x]

nmax = 30; CoefficientList[Series[Exp[Sum[(-1)^(k + 1) (1 - x)^k x^k/(k ((1 - x)^k - x^k)^2), {k, 1, nmax}]], {x, 0, nmax}], x]

CROSSREFS

Cf. A026007, A129519, A294502, A320563.

Sequence in context: A092756 A027205 A026937 * A094306 A257596 A261336

Adjacent sequences: A320561 A320562 A320563 * A320565 A320566 A320567

KEYWORD

nonn

AUTHOR

Ilya Gutkovskiy, Oct 15 2018

STATUS

approved