cp's OEIS Frontend

A299768 Triangle read by rows: T(n,k) = sum of all squares of the parts k in all partitions of n, with n >= 1, 1 <= k <= n.

%I A299768 #47 Jun 07 2018 21:58:44
%S A299768 1,2,4,4,4,9,7,12,9,16,12,16,18,16,25,19,32,36,32,25,36,30,44,54,48,
%T A299768 50,36,49,45,76,81,96,75,72,49,64,67,104,135,128,125,108,98,64,81,97,
%U A299768 164,189,208,200,180,147,128,81,100,139,224,279,288,300,252,245,192,162,100,121
%N A299768 Triangle read by rows: T(n,k) = sum of all squares of the parts k in all partitions of n, with n >= 1, 1 <= k <= n.
%H A299768 Alois P. Heinz, <a href="/A299768/b299768.txt">Rows n = 1..200, flattened</a>
%F A299768 T(n,k) = (k^2)*A066633(n,k) = k*A138785(n,k). - _Omar E. Pol_, Jun 07 2018
%e A299768 Triangle begins:
%e A299768    1;
%e A299768    2,  4;
%e A299768    4,  4,  9;
%e A299768    7, 12,  9, 16;
%e A299768   12, 16, 18, 16, 25,
%e A299768   19, 32, 36, 32, 25, 36;
%e A299768   30, 44, 54, 48, 50, 36, 49;
%e A299768 ...
%e A299768 For n = 4 the partitions of 4 are [4], [2, 2], [3, 1], [2, 1, 1], [1, 1, 1, 1], so the squares of the parts are respectively [16], [4, 4], [9, 1], [4, 1, 1], [1, 1, 1, 1]. The sum of the squares of the parts 1 is 1 + 1 + 1 + 1 + 1 + 1 + 1 = 7. The sum of the squares of the parts 2 is 4 + 4 + 4 = 12. The sum of the squares of the parts 3 is 9. The sum of the squares of the parts 4 is 16. So the fourth row of triangle is [7, 12, 9, 16].
%p A299768 b:= proc(n, i) option remember; `if`(n=0 or i=1, 1+n*x, b(n, i-1)+
%p A299768       (p-> p+(coeff(p, x, 0)*i^2)*x^i)(b(n-i, min(n-i, i))))
%p A299768     end:
%p A299768 T:= n-> (p-> seq(coeff(p, x, i), i=1..n))(b(n$2)):
%p A299768 seq(T(n), n=1..14);  # _Alois P. Heinz_, Mar 20 2018
%t A299768 b[n_, i_] := b[n, i] = If[n == 0 || i == 1, 1 + n*x, b[n, i - 1] + # + (Coefficient[#, x, 0]*i^2*x^i)&[b[n - i, Min[n - i, i]]]];
%t A299768 T[n_] := Table[Coefficient[#, x, i], {i, 1, n}]&[b[n, n]];
%t A299768 Table[T[n], {n, 1, 14}] // Flatten (* _Jean-François Alcover_, May 22 2018, after _Alois P. Heinz_ *)
%o A299768 (PARI) row(n) = {v = vector(n); forpart(p=n, for(k=1, #p, v[p[k]] += p[k]^2;);); v;} \\ _Michel Marcus_, Mar 20 2018
%Y A299768 Column 1 is A000070.
%Y A299768 Leading diagonal is A000290, n >= 1.
%Y A299768 Row sums give A066183.
%Y A299768 Both A180681 and A206561 have the same row sums as this triangle.
%Y A299768 Cf. A066633, A138785.
%K A299768 nonn,tabl
%O A299768 1,2
%A A299768 _Omar E. Pol_, Mar 19 2018
%E A299768 More terms from _Michel Marcus_, Mar 20 2018