A128196 - OEIS

%I #9 Nov 16 2022 05:06:01

%S 1,3,13,73,527,4775,52589,683785,10257031,174370039,3313031765,

%T 69573669113,1600194393695,40004859850567,1080131215981693,

%U 31323805263501865,971037963168623351,32044252784564701655,1121548847459764820069,41497307356011298866841,1618394986884440656855375

%N a(n) = (2*n - 1)*a(n - 1) + 2^n for n >= 1, a(0) = 1.

%C A weighted sum of quotients of double factorials.

%C a(n) are the row sum of triangle A126063.

%H P. Luschny, <a href="http://www.luschny.de/math/seq/variations.html">Variants of Variations</a>.

%F a(n) = (2n)!/(n! 2^n) Sum(k=0..n, 4^k k!/(2k)!)

%F a(n) = 2^n Gamma(n+1/2) Sum(k=0..n, 1/Gamma(k+1/2))

%F a(n) = Sum(k=0..n, 2^k n!!/k!!) [n!! defined as A001147(n), Gottfried Helms]

%F a(n) = Sum(k=0..n, 2^(2k-n)((n+1)! Catalan(n))/((k+1)! Catalan(k))) [Catalan(n) A000108]

%F a(n) = Sum(k=0..n, 2^(2k-n) QuadFact(n)/QuadFact(k)) [QuadFact(n) A001813]

%F a(n) = Sum(k=0..n, 2^(2k-n) (-1)^(n-k) A097388(n)/A097388(k) )

%F a(n) = A001147(n) Sum(k=0..n, 2^k / A001147(k))

%F a(n) = A128195(n)/A005408(n)

%F a(n) = A128195(n-1)+A000079(n) (if n>0)

%F Recursive form: a(n) = (2n-1)*a(n-1) + 2^n; a(0) = 1 [Gottfried Helms]

%F Note: The following constants will be used in the next formulas.

%F K = (1-exp(1)*Gamma(1/2,1))/Gamma(1/2)

%F M = sqrt(2)(1+exp(1)(Gamma(1/2)-Gamma(1/2,1)))

%F Generalized form: For x>0

%F a(x) = 2^x(exp(1)*Gamma(x+1/2,1) + K*Gamma(x+1/2))

%F Asymptotic formula:

%F a(n) ~ 2^n*(1+(exp(1)+K)*(n-1/2)!)

%F a(n) ~ M(2exp(-1)(n-1/(24*n+19/10*1/n)))^n

%p a := n -> `if`(n=0,1,(2*n-1)*a(n-1)+2^n);

%t a[n_] := Sum[2^k*((2*n-1)!!/(2*k-1)!!), {k, 0, n}]; Table[a[n], {n, 0, 14}] (* _Jean-François Alcover_, Jun 28 2013 *)

%Y Cf. A128195, A001147, A126063.

%K easy,nonn

%O 0,2

%A _Peter Luschny_, Feb 26 2007