OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..400
FORMULA
E.g.f.: exp( Sum_{n>=1} (2*x)^n/n * Sum_{k=0..n-1} 1/C(n-1,k) ).
a(n) ~ n!*n*2^n * (1 - 2*log(n)/n). - Vaclav Kotesovec, Jun 27 2013
EXAMPLE
E.g.f.: A(x) = 1 + 2*x + 12*x^2/2! + 96*x^3/3! + 976*x^4/4! + 12000*x^5/5! +...
where the logarithm involves sums of reciprocal binomial coefficients:
log(A(x)) = 2*x*(1) + (2*x)^2/2*(1 + 1) + (2*x)^3/3*(1 + 1/2 + 1) + (2*x)^4/4*(1 + 1/3 + 1/3 + 1) + (2*x)^5/5*(1 + 1/4 + 1/6 + 1/4 + 1) + (2*x)^6/6*(1 + 1/5 + 1/10 + 1/10 + 1/5 + 1) +...
Explicitly, the logarithm begins:
log(A(x)) = 2*x + 8*x^2/2! + 40*x^3/3! + 256*x^4/4! + 2048*x^5/5! + 19968*x^6/6! +...
in which the coefficients equal 2*A126674(n).
MATHEMATICA
CoefficientList[Series[(1-2*x)^(-1/(1-x)), {x, 0, 20}], x]* Range[0, 20]! (* Vaclav Kotesovec, Jun 27 2013 *)
PROG
(PARI) {a(n)=n!*polcoeff(exp(sum(m=1, n, 2^m*x^m/m*sum(k=0, m-1, 1/binomial(m-1, k)))+x*O(x^n)), n)}
(PARI) {a(n)=n!*polcoeff((1-2*x+x*O(x^n))^(-1/(1-x)), n)}
(Magma)
m:=50;
f:= func< x | Exp((&+[(&+[ 1/Binomial(n-1, k): k in [0..n-1]])*(2*x)^n/n: n in [1..m+2]])) >;
R<x>:=PowerSeriesRing(Rationals(), m);
Coefficients(R!(Laplace( f(x) ))); // G. C. Greubel, Feb 02 2023
(SageMath)
m=50
def f(x): return exp(sum(sum( 1/binomial(n-1, k) for k in range(n))*(2*x)^n/n for n in range(1, m+2)))
def A193425_list(prec):
P.<x> = PowerSeriesRing(QQ, prec)
return P( f(x) ).egf_to_ogf().list()
A193425_list(m) # G. C. Greubel, Feb 02 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jul 27 2011
STATUS
approved