A305101 G.f.: Sum_{k>=1} x^k/(1+x^k) * Product_{k>=1} (1+x^k)/(1-x^k).
0, 1, 2, 6, 11, 22, 40, 70, 116, 191, 304, 474, 726, 1094, 1624, 2384, 3453, 4950, 7030, 9890, 13798, 19108, 26264, 35858, 48652, 65615, 87996, 117396, 155826, 205854, 270728, 354506, 462306, 600544, 777184, 1002180, 1287889, 1649578, 2106152, 2680924
Offset: 0
Keywords
Programs
-
Mathematica
nmax = 40; CoefficientList[Series[Sum[x^k/(1+x^k), {k, 1, nmax}] * Product[(1+x^k)/(1-x^k), {k, 1, nmax}], {x, 0, nmax}], x]
-
PARI
my(N=44, q='q+O('q^N)); Vec( prod(k=1,N, (1+q^k)/(1-q^k)) * sum(k=1,N, 1*q^k/(1+q^k)) ) \\ Joerg Arndt, Jun 18 2020
Comments