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

%I #11 Sep 19 2018 10:48:21

%S 1,2,3,8,10,12,14,24,36,40,44,60,78,84,90,128,153,180,190,240,273,308,

%T 322,384,475,520,567,644,754,810,868,992,1122,1258,1330,1548,1702,

%U 1862,1950,2200,2460,2646,2838,3124,3510,3726,3948,4320,4802,5200,5457,6032,6572,7128,7425,8064,8778,9454,9971,10680

%N Expansion of Sum_{k>=1} k^2*x^k^2/(1 - x^k^2) * Product_{k>=1} 1/(1 - x^k^2).

%C Sum of all parts of all partitions of n into squares.

%C Convolution of the sequences A001156 and A035316.

%H Alois P. Heinz, <a href="/A276559/b276559.txt">Table of n, a(n) for n = 1..10000</a>

%H <a href="/index/Par#partN">Index entries for related partition-counting sequences</a>

%F G.f.: Sum_{k>=1} k^2*x^k^2/(1 - x^k^2) * Product_{k>=1} 1/(1 - x^k^2).

%F G.f.: x*f'(x), where f(x) = Product_{k>=1} 1/(1 - x^k^2).

%F a(n) = n * A001156(n).

%F a(n) = n * Sum_{k=1..n} A243148(n,k). - _Alois P. Heinz_, Sep 19 2018

%e a(8) = 24 because we have [4, 4], [4, 1, 1, 1, 1], [1, 1, 1, 1, 1, 1, 1, 1] and 3*8 = 24.

%p b:= proc(n, i) option remember; `if`(n=0 or i=1, [1, n], (s->

%p `if`(s>n, 0, (p->p+[0, p[1]*s])(b(n-s, i))))(i^2)+b(n, i-1))

%p end:

%p a:= n-> b(n, isqrt(n))[2]:

%p seq(a(n), n=1..70); # _Alois P. Heinz_, Sep 19 2018

%t nmax = 60; Rest[CoefficientList[Series[Sum[k^2 x^k^2/(1 - x^k^2), {k, 1, nmax}] Product[1/(1 - x^k^2), {k, 1, nmax}], {x, 0, nmax}], x]]

%t nmax = 60; Rest[CoefficientList[Series[x D[Product[1/(1 - x^k^2), {k, 1, nmax}], x], {x, 0, nmax}], x]]

%Y Cf. A000290, A001156, A035316, A066186, A243148, A281541.

%K nonn

%O 1,2

%A _Ilya Gutkovskiy_, Apr 10 2017