cp's OEIS Frontend

A326985 G.f.: B(x)*B(2*x^2)*B(3*x^3)*..., where B(x) is g.f. of A000312.

%I A326985 #8 Aug 23 2019 18:21:42
%S A326985 1,1,6,32,287,3222,47606,831488,16890792,389286222,10037183606,
%T A326985 286154919078,8937624574652,303483905672078,11130904101218094,
%U A326985 438532313635906858,18470060947222927499,828155619735377936654,39384843256547964375436,1980138439071577626157382
%N A326985 G.f.: B(x)*B(2*x^2)*B(3*x^3)*..., where B(x) is g.f. of A000312.
%H A326985 Alois P. Heinz, <a href="/A326985/b326985.txt">Table of n, a(n) for n = 0..386</a>
%F A326985 a(n) ~ n^n.
%p A326985 B:= proc(n) option remember; n^n end:
%p A326985 b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i=1,
%p A326985       B(n), add(b(j, 1)*i^j*b(n-i*j, i-1), j=0..n/i)))
%p A326985     end:
%p A326985 a:= n-> b(n$2):
%p A326985 seq(a(n), n=0..20);  # _Alois P. Heinz_, Aug 23 2019
%t A326985 nmax = 20; CoefficientList[Series[Product[1+Sum[k^k*j^k*x^(j*k), {k, 1, nmax/j}], {j, 1, nmax}], {x, 0, nmax}], x]
%Y A326985 Cf. A000312, A110143, A326986.
%K A326985 nonn
%O A326985 0,3
%A A326985 _Vaclav Kotesovec_, Aug 10 2019