A156616 G.f.: Product_{n>0} ((1+x^n)/(1-x^n))^n.

1, 2, 6, 16, 38, 88, 196, 420, 878, 1794, 3584, 7032, 13572, 25792, 48352, 89512, 163774, 296444, 531234, 943072, 1659560, 2896376, 5015700, 8622108, 14718652, 24960138, 42062200, 70458160, 117349856, 194381704, 320295312, 525123604

Offset

0,2

Comments

Generating function for a sum over strict plane partitions weighted with 2 powered to their number of connected components.

The inverse Euler transform is apparently 2, 3, 6, 6, 10, 9, 14, 12, 18, 15, 22, 18, 26, 21, ..., A016825 interlaced with A008585. - R. J. Mathar, Apr 23 2009

In general, for m >= 1, if g.f. = Product_{k>=1} ((1+x^k)/(1-x^k))^(m*k), then a(n) ~ exp(m/12 + 3/2 * (7*m*Zeta(3)/2)^(1/3) * n^(2/3)) * m^(1/6 + m/36) * (7*Zeta(3))^(1/6 + m/36) / (A^m * 2^(2/3 + m/9) * sqrt(3*Pi) * n^(2/3 + m/36)), where Zeta(3) = A002117 and A = A074962 is the Glaisher-Kinkelin constant. - Vaclav Kotesovec, Aug 17 2015

In general, for m >= 0, if g.f. = Product_{k>=1} ((1+x^k)/(1-x^k))^(k^m), then a(n) ~ ((2^(m+2)-1) * Gamma(m+2) * Zeta(m+2) / (2^(2*m+3) * n))^((1-2*Zeta(-m))/(2*m+4)) * exp((m+2)/(m+1) * ((2^(m+2)-1) * n^(m+1) * Gamma(m+2) * Zeta(m+2) / 2^(m+1))^(1/(m+2)) + Zeta'(-m)) / sqrt((m+2)*Pi*n). - Vaclav Kotesovec, Aug 19 2015

Links

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from Vaclav Kotesovec)

Ali H. Al-Saedi, Congruences for restricted plane overpartitions modulo 4 and 8, Raman. J. 48 (2) (2019) 251

Vaclav Kotesovec, A method of finding the asymptotics of q-series based on the convolution of generating functions, arXiv:1509.08708 [math.CO], Sep 30 2015, p. 18.

Mirjana Vuletic, A generalization of MacMahon's formula, Trans. Am. Math. Soc. 361 (2009) 2789-2804.

Formula

Convolve A000219 with A026007.

O.g.f.: exp( Sum_{n>=1} (sigma_2(2n) - sigma_2(n))/2 *x^n/n ), where sigma_2(n) is the sum of squares of divisors of n (A001157). - Paul D. Hanna, May 01 2010

a(n) ~ exp(1/12 + 3 * 2^(-4/3) * (7*Zeta(3))^(1/3) * n^(2/3)) * (7*Zeta(3))^(7/36) / (A * 2^(7/9) * sqrt(3*Pi) * n^(25/36)), where Zeta(3) = A002117 and A = A074962 is the Glaisher-Kinkelin constant. - Vaclav Kotesovec, Aug 17 2015

a(0) = 1, a(n) = (2/n)*Sum_{k=1..n} A076577(k)*a(n-k) for n > 0. - Seiichi Manyama, Apr 30 2017

G.f.: A(x) = exp( 2*Sum_{n >= 0} x^(2*n+1)/((2*n+1)*(1 - x^(2*n+1))^2) ). Cf. A000122 and A302237. - Peter Bala, Dec 23 2021

Mathematica

nmax = 40; CoefficientList[Series[Product[((1+x^k)/(1-x^k))^k, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 17 2015 *)

Prog

(PARI) {a(n)=polcoeff(exp(sum(m=1, n, (sigma(2*m, 2)-sigma(m, 2))/2*x^m/m)+x*O(x^n)), n)} \\ Paul D. Hanna, May 01 2010

Crossrefs

Cf. A000219, A015128, A026007, A261384, A261386, A261389.

Cf. A206622, A206623, A206624, A261519, A261520.

Cf. A285462, A285447, A285460, A285461.

Cf. A285446, A285458, A285459.

Cf. A306081.

Sequence in context: A097813 A167821 A093041 * A330886 A265758 A265107

Adjacent sequences: A156613 A156614 A156615 * A156617 A156618 A156619

Keyword

easy,nonn

Author

R. J. Mathar, Feb 11 2009

Status

approved