%I #10 Mar 23 2018 05:31:55
%S 1,1,3,10,39,151,588,2304,9111,36307,145553,586246,2370264,9614242,
%T 39105580,159444160,651468967,2666771488,10934393619,44899828056,
%U 184616878289,760010818689,3132147583744,12921037206764,53351800567200,220478125956426,911839751015196,3773836780169050
%N a(n) = [x^n] Product_{k>=1} 1/(1 - x^(k^2))^n.
%C Number of partitions of n into squares of n kinds.
%H Vaclav Kotesovec, <a href="/A300974/b300974.txt">Table of n, a(n) for n = 0..500</a>
%H <a href="/index/Par#part">Index entries for sequences related to partitions</a>
%H <a href="/index/Su#ssq">Index entries for sequences related to sums of squares</a>
%F From _Vaclav Kotesovec_, Mar 23 2018: (Start)
%F a(n) ~ c * d^n / sqrt(n), where
%F d = 4.216358447600641565890184638418336163396695730036... and
%F c = 0.26442245016754864773722176155288663999776... (End)
%p a:= proc(m) option remember; local b; b:= proc(n, i)
%p option remember; `if`(n=0, 1, `if`(i<1, 0, add(
%p binomial(m+j-1, j)*b(n-i^2*j, i-1), j=0..n/i^2)))
%p end: b(n, isqrt(n))
%p end:
%p seq(a(n), n=0..30); # _Alois P. Heinz_, Mar 17 2018
%t Table[SeriesCoefficient[Product[1/(1 - x^k^2)^n, {k, 1, n}], {x, 0, n}], {n, 0, 27}]
%Y Cf. A000290, A001156, A008485, A023871, A240944, A279225, A285047, A300975, A301518.
%K nonn
%O 0,3
%A _Ilya Gutkovskiy_, Mar 17 2018
|