|
%I #17 Apr 17 2019 08:54:16
%S 1,2,2,2,3,4,4,4,4,5,6,6,6,7,8,8,9,10,10,10,11,12,12,12,12,14,16,16,
%T 16,18,20,20,20,20,21,22,23,24,25,26,27,29,30,30,30,32,34,34,34,36,39,
%U 40,41,43,45,46,47,48,49,50,50,52,55,56,57,61,64,64,65
%N Number of sets of nonzero squares with sum <= n
%C Partial sums of A033461.
%H Robert Israel, <a href="/A248801/b248801.txt">Table of n, a(n) for n = 0..10000</a>
%F G.f.: (1-x)^(-1) * product(k>=1, 1 + x^(k^2)).
%e For n=5 the sets are {}, {1^2}, {2^2}, {1^2, 2^2} so a(5) = 4.
%p N:= 200: # to get a(0) to a(N)
%p g:= (1-x)^(-1)*mul(1 + x^(m^2), m=1 .. floor(sqrt(N))):
%p S:= series(g, x, N+1):
%p seq(coeff(S,x,j),j=0..N);
%t CoefficientList[Series[(1 - x)^(-1) Product[1 + x^(k^2), {k, 50}], {x, 0, 50}], x] (* _Vincenzo Librandi_, Oct 15 2014 *)
%t b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i^2>n, 0, b[n-i^2, i-1]]]]; Table[b[n, Floor[Sqrt[n]]], {n, 0, 100}] // Accumulate (* _Jean-François Alcover_, Apr 17 2019, after _Alois P. Heinz_ in A033461 *)
%Y Cf. A033461.
%K nonn
%O 0,2
%A _Robert Israel_, Oct 14 2014
|
