%I #18 Oct 05 2024 09:40:24
%S 1,9,27,73,368,343,585,3825,6615,3375,4681,36394,88536,86614,29791,
%T 37449,332883,1024002,1449198,970677,250047,299593,2979420,10970133,
%U 20078192,19714083,9974580,2048383,2396745,26298405,112122225,250021125,320944275,239783895,97221555,16581375
%N Triangle, read by rows, where T(n,k) equals the sum of cubes of numbers < 2^n having exactly k ones in their binary expansion.
%C Compare to triangle A110200 (sum of squares).
%H Paul D. Hanna, <a href="/A110205/b110205.txt">Rows n = 1..45, flattened.</a>
%F T(n, k) = (8^n-1)/7*C(n-3, k-1) + ((2^n-1)*(4^n-1)-(8^n-1)/7)*C(n-3, k-2) + (2^n-1)^3*C(n-3, k-3).
%F G.f. for row n: ((8^n-1)/7 + ((2^n-1)*(4^n-1)-(8^n-1)/7)*x + (2^n-1)^3*x^2)*(1+x)^(n-3).
%e Row 4 is formed by sums of cubes of numbers < 2^4:
%e T(4,1) = 1^3 + 2^3 + 4^3 + 8^3 = 585;
%e T(4,2) = 3^3 + 5^3 + 6^3 + 9^3 + 10^3 + 12^3 = 3825;
%e T(4,3) = 7^3 + 11^3 + 13^3 + 14^3 = 6615;
%e T(4,4) = 15^3 = 3375.
%e Triangle begins:
%e 1;
%e 9, 27;
%e 73, 368, 343;
%e 585, 3825, 6615, 3375;
%e 4681, 36394, 88536, 86614, 29791;
%e 37449, 332883, 1024002, 1449198, 970677, 250047;
%e 299593, 2979420, 10970133, 20078192, 19714083, 9974580, 2048383;
%e 2396745, 26298405, 112122225, 250021125, 320944275, 239783895, 97221555, 16581375; ...
%e Row g.f.s are:
%e row 1: (1 + 2*x + 1*x^2)/(1+x)^2;
%e row 2: (9 + 36*x + 27*x^2)/(1+x);
%e row 3: (73 + 368*x + 343*x^2);
%e row 4: (585 + 3240*x + 3375*x^2)*(1+x).
%e G.f. for row n is:
%e ((8^n-1)/7 + ((2^n-1)*(4^n-1)-(8^n-1)/7)*x + (2^n-1)^3*x^2)*(1+x)^(n-3).
%t b[n_, k_]= Binomial[n-3, k];
%t T[n_, k_]:= (8^n-1)/7*(b[n,k-1] -b[n,k-2]) + (2^n-1)^2*((2^n+1)*b[n,k-2] + (2^n-1)*b[n,k-3]);
%t A110205[n_, k_]:= If[n<3, T[n,k]/2, T[n,k]];
%t Table[A110205[n,k], {n,12}, {k,n}]//Flatten (* _G. C. Greubel_, Oct 03 2024 *)
%o (PARI) T(n,k)=(8^n-1)/7*binomial(n-3,k-1)+((2^n-1)*(4^n-1)-(8^n-1)/7)*binomial(n-3,k-2) +(2^n-1)^3*binomial(n-3,k-3)
%o (PARI) /* Sum of cubes of numbers<2^n with k 1-bits: */
%o T(n,k)=local(B=vector(n+1));if(n<k || k<1,0, for(m=1,2^n-1, B[1+sum(i=1,#binary(m),(binary(m))[i])]+=m^3);B[k+1])
%o (Magma)
%o b:= func< n,k | Binomial(n-3, k) >;
%o A110205:= func< n,k | (8^n-1)/7*(b(n,k-1) -b(n,k-2)) +(2^n-1)^2*((2^n+1)*b(n,k-2) +(2^n-1)*b(n,k-3)) >;
%o [A110205(n,k): k in [1..n], n in [1..12]]; // _G. C. Greubel_, Oct 03 2024
%o (SageMath)
%o def b(n,k): return binomial(n-3, k)
%o def A110205(n,k): return (8^n-1)/7*(b(n,k-1) - b(n,k-2)) + (2^n-1)^2*((2^n+1)*b(n,k-2) + (2^n-1)*b(n,k-3))
%o flatten([[A110205(n,k) for k in range(1,n+1)] for n in range(1,13)]) # _G. C. Greubel_, Oct 03 2024
%Y Cf. A110206 (row sums), A110207 (central terms), A023001 (column 1).
%K nonn,tabl
%O 1,2
%A _Paul D. Hanna_, Jul 16 2005