A156616 -id:A156616 - cp's OEIS frontend

A015128 Number of overpartitions of n: an overpartition of n is an ordered sequence of nonincreasing integers that sum to n, where the first occurrence of each integer may be overlined.

1, 2, 4, 8, 14, 24, 40, 64, 100, 154, 232, 344, 504, 728, 1040, 1472, 2062, 2864, 3948, 5400, 7336, 9904, 13288, 17728, 23528, 31066, 40824, 53408, 69568, 90248, 116624, 150144, 192612, 246256, 313808, 398640, 504886, 637592, 802936, 1008448

Offset: 0

N. J. A. Sloane

The over-partition function.
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Also the number of jagged partitions of n.
According to Ramanujan (1913) a(n) is close to (cosh(x)-sinh(x)/x)/(4*n) where x=Pi*sqrt(n). - Michael Somos, Mar 17 2003
Number of partitions of 2n with all odd parts occurring with even multiplicities. There is no restriction on the even parts. Cf. A006950, A046682. - Mamuka Jibladze, Sep 05 2003
Number of partitions of n where there are two kinds of odd parts. - Joerg Arndt, Jul 30 2011. Or, in Gosper's words, partitions into red integers and blue odd integers. - N. J. A. Sloane, Jul 04 2016.
Coincides with the sequence of numbers of nilpotent conjugacy classes in the Lie algebras sp(n), n=0,1,2,3,... (the case n=0 being degenerate). A006950, this sequence and A000041 together cover the nilpotent conjugacy classes in the classical A,B,C,D series of Lie algebras. - Alexander Elashvili, Sep 08 2003
Also, number of 01-partitions of n. A 01-partition of n is a weakly decreasing sequence of m nonnegative integers n(i) such that sum(n(i))=n, n(m)>0, n(j)>=n(j+1)-1 and n(j)>=n(j+2). They are special cases of jagged partitions.
a(8n+7) is divisible by 64 (from Fortin/Jacob/Mathieu paper).
Smallest sequence of even numbers (except a(0)) which is the Euler transform of a sequence of positive integers. - Franklin T. Adams-Watters, Oct 16 2006
Convolution of A000041 and A000009. - Vladeta Jovovic, Nov 26 2002
Equals A022567 convolved with A035363. - Gary W. Adamson, Jun 09 2009
Equals the infinite product [1,2,2,2,...] * [1,0,2,0,2,0,2,...] * [1,0,0,2,0,0,2,0,0,2,...] * ... . - Gary W. Adamson, Jul 05 2009
Equals A182818 convolved with A010815. - Gary W. Adamson, Jul 20 2012
Partial sums of A211971. - Omar E. Pol, Jan 09 2014
Also 1 together with the row sums of A235790. - Omar E. Pol, Jan 19 2014
Antidiagonal sums of A284592. - Peter Bala, Mar 30 2017
The overlining method is equivalent to enumerating the k-subsets of the distinct parts of the i-th partition. - Richard Joseph Boland, Sep 02 2021

			G.f. = 1 + 2*q + 4*q^2 + 8*q^3 + 14*q^4 + 24*q^5 + 40*q^6 + 64*q^7 + 100*q^8 + ...
For n = 4 the 14 overpartitions of 4 are [4], [4'], [2, 2], [2', 2], [3, 1], [3', 1], [3, 1'], [3', 1'], [2, 1, 1], [2', 1, 1], [2, 1', 1], [2', 1', 1], [1, 1, 1, 1], [1', 1, 1, 1]. - _Omar E. Pol_, Jan 19 2014

J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 103.
R. W. Gosper, Experiments and discoveries in q-trigonometry, in Symbolic Computation, Number Theory, Special Functions, Physics and Combinatorics. Editors: F. G. Garvan and M. E. H. Ismail. Kluwer, Dordrecht, Netherlands, 2001, pp. 79-105. See the function g(q).
James R. Newman, The World of Mathematics, Simon and Schuster, 1956, Vol. I p. 372.

Vaclav Kotesovec, Table of n, a(n) for n = 0..10000 (first 1000 terms from T. D. Noe)
Gert Almkvist, Asymptotic formulas and generalized Dedekind sums, Exper. Math., 7 (No. 4, 1998), pp. 343-359.
Brennan Benfield and Arindam Roy, Log-concavity And The Multiplicative Properties of Restricted Partition Functions, arXiv:2404.03153 [math.NT], 2024.
Noureddine Chair, Partition identities from Partial Supersymmetry, arXiv:hep-th/0409011, 2004.
Shi-Chao Chen, On the number of overpartitions into odd parts, Discrete Math. 325 (2014), 32--37. MR3181230.
William Y.C. Chen and Ernest X.W. Xia, Proof of a conjecture of Hirschhorn and Sellers on overpartitions, arXiv:1307.4155 [math.CO], 2013; Acta Arith. 163 (2014), no. 1, 59--69. MR3194057
Sylvie Corteel, Particle seas and basic hypergeometric series, Advances Appl. Math., 31 (2003), 199-214.
Sylvie Corteel and Jeremy Lovejoy, Frobenius partitions and the combinatorics of Ramanujan's 1 psi 1 summation, J. Combin. Theory A 97 (2002), 177-183.
Sylvie Corteel and Jeremy Lovejoy, Overpartitions, Trans. Amer. Math. Soc., 356 (2004), 1623-1635.
Brian Drake, Limits of areas under lattice paths, Discrete Math. 309 (2009), no. 12, 3936-3953.
Benjamin Engel, Log-concavity of the overpartition function, The Ramanujan Journal, Vol. 43, No. 2 (2017), pp. 229-241; arXiv preprint, arXiv:1412.4603 [math.NT], 2014.
Alex Fink, Richard K. Guy and Mark Krusemeyer, Partitions with parts occurring at most thrice, Contributions to Discrete Mathematics, Vol 3, No 2 (2008).
J.-F. Fortin, P. Jacob and P. Mathieu, Jagged partitions, The Ramanujan Journal, Vol. 10, No. 2 (2005), pp. 215-235; arXiv preprint, arXiv:math/0310079 [math.CO], 2003-2005.
Frank Garvan, Table of a(n) for n = 1..10000.
R. W. Gosper, Experiments and discoveries in q-trigonometry, in F. G. Garvan and M. E. H. Ismail (eds.), Symbolic computation, number theory, special functions, physics and combinatorics, Springer, Boston, MA, 2001, pp. 79-105; preprint.
R. W. Gosper, q-Trigonometry: Some Prefatory Afterthoughts
William J. Keith, Restricted k-color partitions, The Ramanujan Journal, Vol. 40, No. 1 (2016), pp. 71-92; arXiv preprint, arXiv:1408.4089 [math.CO], 2014.
Byungchan Kim, A short note on the overpartition function, Discr. Math., 309 (2009), 2528-2532.
Byungchan Kim, Overpartition pairs modulo powers of 2, Discrete Math., 311 (2011), 835-840.
Vaclav Kotesovec, A method of finding the asymptotics of q-series based on the convolution of generating functions, arXiv:1509.08708 [math.CO], 2015-2016.
Jeremy Lovejoy, Gordon's theorem for overpartitions, J. Combin. Theory A 103 (2003), 393-401.
Karl Mahlburg, The overpartition function modulo small powers of 2, Discr. Math., 286 (2004), 263-267.
Mircea Merca, Overpartitions in terms of 2-adic valuation, Aequat. Math. (2024). See p. 9.
Igor Pak, Partition bijections, a survey, The Ramanujan Journal, Vol. 12, No. 1 (2006), pp. 5-75; alternative link.
Maxie D. Schmidt, Exact Formulas for the Generalized Sum-of-Divisors Functions, arXiv:1705.03488 [math.NT], 2017-2019. See Example 4.1, p. 13.
Michael Somos, Introduction to Ramanujan theta functions
Liuquan Wang, Another Proof of a Conjecture by Hirschhorn and Sellers on Overpartitions, J. Int. Seq. 17 (2014) # 14.9.8.
Eric Weisstein's World of Mathematics, Modular Equation.
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions.
Michael P. Zaletel and Roger S. K. Mong, Exact matrix product states for quantum Hall wave functions, Physical Review B, Vol. 86, No. 24 (2012), 245305; arXiv preprint, arXiv:1208.4862 [cond-mat.str-el], 2012. - From _N. J. A. Sloane_, Dec 25 2012

Cf. A022567, A035363, A002448, A106507, A156616, A261386, A265835, A014968, A284592, A303662.
See A004402 for a version with signs.
Column k=2 of A321884.
Cf. A002513.

# JacobiTheta4 is defined in A002448.
A015128List(len) = JacobiTheta4(len, -1)
A015128List(40) |> println # Peter Luschny, Mar 12 2018

mul((1+x^n)/(1-x^n),n=1..256): seq(coeff(series(%,x,n+1),x,n), n=0..40);
# second Maple program:
b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
      b(n, i-1) +2*add(b(n-i*j, i-1), j=1..n/i)))
    end:
a:= n-> b(n$2):
seq(a(n), n=0..40);  # Alois P. Heinz, Feb 10 2014
a_list := proc(len) series(1/JacobiTheta4(0,x),x,len+1); seq(coeff(%,x,j),j=0..len) end: a_list(39); # Peter Luschny, Mar 14 2017

max = 39; f[x_] := Exp[Sum[(DivisorSigma[1, 2*n] - DivisorSigma[1, n])*(x^n/n), {n, 1, max}]]; CoefficientList[ Series[f[x], {x, 0, max}], x] (* Jean-François Alcover, Jun 11 2012, after Joerg Arndt *)
a[ n_] := SeriesCoefficient[ QHypergeometricPFQ[ {-1}, {}, x, x], {x, 0, n}]; (* Michael Somos, Mar 11 2014 *)
QP = QPochhammer; s = QP[q^2]/QP[q]^2 + O[q]^40; CoefficientList[s + O[q]^100, q] (* Jean-François Alcover, Nov 25 2015, after Michael Somos *)
Table[Sum[PartitionsP[n-k]*PartitionsQ[k], {k, 0, n}], {n, 0, 50}] (* Vaclav Kotesovec, Nov 28 2015 *)
(QPochhammer[-x, x]/QPochhammer[x, x] + O[x]^50)[[3]] (* Vladimir Reshetnikov, Nov 12 2016 *)
nmax = 100; p = ConstantArray[0, nmax+1]; p[[1]] = 1; Do[p[[n+1]] = 0; k = 1; While[n + 1 - k^2 > 0, p[[n+1]] += (-1)^(k+1)*p[[n + 1 - k^2]]; k++;]; p[[n+1]] = 2*p[[n+1]];, {n, 1, nmax}]; p (* Vaclav Kotesovec, Apr 11 2017 *)
a[ n_] := SeriesCoefficient[ 1 / EllipticTheta[ 4, 0, x], {x, 0, n}]; (* Michael Somos, Nov 15 2018 *)
a[n_] := Sum[2^Length[Union[IntegerPartitions[n][[i]]]], {i, 1, PartitionsP[n]}]; (* Richard Joseph Boland, Sep 02 2021 *)
n = 39; CoefficientList[Product[(1 + x^k)/(1 - x^k), {k, 1, n}] + O[x]^(n + 1), x] (* Oliver Seipel, Sep 19 2021 *)

{a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A) / eta(x + A)^2, n))}; /* Michael Somos, Nov 01 2008 */

{a(n)=polcoeff(exp(sum(m=1,n\2+1,2*x^(2*m-1)/(1-x^(2*m-1)+x*O(x^n))/(2*m-1))),n)} /* Paul D. Hanna, Aug 06 2009 */

N=66; x='x+O('x^N); gf=exp(sum(n=1,N,(sigma(2*n)-sigma(n))*x^n/n));Vec(gf) /* Joerg Arndt, Jul 30 2011 */

lista(nn) = {q='q+O('q^nn); Vec(eta(q^2)/eta(q)^2)} \\ Altug Alkan, Mar 20 2018

# uses[EulerTransform from A166861]
a = BinaryRecurrenceSequence(0, 1, 1, 2)
b = EulerTransform(a)
print([b(n) for n in range(40)]) # Peter Luschny, Nov 11 2020

Euler transform of period 2 sequence [2, 1, ...]. - Michael Somos, Mar 17 2003
G.f.: Product_{m>=1} (1 + q^m)/(1 - q^m).
G.f.: 1 / (Sum_{m=-inf..inf} (-q)^(m^2)) = 1/theta_4(q).
G.f.: 1 / Product_{m>=1} (1 - q^(2*m)) * (1 - q^(2*m-1))^2.
G.f.: exp( Sum_{n>=1} 2*x^(2*n-1)/(1 - x^(2*n-1))/(2*n-1) ). - Paul D. Hanna, Aug 06 2009
G.f.: exp( Sum_{n>=1} (sigma(2*n) - sigma(n))*x^n/n ). - Joerg Arndt, Jul 30 2011
G.f.: Product_{n>=0} theta_3(q^(2^n))^(2^n). - Joerg Arndt, Aug 03 2011
A004402(n) = (-1)^n * a(n). - Michael Somos, Mar 17 2003
Expansion of eta(q^2) / eta(q)^2 in powers of q. - Michael Somos, Nov 01 2008
Expansion of 1 / phi(-q) in powers of q where phi() is a Ramanujan theta function. - Michael Somos, Nov 01 2008
Convolution inverse of A002448. - Michael Somos, Nov 01 2008
Recurrence: a(n) = 2*Sum_{m>=1} (-1)^(m+1) * a(n-m^2).
a(n) = (1/n)*Sum_{k=1..n} (sigma(2*k) - sigma(k))*a(n-k). - Vladeta Jovovic, Dec 05 2004
G.f.: Product_{i>=1} (1 + x^i)^A001511(2i) (see A000041). - Jon Perry, Jun 06 2004
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = w^4 * (u^4 + v^4) - 2 * u^2 * v^6. - Michael Somos, Nov 01 2008
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^3), A(x^6)) where f(u1, u2, u3, u6) = u6^3 * (u1^2 + u3^2) - 2 * u1 * u2 * u3^3. - Michael Somos, Nov 01 2008
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^3), A(x^6)) where f(u1, u2, u3, u6) = u2^3 * (u3^2 - 3 * u1^2) + 2 * u1^3 * u3 * u6. - Michael Somos, Nov 01 2008
G.f. is a period 1 Fourier series which satisfies f(-1 / (16 t)) = 32^(-1/2) (t/i)^(-1/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A106507. - Michael Somos, Nov 01 2008
a(n) = 2*A014968(n), n >= 1. - Omar E. Pol, Jan 19 2014
a(n) ~ Pi * BesselI(3/2, Pi*sqrt(n)) / (4*sqrt(2)*n^(3/4)). - Vaclav Kotesovec, Jan 11 2017
Let T(n,k) = the number of partitions of n with parts 1 through k of two kinds, T(n,0) = A000041(n), the number of partitions of n. Then a(n) = T(n,0) + T(n-1,1) + T(n-3,2) + T(n-6,3) + T(n-10,4) + T(n-15,5) + ... .   Gregory L. Simay, May 29 2019
For n >= 1, a(n) = Sum_{k>=1} 2^k * A116608(n,k). - Gregory L. Simay, Jun 01 2019
Sum_{n>=1} 1/a(n) = A303662. - Amiram Eldar, Nov 15 2020
a(n) = Sum_{i=1..p(n)} 2^(d(n,i)), where d(n,i) is the number of distinct parts in the i-th partition of n. - Richard Joseph Boland, Sep 02 2021
G.f.: A(x) = exp( Sum_{n >= 1} x^n*(2 + x^n)/(n*(1 - x^(2*n))) ). - Peter Bala, Dec 23 2021
G.f. A(q) satisfies (3*A(q)/A(q^9) - 1)^3 = 9*A(q)^4/A(q^3)^4 - 1. - Paul D. Hanna, Oct 14 2024

Minor edits by Vaclav Kotesovec, Sep 13 2014

A073592 Euler transform of negative integers.

Original entry on oeis.org

1, -1, -2, -1, 0, 4, 4, 7, 3, -2, -9, -17, -25, -24, -13, -1, 32, 61, 97, 111, 112, 74, 8, -108, -243, -392, -512, -569, -542, -358, -33, 473, 1078, 1788, 2395, 2865, 2955, 2569, 1496, -245, -2751, -5783, -9121, -12299, -14739, -15806, -14719, -10930, -3813, 6593, 20284, 36139, 53081, 68620, 80539

Offset: 0

Vladeta Jovovic, Aug 28 2002

sign

1/A(x) is g.f. for A000219.

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from Vaclav Kotesovec)
E. M. Wright, Coefficients of a reciprocal generating function, Quart. J. Math. 17 (1) (1966) 39-43, ADS Abstracts.
N. J. A. Sloane, Transforms

Column k=1 of A283272.

Cf. A000219, A001157, A001478, A026007, A156616, A255528.

a:= proc(n) option remember; `if`(n=0, 1, -add(
      numtheory[sigma][2](j)*a(n-j), j=1..n)/n)
    end:
seq(a(n), n=0..60);  # Alois P. Heinz, Mar 12 2015

nmax=50; CoefficientList[Series[Exp[Sum[-x^k/(k*(1-x^k)^2),{k,1,nmax}]],{x,0,nmax}],x] (* Vaclav Kotesovec, Mar 02 2015 *)
a[n_]:= a[n] = -1/n*Sum[DivisorSigma[2,k]*a[n-k],{k,1,n}]; a[0]=1; Table[a[n],{n,0,100}] (* Vaclav Kotesovec, Mar 02 2015 *)

# uses[EulerTransform from A166861]
b = EulerTransform(lambda n: -n)
print([b(n) for n in range(55)]) # Peter Luschny, Nov 11 2020

G.f.: Product_{k>0} (1-x^k)^k.
a(n) = -1/n*Sum_{k=1..n} sigma[2](k)*a(n-k).
G.f.: exp( Sum_{n>=1} -sigma_2(n)*x^n/n ). - Seiichi Manyama, Mar 04 2017

A270919 Coefficient of x^n in Product_{k>=1} ((1 + x^k) / (1 - x^k))^n.

Original entry on oeis.org

1, 2, 12, 80, 552, 3912, 28224, 206208, 1520784, 11297546, 84413912, 633713808, 4776117216, 36115518376, 273868321536, 2081866609920, 15859616674336, 121046064563376, 925411686479820, 7085465166635440, 54323193841192752, 416993869451825424, 3204447137019290944

Offset: 0

Vaclav Kotesovec, Mar 25 2016

nonn

From Peter Bala, Apr 18 2023: (Start)
The Gauss congruences a(n*p^k) == a(n*p^(k-1)) (mod p^k) hold for all primes p and all positive integers n and k.
Conjecture: the supercongruence a(p) == 2*p + 2 (mod p^2) holds for all primes p. Cf. A291697. (End)

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

Main diagonal of A288515.

Cf. A008485, A008705, A156616, A270913, A270923, A270924, A291697, A309986, A362408.

Table[SeriesCoefficient[Product[((1+x^k)/(1-x^k))^n, {k, 1, n}], {x, 0, n}], {n, 0, 25}]
Table[SeriesCoefficient[(QPochhammer[-1, x]/QPochhammer[x, x])^n, {x, 0, n}]/2^n, {n, 0, 25}]
(* Calculation of constants {d,c}: *) eq = FindRoot[{2*s*QPochhammer[r*s] == QPochhammer[-1, r*s], (Log[1 - r*s] + QPolyGamma[0, 1, r*s])/Log[r*s] + r*((Derivative[0, 1][QPochhammer][-1, r*s] - 2*s*Derivative[0, 1][QPochhammer][r*s, r*s]) / (2*QPochhammer[r*s])) == 1}, {r, 1/8}, {s, 2}, WorkingPrecision -> 1000]; {N[1/r /. eq, 120], val = Sqrt[(1 - r*s)*Log[r*s]^2*(QPochhammer[r*s] / (Pi*(-r*s*(-1 + r*s) * Log[r*s]*(4*(2*ArcTanh[1 - 2*r*s] + QPolyGamma[0, 1, r*s])* Derivative[0, 1][QPochhammer][r*s, r*s] + r*Log[r*s]*(Derivative[0, 2][QPochhammer][-1, r*s] - 2*s*Derivative[0, 2][QPochhammer][r*s, r*s])) + 2*QPochhammer[r*s] * (4*r*s*ArcTanh[1 - 2*r*s] + 2*(-1 + (-1 + r*s)*ArcTanh[1 - 2*r*s])*Log[1 - r*s] - (-1 + r*s)*(-2 + Log[r*s] - 2*Log[1 - r*s])*QPolyGamma[0, 1, r*s] + (-1 + r*s) * QPolyGamma[0, 1, r*s]^2 + (-1 + r*s)*(QPolyGamma[1, 1, r*s] - 2*r*s*Log[r*s]* Derivative[0, 0, 1][QPolyGamma][0, 1, r*s])))))] /. eq; N[Chop[val], -Floor[Log[10, Abs[Im[val]]]] - 3]} (* Vaclav Kotesovec, Oct 03 2023 *)

a(n) ~ c * d^n / sqrt(n), where d = 7.862983395705905261519347909953827161057584... and c = 0.299856802806668079413694689903953367699319...

a(n) = [x^n] 1/theta_4(x)^n, where theta_4() is the Jacobi theta function. - Ilya Gutkovskiy, Nov 03 2017

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

Original entry on oeis.org

1, 2, 10, 36, 118, 376, 1148, 3376, 9654, 26894, 73192, 195188, 510948, 1315048, 3332720, 8326448, 20529526, 49998884, 120379574, 286726340, 676057144, 1578880480, 3654180236, 8385122192, 19085029540, 43103203626, 96630606968, 215105226728, 475608824400

Offset: 0

Paul D. Hanna, Feb 10 2012

nonn

Compare g.f. to: Product_{n>0} (1+x^n)/(1-x^n) = exp( Sum_{n>=1} (sigma(2*n) - sigma(n))*x^n/n ) which equals 1/theta_4(x) = 1/(1 + 2*Sum_{n>=1} (-x)^(n^2)).
Convolution of A023871 and A027998. - Vaclav Kotesovec, Aug 19 2015
In general, if g.f. = Product_{k>=1} ((1 + x^k)/(1 - x^k))^(c2*k^2 + c1*k + c0) and c2>0, then a(n) ~ exp(Pi * 2^(5/4) * c2^(1/4) * n^(3/4) / 3 + 7*c1 * Zeta(3) * sqrt(n) / (Pi^2 * sqrt(2*c2)) + (c0*Pi / (2^(5/4) * c2^(1/4)) - 49*c1^2 * Zeta(3)^2 / (2^(5/4) * c2^(5/4) * Pi^5)) * n^(1/4) + 22411 * c1^3 * Zeta(3)^3 / (196 * c2^2 * Pi^8) - 7*c0*c1 * Zeta(3) / (4*c2 * Pi^2) - c2 * Zeta(3) / (4*Pi^2) + c1/12) * Pi^(c1/12) * c2^(1/8 + c0/8 + c1/48) / (A^c1 * 2^(15/8 + 11*c0/8 + 7*c1/48) * n^(5/8 + c0/8 + c1/48)), where A is the Glaisher-Kinkelin constant A074962. - Vaclav Kotesovec, Nov 08 2017
Let A(x) denote the g.f. and let m be an integer. Define a sequence by u(n) = [x^n] A(x)^(m*n). We conjecture that the supercongruence u(n*p^r) == u(n*p^(r-1)) (mod p^(3*r)) holds for all positive integers n and r and all primes p >= 5. Cf. A380582. - Peter Bala, Jan 21 2025

			G.f.: A(x) = 1 + 2*x + 10*x^2 + 36*x^3 + 118*x^4 + 376*x^5 + 1148*x^6 +...
where A(x) = (1+x)/(1-x) * (1+x^2)^4/(1-x^2)^4 * (1+x^3)^9/(1-x^3)^9 *...
Also, A(x) = Euler transform of [2,7,18,28,50,63,98,112,162,175,...]:
A(x) = 1/((1-x)^2*(1-x^2)^7*(1-x^3)^18*(1-x^4)^28*(1-x^5)^50*(1-x^6)^63*...).

Seiichi Manyama, Table of n, a(n) for n = 0..9042 (terms 0..1000 from Vaclav Kotesovec)
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. 23.

Cf. A156616, A206623, A206624, A001158 (sigma_3), A380582.

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

{a(n)=polcoeff(prod(m=1,n+1,((1+x^m)/(1-x^m+x*O(x^n)))^(m^2)),n)}

{a(n)=polcoeff(exp(sum(m=1, n, (sigma(2*m, 3)-sigma(m, 3))/4*x^m/m)+x*O(x^n)), n)}

{a(n)=local(InvEulerGF=x*(2+7*x+12*x^2+7*x^3+2*x^4)/(1-x^2+x*O(x^n))^3);polcoeff(1/prod(k=1,n,(1-x^k+x*O(x^n))^polcoeff(InvEulerGF,k)),n)}
for(n=0,35,print1(a(n),", "))

G.f.: exp( Sum_{n>=1} (sigma_3(2*n) - sigma_3(n))/4 * x^n/n ), where sigma_3(n) is the sum of cubes of divisors of n (A001158).
The inverse Euler transform has g.f.: x*(2 + 7*x + 12*x^2 + 7*x^3 + 2*x^4)/(1-x^2)^3.
a(n) ~ exp(2^(5/4)*Pi*n^(3/4)/3 - Zeta(3)/(4*Pi^2)) / (2^(15/8) * n^(5/8)), where Zeta(3) = A002117. - Vaclav Kotesovec, Aug 19 2015
a(0) = 1, a(n) = (2/n)*Sum_{k=1..n} A007331(k)*a(n-k) for n > 0. - Seiichi Manyama, Apr 30 2017

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

Original entry on oeis.org

1, 4, 16, 56, 176, 520, 1456, 3896, 10048, 25100, 60960, 144440, 334752, 760456, 1696464, 3722224, 8043040, 17135624, 36031104, 74840568, 153680928, 312198160, 627828272, 1250540024, 2468443296, 4830809868, 9377190336, 18061370288, 34531009760, 65552873736

Offset: 0

Vaclav Kotesovec, Aug 17 2015

nonn

Convolution of A161870 and A026011.

Seiichi Manyama, Table of n, a(n) for n = 0..10000
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. 19.

Cf. A015128, A026011, A156616, A161870, A216406, A261384, A261389.

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

a(n) ~ exp(1/6 + 3/2*(7*Zeta(3))^(1/3) * n^(2/3)) * (7*Zeta(3))^(2/9) / (A^2 * 2^(2/3) * n^(13/18) * sqrt(3*Pi)), where Zeta(3) = A002117 and A = A074962 is the Glaisher-Kinkelin constant.

G.f.: exp(Sum_{k>=1} (sigma_2(2*k) - sigma_2(k))*x^k/k). - Ilya Gutkovskiy, Apr 14 2019

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

Original entry on oeis.org

1, 2, 18, 88, 398, 1768, 7508, 30644, 121310, 467234, 1756080, 6457168, 23274788, 82381584, 286760344, 982874120, 3320800590, 11070619228, 36446345198, 118581503192, 381552358872, 1214868568728, 3829841265428, 11959828895612, 37013411304892, 113570015855642

Offset: 0

Paul D. Hanna, Feb 12 2012

nonn

Convolution of A023872 and A248882. - Vaclav Kotesovec, Aug 19 2015

			G.f.: A(x) = 1 + 2*x + 18*x^2 + 88*x^3 + 398*x^4 + 1768*x^5 + 7508*x^6 +...
where A(x) = (1+x)/(1-x) * (1+x^2)^8/(1-x^2)^8 * (1+x^3)^27/(1-x^3)^27 *...
Also, A(x) = Euler transform of [2,15,54,120,250,405,686,960,1458,...]:
A(x) = 1/((1-x)^2*(1-x^2)^15*(1-x^3)^54*(1-x^4)^120*(1-x^5)^250*(1-x^6)^405*...).

Seiichi Manyama, Table of n, a(n) for n = 0..4595 (terms 0..1000 from Vaclav Kotesovec)
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. 23.

Cf. A156616, A206622, A206624, A001159 (sigma_4).

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

{a(n)=polcoeff(prod(m=1,n+1,((1+x^m)/(1-x^m+x*O(x^n)))^(m^3)),n)}

{a(n)=polcoeff(exp(sum(m=1, n, (sigma(2*m, 4)-sigma(m, 4))/8*x^m/m)+x*O(x^n)), n)}

{a(n)=local(InvEulerGF=x*(2+15*x+46*x^2+60*x^3+46*x^4+15*x^5+2*x^6)/(1-x^2+x*O(x^n))^4);polcoeff(1/prod(k=1,n,(1-x^k+x*O(x^n))^polcoeff(InvEulerGF,k)),n)}
for(n=0,30,print1(a(n),", "))

G.f.: exp( Sum_{n>=1} (sigma_4(2*n) - sigma_4(n))/8 * x^n/n ), where sigma_4(n) is the sum of 4th powers of divisors of n (A001159).
Inverse Euler transform has g.f.: x*(2 + 15*x + 46*x^2 + 60*x^3 + 46*x^4 + 15*x^5 + 2*x^6)/(1-x^2)^4.
a(n) ~ (93*Zeta(5))^(59/600) * exp(5/4 * (93*Zeta(5)/2)^(1/5) * n^(4/5) + Zeta'(-3)) / (2^(59/100) * sqrt(5*Pi) * n^(359/600)), where Zeta(5) = A013663, Zeta'(-3) = A259068. - Vaclav Kotesovec, Aug 19 2015

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

Original entry on oeis.org

1, 2, 34, 228, 1414, 8872, 52876, 301136, 1662614, 8929406, 46738920, 239036116, 1197187780, 5882369976, 28397283056, 134864166352, 630819797174, 2908948327780, 13236421303742, 59477002686404, 264104800719672, 1159649708139680, 5037895127964316

Offset: 0

Paul D. Hanna, Feb 12 2012

nonn

Convolution of A023873 and A248883. - Vaclav Kotesovec, Aug 19 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
If m is even and m >= 2, then can be simplified as: a(n) ~ ((2^(m+2)-1) * Gamma(m+2) * Zeta(m+2) / (2^(2*m+3) * n))^(1/(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)) + (-1)^(m/2) * Gamma(m+1) * Zeta(m+1) / (2^(m+1) * Pi^m)) / sqrt((m+2)*Pi*n). - Vaclav Kotesovec, Aug 19 2015

			G.f.: A(x) = 1 + 2*x + 18*x^2 + 88*x^3 + 398*x^4 + 1768*x^5 + 7508*x^6 +...
where A(x) = (1+x)/(1-x) * (1+x^2)^16/(1-x^2)^16 * (1+x^3)^81/(1-x^3)^81 *...
Also, A(x) = Euler transform of [2,31,162,496,1250,2511,4802,7936,...]:
A(x) = 1/((1-x)^2*(1-x^2)^31*(1-x^3)^162*(1-x^4)^496*(1-x^5)^1250*(1-x^6)^2511*...).

Seiichi Manyama, Table of n, a(n) for n = 0..2916 (terms 0..1000 from Vaclav Kotesovec)
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. 23.

Cf. A015128 (m=0), A156616 (m=1), A206622 (m=2), A206623 (m=3), A001160 (sigma_5).

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

{a(n)=polcoeff(prod(m=1,n+1,((1+x^m)/(1-x^m+x*O(x^n)))^(m^4)),n)}

{a(n)=polcoeff(exp(sum(m=1, n, (sigma(2*m, 5)-sigma(m, 5))/16*x^m/m)+x*O(x^n)), n)}

{a(n)=local(InvEulerGF=x*(2+31*x+152*x^2+341*x^3+460*x^4+341*x^5+152*x^6+31*x^7+2*x^8)/(1-x^2+x*O(x^n))^5); polcoeff(1/prod(k=1,n,(1-x^k+x*O(x^n))^polcoeff(InvEulerGF,k)),n)}
for(n=0,30,print1(a(n),", "))

G.f.: exp( Sum_{n>=1} (sigma_5(2*n) - sigma_5(n))/16 * x^n/n ), where sigma_5(n) is the sum of 5th powers of divisors of n (A001160).
Inverse Euler transform has g.f.: x*(2 + 31*x + 152*x^2 + 341*x^3 + 460*x^4 + 341*x^5 + 152*x^6 + 31*x^7 + 2*x^8)/(1-x^2)^5.
a(n) ~ exp(3*2^(2/3)*Pi*n^(5/6)/5 + 3*Zeta(5)/(4*Pi^4)) / (2^(7/6) * 3^(1/2) * n^(7/12)), where Zeta(5) = A013663. - Vaclav Kotesovec, Aug 19 2015
a(0) = 1, a(n) = (2/n)*Sum_{k=1..n} A096960(k)*a(n-k) for n > 0. - Seiichi Manyama, Apr 30 2017

cp's OEIS Frontend

A294504 Binomial transform of A156616.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Formula

A294505 Inverse binomial transform of A156616.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Formula

A015128 Number of overpartitions of n: an overpartition of n is an ordered sequence of nonincreasing integers that sum to n, where the first occurrence of each integer may be overlined.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Julia

Maple

Mathematica

PARI

PARI

PARI

PARI

SageMath

Formula

Extensions

A073592 Euler transform of negative integers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

SageMath

Formula

A270919 Coefficient of x^n in Product_{k>=1} ((1 + x^k) / (1 - x^k))^n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

PARI

Formula

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

Views

Author

Keywords

Comments