A243148 - OEIS

%I #57 Jul 25 2022 08:33:15

%S 1,0,1,0,0,1,0,0,0,1,0,1,0,0,1,0,0,1,0,0,1,0,0,0,1,0,0,1,0,0,0,0,1,0,

%T 0,1,0,0,1,0,0,1,0,0,1,0,1,0,1,0,0,1,0,0,1,0,0,1,0,1,0,0,1,0,0,1,0,0,

%U 0,1,0,1,0,0,1,0,0,1,0,0,0,1,1,0,1,0,0,1,0,0,1,0,0,1,0,1,1,0,1,0,0,1,0,0,1,0,0,0,1,0,1,1,0,1,0,0,1,0,0,1

%N Triangle read by rows: T(n,k) = number of partitions of n into k nonzero squares; n >= 0, 0 <= k <= n.

%H Alois P. Heinz, <a href="/A243148/b243148.txt">Rows n = 0..200, flattened</a>

%F T(n,k) = [x^n y^k] 1/Product_{j>=1} (1-y*x^A000290(j)).

%F Sum_{k=1..n} k * T(n,k) = A281541(n).

%F Sum_{k=1..n} n * T(n,k) = A276559(n).

%F Sum_{k=0..n} (-1)^k * T(n,k) = A292520(n).

%e T(20,5) = 2 = #{ (16,1,1,1,1), (4,4,4,4,4) } since 20 = 4^2 + 4 * 1^2 = 5 * 2^2.

%e Triangle T(n,k) begins:

%e   1;

%e   0, 1;

%e   0, 0, 1;

%e   0, 0, 0, 1;

%e   0, 1, 0, 0, 1;

%e   0, 0, 1, 0, 0, 1;

%e   0, 0, 0, 1, 0, 0, 1;

%e   0, 0, 0, 0, 1, 0, 0, 1;

%e   0, 0, 1, 0, 0, 1, 0, 0, 1;

%e   0, 1, 0, 1, 0, 0, 1, 0, 0, 1;

%e   0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1;

%e   0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1;

%e   0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 1;

%e   (...)

%p b:= proc(n, i, t) option remember; `if`(n=0, `if`(t=0, 1, 0),

%p       `if`(i<1 or t<1, 0, b(n, i-1, t)+

%p       `if`(i^2>n, 0, b(n-i^2, i, t-1))))

%p     end:

%p T:= (n, k)-> b(n, isqrt(n), k):

%p seq(seq(T(n, k), k=0..n), n=0..14);

%p # second Maple program:

%p b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,

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

%p     end:

%p T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(n, isqrt(n))):

%p seq(T(n), n=0..14);  # _Alois P. Heinz_, Oct 30 2021

%t b[n_, i_, k_, t_] := b[n, i, k, t] = If[n == 0, If[t == 0, 1, 0], If[i < 1 || t < 1, 0, b[n, i-1, k, t] + If[i^2 > n, 0, b[n-i^2, i, k, t-1]]]]; T[n_, k_] := b[n, Sqrt[n] // Floor, k, k]; Table[Table[T[n, k], {k, 0, n}], {n, 0, 14}] // Flatten (* _Jean-François Alcover_, Jun 06 2014, after _Alois P. Heinz_ *)

%t T[n_, k_] := Count[PowersRepresentations[n, k, 2], r_ /; FreeQ[r, 0]]; T[0, 0] = 1; Table[T[n, k], {n, 0, 14}, {k, 0, n}] // Flatten (* _Jean-François Alcover_, Feb 19 2016 *)

%o (PARI) T(n,k,L=n)=if(n>k*L^2, 0, k>n-3, k==n, k<2, issquare(n,&n) && n<=L*k, k>n-6, n-k==3, L=min(L,sqrtint(n-k+1)); sum(r=0,min(n\L^2,k-1),T(n-r*L^2,k-r,L-1), n==k*L^2)) \\ _M. F. Hasler_, Aug 03 2020

%Y Columns k = 0..10 give: A000007, A010052 (for n>0), A025426, A025427, A025428, A025429, A025430, A025431, A025432, A025433, A025434.

%Y Row sums give A001156.

%Y T(2n,n) gives A111178.

%Y T(n^2,n) gives A319435.

%Y Cf. A000290, A276559, A281541, A292520, A341040.

%Y T(n,k) = 1 for n in A025284, A025321, A025357, A294675, A295670, A295797 (for k = 2..7, respectively).

%K nonn,tabl

%O 0,216

%A _Alois P. Heinz_, May 30 2014