De hecho hay una relación con los números de Bernoulli de segunda clase. Si estoy en lo correcto, este es el resultado :
Teorema : la generación de la función de f_n(x)
\sum_{n=0}^{+\infty} f_n(x)z^n = \dfrac{(1+z)^{x-1}-1}{\log (1+z)}
Corrolary :
f_n(x) = \sum_{k=0}^n {{x-1}\choose{k+1}} \dfrac{b_{n-k}}{(n-k)!}
Voy a presentar aquí principalmente a la "formal" pasos para llegar a este resultado, voy a omitir los detalles de la convergencia, la derivación bajo integral y de todos. Inicio con :
f_n(x)=\dfrac{1}{n!} \int_1^x \frac{\Gamma(t)}{\Gamma(t-n)} \text{d}t=\dfrac{1}{n!} \int_1^x (t-1)(t-2) \cdots (t-n) \text{d}t
Ahora introduce g(t,y)=y^{t-1}. Tenemos \dfrac{\partial^n g}{\partial y^n}(t,y)=(t-1)(t-2) \cdots (t-n)y^{t-n-1}. Por lo tanto :
\begin{array}{rcl}
f_n(x,y) & = & \displaystyle \frac{1}{n!} \int_1^x (t-1)(t-2) \cdots (t-n)y^{t-n-1} \text{d}t\\
& = & \displaystyle \frac{1}{n!} \dfrac{\partial^n \;}{\partial y^n} \left( \int_1^x y^{t-1} \text{d}t\right) \\
& = & \displaystyle \frac{1}{n!} \dfrac{\partial^n \;}{\partial y^n} \left( \frac{y^{x-1}-1}{\log(y)}\right) \\
& = &
\end{array}
Tenemos f_n(x)=f_n(x,1). Por lo tanto, obtenemos (modulo algún buen argumento para la expansión de Taylor) :
\boxed{\sum_{n=0}^{+\infty} f_n(x)z^n = \sum_{n=0}^{+\infty} \displaystyle \left. \frac{1}{n!} \dfrac{\partial^n \;}{\partial y^n} \left( \frac{y^{x-1}-1}{\log(y)}\right) \right|_{y=1} z^n = \dfrac{(1+z)^{x-1}-1}{\log (1+z)}}
Esto demuestra el teorema. El corolario es una consecuencia directa de la siguiente serie exapnsion :
\displaystyle \dfrac{(1+z)^{x-1}-1}{z} = \sum_{n=0}^{+\infty} {{x-1}\choose{n+1}}z^n \; \; \; \; \; \; \; \text{and} \; \; \; \; \; \; \; \displaystyle \frac{z}{\log(1+z)} = \sum_{n=0}^{+\infty} \frac{b_n}{n!}z^n
f_n(x) es el dado por el producto de Cauchy :
\boxed{f_n(x) = \sum_{k=0}^n {{x-1}\choose{k+1}} \dfrac{b_{n-k}}{(n-k)!}}
No estoy seguro de que podemos hacer más que esto.
Por curiosidad, ¿cómo has encontrado con este problema ?