%I #6 Dec 16 2016 14:03:38
%S 1,4,8,27,59,188,384,427,2525,71429,80727,357492,13200877,65161876,
%T 7439912342,15555881542,71559279848,116275866868,345574982189,
%U 737460049244,9183275685671,12641946167319,126181443702371
%N Engel expansion of the Glaisher-Kinkelin constant A074962.
%C See A006784 for more details on the Engel expansion.
%e 1.2824271291006... = 1/1 + 1/(1*4) + 1/(1*4*8) + 1/(1*4*8*27) + ...
%t EngelExp[A_, n_] := Join[Array[1 &, Floor[A]], First@Transpose@NestList[{Ceiling[1/Expand[#[[1]] #[[2]] - 1]], Expand[#[[1]] #[[2]] - 1]} &, {Ceiling[1/(A - Floor[A])], A - Floor[A]}, n - 1]];
%t EngelExp[Glaisher, 22]
%Y Cf. A074962.
%K nonn
%O 1,2
%A _Benedict W. J. Irwin_, Dec 16 2016
|