OFFSET
0,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..430
FORMULA
E.g.f.: exp( Sum_{n>=1} x^n/n * Sum_{k=0..n-1} 1/C(n-1,k) ).
a(n) ~ n! * (n-2*log(n)). - Vaclav Kotesovec, Jul 13 2014
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)) = x*(1) + x^2/2*(1 + 1) + x^3/3*(1 + 1/2 + 1) + x^4/4*(1 + 1/3 + 1/3 + 1) + x^5/5*(1 + 1/4 + 1/6 + 1/4 + 1) + x^6/6*(1 + 1/5 + 1/10 + 1/10 + 1/5 + 1) + ...
Explicitly, the logarithm begins:
log(A(x)) = x + 2*x^2/2! + 5*x^3/3! + 16*x^4/4! + 64*x^5/5! + 312*x^6/6! + 1812*x^7/7! + 12288*x^8/8! + ... + A003149(n-1)*x^n/n! + ...
MATHEMATICA
CoefficientList[Series[1/(1-x)^(2/(2-x)), {x, 0, 20}], x]*Range[0, 20]! (* Vaclav Kotesovec, Jul 13 2014 *)
PROG
(PARI) {a(n)=n!*polcoeff(exp(sum(m=1, n, x^m/m*sum(k=0, m-1, 1/binomial(m-1, k))) +x*O(x^n)), n)}
for(n=0, 25, print1(a(n), ", "))
(PARI) {a(n)=n!*polcoeff(1/(1-x+x*O(x^n))^(2/(2-x)), n)}
for(n=0, 25, print1(a(n), ", "))
(Magma) R<x>:=PowerSeriesRing(Rationals(), 50); Coefficients(R!(Laplace( 1/(1-x)^(2/(2-x)) ))); // G. C. Greubel, Jul 12 2023
(SageMath)
m=50
def f(x): return exp(sum(sum( 1/binomial(n-1, k) for k in range(n))*x^n/n for n in range(1, m+2)))
def A235802_list(prec):
P.<x> = PowerSeriesRing(QQ, prec)
return P( f(x) ).egf_to_ogf().list()
A235802_list(m) # G. C. Greubel, Jul 12 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jan 15 2014
STATUS
approved