0,3

Define c(n)=A067764(n)/A067653(n). For a given sequence s(n) consider P[s(n)](z):=e^(-z/(1-z))*Sum_{k>=0} s(k)c(k)z^k. Regarding complex valued abelian limitation the following holds true: if s(n) is convergent (to the limit s) then lim P[s(n)](z)=s as z tends to +1 in a certain subdomain D of the unit circle. There are two constraints: (1) D contains the line [0,1[. (2) There is a d>0 such that the intersection of {w|Re(w)>1-d} and D is a nonempty subset of a generalized Stolz set defined by {w||Im(w)|<=t*(1-Re(w))^(3/2)}, t<1. If z tends to +1 from outside such a domain that limit doesn't exist in general. - Hieronymus Fischer, Oct 20 2010

There is a significant overlap between the terms of this sequence and the terms of factorials (A000142): 44 of the first 100 terms of this sequence are also factorials. - Harvey P. Dale, Oct 27 2011

If there is no cancellation in the defining sum, a(n) = n!. In particular, a(p) = p! for prime p. In any case a(n) | n!. - Charles R Greathouse IV, Oct 27 2011

H. Fischer, Eine Theorie komplexwertiger Abelscher Limitierungsmethoden (A theory of complex valued abelian limitation methods), Dissertation (1987), pp. 29-32.

K. Knopp, Theory and application of infinite series, Dover, p. 547.

O. Perron, Über das infinitäre Verhalten der Koeffizienten einer gewissen Potenzreihe, Archiv d. Math. u. Phys. (3), Vol. 22, pp. 329-340, 1914.

K. Zeller, W. Beekmann, Theorie der Limitierungsverfahren, Springer-Verlag, Berlin (1970).

Alois P. Heinz, Table of n, a(n) for n = 0..450

D. Borwein, On methods of summability based on power series, Proc. Royal Soc. Edinburgh, Sect. A Volume 64 / Issue 04 / January 1957, pp 342-349.

Richard P. Brent, M. L. Glasser, and Anthony J. Guttmann, A Conjectured Integer Sequence Arising From the Exponential Integral, arXiv:1812.00316 [math.NT], 2018.

a(n) is the denominator of Sum_{i=1..n} binomial(n-1, i-1)/i!.

exp(x/(1-x)) = 1 + x + (3/2)*x^2 + (13/6)*x^3 + (73/24)*x^4 + (167/40)*x^5 + (4051/720)*x^6 + (37633/5040)*x^7 + (43817/4480)*x^8 + (4596553/362880)*x^9 + ... .

b:= proc(n) option remember; `if`(n=0, 1,

add((n-k)*b(k), k=0..n-1)/n)

end:

a:= n-> denom(b(n)):

seq(a(n), n=0..25); # Alois P. Heinz, May 12 2016

Denominator[Rest[CoefficientList[Series[Exp[x/(1-x)], {x, 0, 20}], x]]] (* Harvey P. Dale, Oct 26 2011 *)

r[n_] := If[n == 0, 1, Hypergeometric1F1[1 - n, 2, -1]]; Table[Denominator@ r[n], {n, 0, 22}] (* Peter Luschny, Feb 02 2019 *)

(PARI) apply(x->denominator(x), Vec(exp(x/(1-x)))) \\ Charles R Greathouse IV, Oct 27 2011

(Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(Exp(x/(1-x)))); [Denominators(b[n]): n in [1..m]]; // G. C. Greubel, Dec 04 2018

(SageMath) [denominator(sum(binomial(n-1, j-1)/factorial(j) for j in (1..n))) for n in range(30)] # G. C. Greubel, Dec 04 2018

Cf. A067764.

Sequence in context: A257546 A274038 A143383 * A090755 A192196 A000496

Adjacent sequences: A067650 A067651 A067652 * A067654 A067655 A067656

nonn,frac

Benoit Cloitre, Feb 03 2002

a(0)=1 prepended by Alois P. Heinz, May 12 2016

approved