|
%I #29 Mar 17 2024 10:45:48
%S 1,1,4,14,51,188,702,2644,10026,38223,146359,562456,2168134,8379539,
%T 32459199,125984039,489837300,1907490728,7438346255,29042470132,
%U 113522618066,444199913556,1739735079466,6819657196928,26753893533257,105034060120469,412637434996367
%N Number of partitions of n with n kinds of 1.
%H Alois P. Heinz, <a href="/A292463/b292463.txt">Table of n, a(n) for n = 0..1664</a>
%F a(n) = [x^n] 1/(1-x)^n * 1/Product_{j=2..n} (1-x^j).
%F a(n) is n-th term of the Euler transform of n,1,1,1,... .
%F a(n) ~ c * 4^n / sqrt(n), where c = QPochhammer[-1, 1/2] / (8*sqrt(Pi) * QPochhammer[1/4, 1/4]) = 0.48841139329043831428669851139824427133317... - _Vaclav Kotesovec_, Sep 19 2017
%F Equivalently, c = 1/(4*sqrt(Pi)*QPochhammer(1/2)). - _Vaclav Kotesovec_, Mar 17 2024
%e a(2) = 4: 2, 1a1a, 1a1b, 1b1b.
%p b:= proc(n, i, k) option remember; `if`(n=0 or i<2,
%p binomial(k+n-1, n), add(b(n-i*j, i-1, k), j=0..n/i))
%p end:
%p a:= n-> b(n$3):
%p seq(a(n), n=0..30);
%p # second Maple program:
%p b:= proc(n, k) option remember; `if`(n=0, 1, add(
%p (numtheory[sigma](j)+k-1)*b(n-j, k), j=1..n)/n)
%p end:
%p a:= n-> b(n$2):
%p seq(a(n), n=0..30);
%p # third Maple program:
%p b:= proc(n, k) option remember; `if`(n=0, 1, `if`(k=1,
%p combinat[numbpart](n), b(n-1, k) +b(n, k-1)))
%p end:
%p a:= n-> b(n$2):
%p seq(a(n), n=0..30);
%t Table[SeriesCoefficient[1/(1-x)^(n-1) * Product[1/(1-x^k), {k,1,n}], {x,0,n}], {n,0,30}] (* _Vaclav Kotesovec_, Sep 19 2017 *)
%Y Main diagonal of A292508.
%Y Cf. A000041, A014329, A292462, A292503, A292507, A292541.
%K nonn
%O 0,3
%A _Alois P. Heinz_, Sep 16 2017
|
