%I #14 Feb 19 2019 09:21:18
%S 1,0,1,0,0,1,0,1,0,0,0,1,0,1,0,0,0,1,0,0,0,0,1,0,1,0,1,0,1,0,0,0,0,0,
%T 1,0,0,0,0,0,1,0,1,0,1,0,1,0,1,0,0,0,1,0,0,0,0,0,0,0,1,0,0,0,0,0,0,1,
%U 0,1,0,1,0,2,0,1,0,0,0,2,0,1,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,1,0,1,0,1,0,2,0,1,0,1,0,2,0,1,0,0,0,1,0,1,0,1,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,1
%N Triangle read by rows: T(n,k) is the number of partitions of n for which the sum of the squares of the parts is k (n>=0, k>=0).
%C Number of entries in row n = 1 + n^2.
%C Sum of entries in row n = A000041(n).
%C Sum(k*T(n,k), k>=0) = A066183(n).
%H Guo-Niu Han, <a href="http://arxiv.org/abs/0804.1849">An explicit expansion formula for the powers of the Euler Product in terms of partition hook lengths</a>, arXiv:0804.1849 [math.CO], 2008.
%F G.f.: G(t,x) = 1/Product_{k>=1} (1 - t^{k^2}*x^k).
%e Row 3 is 0,0,0,1,0,1,0,0,0,1 because in the partitions of 3, namely [1,1,1], [2,1], [3], the sums of the squares of the parts are 3, 5, and 9, respectively.
%e Triangle starts:
%e 1;
%e 0,1;
%e 0,0,1,0,1;
%e 0,0,0,1,0,1,0,0,0,1;
%e 0,0,0,0,1,0,1,0,1,0,1,0,0,0,0,0,1.
%p g := 1/(product(1-t^(k^2)*x^k, k = 1 .. 100)): gser := simplify(series(g, x = 0, 15)): for n from 0 to 8 do P[n] := sort(coeff(gser, x, n)) end do: for n from 0 to 8 do seq(coeff(P[n], t, j), j = 0 .. n^2) end do; # yields sequence in triangular form
%t m = 8; CoefficientList[#, t]& /@ CoefficientList[1/Product[(1 - t^(k^2)* x^k), {k, 1, m}] + O[x]^m, x] // Flatten (* _Jean-François Alcover_, Feb 19 2019 *)
%Y Cf. A000041, A066183, A229325 - A229332, A264402, A265247 - A265253.
%K nonn,tabf
%O 0,74
%A _Emeric Deutsch_, Dec 06 2015
|