OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
FORMULA
Expansion of q^(1/24) * eta(q^4)^5 / (eta(q)^5 * eta(q^8)^2) in powers of q. - Michael Somos, Sep 24 2011
a(n) = A000712(2*n).
EXAMPLE
G.f.: 1 + 5*x + 20*x^2 + 65*x^3 + 185*x^4 + 481*x^5 + 1165*x^6 + 2665*x^7 + ...
G.f.: 1/q + 5*q^23 + 20*q^47 + 65*q^71 + 185*q^95 + 481*q^119 + 1165*q^143 + ...
a(1)=5 because we have 2, 2', 11, 1'1 and 1'1'.
MAPLE
with(combinat): A000712:=n-> add(numbpart(k)*numbpart(n-k), k=0..n): seq(A000712(2*n), n=0..32); # Emeric Deutsch, Dec 16 2004
MATHEMATICA
a[n_] := Sum[PartitionsP[k] PartitionsP[2 n - k], {k, 0, 2 n}]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Nov 30 2015, adapted from Maple *)
PROG
(PARI) {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^4 + A)^5 / (eta(x + A)^5 * eta(x^8 + A)^2), n))} /* Michael Somos, Sep 24 2011 */
(PARI) {a(n) = local(A); if( n<0, 0, n = 2*n; A = x * O(x^n); polcoeff( 1 / eta(x + A)^2, n))} /* Michael Somos, Sep 24 2011 */
(Magma)
m:=40;
f:= func< x | (&*[ (1-x^(4*n))^5/((1-x^n)^5*(1-x^(8*n))^2) : n in [1..m+2]]) >;
R<x>:=PowerSeriesRing(Rationals(), m);
Coefficients(R!( f(x) )); // G. C. Greubel, Mar 27 2023
(SageMath)
m=40
def f(x): return product( (1-x^(4*n))^5/((1-x^n)^5*(1-x^(8*n))^2) for n in range(1, m+2) )
def A100535_list(prec):
P.<x> = PowerSeriesRing(QQ, prec)
return P( f(x) ).list()
A100535_list(m) # G. C. Greubel, Mar 27 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Nov 27 2004
EXTENSIONS
More terms from Emeric Deutsch, Dec 16 2004
STATUS
approved