login
E.g.f.: -log(1 - x - Sum_{k>=2} prime(k-1) * x^k / k!).
1

%I #6 Aug 04 2021 16:41:44

%S 1,3,11,59,416,3683,39093,484220,6854176,109150227,1931303809,

%T 37589753206,798135918850,18358887315769,454779141016707,

%U 12070296596154136,341715021307029876,10278722402921420619,327369178071821161755,11005696560250745851048,389469699942038630639524

%N E.g.f.: -log(1 - x - Sum_{k>=2} prime(k-1) * x^k / k!).

%F a(n) = A008578(n) + (1/n) * Sum_{k=1..n-1} binomial(n,k) * A008578(n-k) * k * a(k).

%t nmax = 21; CoefficientList[Series[-Log[1 - x - Sum[Prime[k - 1] x^k/k!, {k, 2, nmax}]], {x, 0, nmax}], x] Range[0, nmax]! // Rest

%Y Cf. A000040, A007446, A007447, A008578, A300632, A336185, A343622.

%K nonn

%O 1,2

%A _Ilya Gutkovskiy_, Aug 04 2021