|
%I #9 Aug 16 2019 12:52:09
%S 1,1,4,10,29,62,176,363,931,2029,4751,10062,23749,48959,109342,230981,
%T 500344,1031667,2223218,4531585,9570395,19523510,40411313,81628389,
%U 168484616,336850254,685112670,1369559157,2757908932,5464925114,10958578421,21574592680
%N Expansion of Product_{k>=1} 1/(1 - (k*(k + 1)/2)*x^k).
%H Alois P. Heinz, <a href="/A305204/b305204.txt">Table of n, a(n) for n = 0..3841</a>
%F G.f.: Product_{k>=1} 1/(1 - A000217(k)*x^k).
%F G.f.: exp(Sum_{k>=1} Sum_{j>=1} (j*(j + 1))^k*x^(j*k)/(k*2^k)).
%p b:= proc(n, i) option remember; `if`(n=0 or i=1,
%p 1, b(n, i-1)+(1+i)*i/2*b(n-i, min(n-i, i)))
%p end:
%p a:= n-> b(n$2):
%p seq(a(n), n=0..33); # _Alois P. Heinz_, Aug 16 2019
%t nmax = 31; CoefficientList[Series[Product[1/(1 - (k (k + 1)/2) x^k), {k, 1, nmax}], {x, 0, nmax}], x]
%t nmax = 31; CoefficientList[Series[Exp[Sum[Sum[(j (j + 1))^k x^(j k)/(k 2^k), {j, 1, nmax}], {k, 1, nmax}]], {x, 0, nmax}], x]
%t a[n_] := a[n] = If[n == 0, 1, Sum[Sum[d^(k/d + 1) ((d + 1)/2)^(k/d), {d, Divisors[k]}] a[n - k], {k, 1, n}]/n]; Table[a[n], {n, 0, 31}]
%Y Cf. A000217, A000294, A006906, A077335, A265836.
%K nonn
%O 0,3
%A _Ilya Gutkovskiy_, May 27 2018
|
