OFFSET
1,2
FORMULA
Denominators of coefficients in expansion of (Ei(x) - log(x) - gamma) / (1 - x), x > 0.
EXAMPLE
1, 5/4, 47/36, 379/288, 9487/7200, 14233/10800, 87179/66150, ...
MATHEMATICA
Table[Sum[1/(k k!), {k, 1, n}], {n, 1, 20}] // Denominator
nmax = 20; Assuming[x > 0, CoefficientList[Series[(ExpIntegralEi[x] - Log[x] - EulerGamma)/(1 - x), {x, 0, nmax}], x]] // Denominator // Rest
PROG
(PARI) a(n) = denominator(sum(k=1, n, 1/(k*k!))); \\ Michel Marcus, May 26 2022
(Python)
from math import factorial
from fractions import Fraction
def A354401(n): return sum(Fraction(1, k*factorial(k)) for k in range(1, n+1)).denominator # Chai Wah Wu, May 27 2022
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Ilya Gutkovskiy, May 25 2022
STATUS
approved