%I M5308 #33 May 08 2018 15:11:54
%S 1,1,53,195,22999,29944523,109535241009,29404527905795295658,
%T 455377030420113432210116914702,
%U 26370812569397719001931992945645578779849,152537496709054809881638897472985990866753853122697839,100043420063777451042472529806266909090824649341814868347109676190691
%N Numerators of numbers occurring in continued fraction connected with expansion of gamma function.
%D M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math.Series 55, Tenth Printing, 1972, p. 258.
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%D H. S. Wall, Analytic Theory of Continued Fractions, Chelsea 1973, p. 365.
%H Vincenzo Librandi, <a href="/A005146/b005146.txt">Table of n, a(n) for n = 0..30</a>
%H M. Abramowitz and I. A. Stegun, eds., <a href="http://www.convertit.com/Go/ConvertIt/Reference/AMS55.ASP">Handbook of Mathematical Functions</a>, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
%H M. Abramowitz and I. A. Stegun, eds., <a href="http://www.convertit.com/Go/ConvertIt/Reference/AMS55.ASP">Handbook of Mathematical Functions</a>, National Bureau of Standards Applied Math.Series 55, Tenth Printing, 1972, p. 258.
%H B. W. Char, <a href="http://www.jstor.org/stable/2006103">On Stieltjes' continued fraction for the gamma function</a>, Math. Comp., 34 (1980), 547-551.
%H Peter Luschny, <a href="/A005146/a005146.txt">Maple program for A005146/A005147</a>
%H Peter Luschny, <a href="http://www.luschny.de/math/factorial/approx/continuedfraction.html">Continued fraction</a>
%H W. Wang, <a href="http://dx.doi.org/10.1016/j.jnt.2015.12.016">Unified approaches to the approximations of the gamma function</a>, J. Number Theory (2016).
%t len = 12; s[p_] := (-1)^p * BernoulliB[2p+2]/(2p+1)/(2p+2); Do[m[n, 1] = 0, {n, 0, len}]; Do[m[n, 2] = s[n+1]/s[n], {n, 0, len-1}]; Do[m[n, k] =
%t If[OddQ[k], m[n+1, k-2]+m[n+1, k-1]-m[n, k-1],
%t m[n+1, k-2]*m[n+1, k-1]/m[n, k-1]], {k, 3, len}, {n, 0,
%t len-k+1}]; Do[m[n, 1] = s[n], {n, 0, len}];
%t Table[m[0, k], {k, 1, len}] // Numerator
%t (* _Jean-François Alcover_, May 24 2011, after _Peter Luschny_ *)
%Y Cf. A005147.
%K nonn,frac,nice
%O 0,3
%A _Simon Plouffe_ and _N. J. A. Sloane_
%E More terms from _Rainer Rosenthal_, Jan 11 2007