login
G.f.: x * Product_{j>=1, k>=1} 1/(1 - a(j)*x^(j*k)).
1

%I #11 May 14 2019 22:05:40

%S 1,1,3,7,19,47,134,357,1031,2912,8612,25007,75378,223884,683915,

%T 2067578,6376800,19503546,60749341,187592661,587938043,1831377952,

%U 5773159368,18092820941,57328904204,180657986051,574735018826,1820143698295,5810522774503,18473074695503

%N G.f.: x * Product_{j>=1, k>=1} 1/(1 - a(j)*x^(j*k)).

%p A:= proc(n) option remember; series(x*`if`(n=1, 1, mul(mul(

%p 1/(1-a(j)*x^(j*k)), k=1..(n-1)/j), j=1..n-1)), x, n+1)

%p end:

%p a:= n-> coeff(A(n), x, n):

%p seq(a(n), n=1..35); # _Alois P. Heinz_, May 14 2019

%t a[n_] := a[n] = SeriesCoefficient[x Product[Product[1/(1 - a[j] x^(j k)), {k, 1, n - 1}], {j, 1, n - 1}], {x, 0, n}]; Table[a[n], {n, 1, 30}]

%Y Cf. A093637, A179467, A308154.

%K nonn

%O 1,3

%A _Ilya Gutkovskiy_, May 14 2019