I hae to prove that $$\int_0^\infty\frac{x^n}{1+e^{x-t}}\mathrm{d}x = \frac{t^{n+1}}{n+1} + o(t^n), \quad\text{ when } t \to \infty,\,n\in\Bbb{R}^+$$ where $o(\cdot)$ is the Little-o notation.
What I have done so far:
Let $F_n(t) = \int_0^\infty\frac{x^n}{1+e^{x-t}}\mathrm{d}x$, then $F_0(t) =\int_0^\infty\frac{1}{1+e^{x-t}}\mathrm{d}x = \log(e^t+1)$. The derivative of $F_n(t)$ is: \begin{align*} F_n'(t) &= \int_0^\infty\frac{e^{x-t}x^n}{(1+e^{x-t})^2}\mathrm{d}x\\ &= -\int_0^\infty x^n\mathrm{d}\frac{1}{1+e^{x-t}}\\ &= \left.\frac{x^n}{1+e^{x-t}}\right|_0^\infty+\int_0^\infty \frac{1}{1+e^{x-t}}\mathrm{d}x^n\\ &= \int_0^\infty \frac{nx^{n-1}}{1+e^{x-t}}\mathrm{d}x\\ &= nF_{n-1}(t) \end{align*}
First, let's split up the integral: $$ \begin{align} \int_0^\infty\frac{x^n}{1+e^{x-t}}\mathrm{d}x &=\int_{-t}^\infty\frac{(x+t)^n}{1+e^x}\mathrm{d}x\\ &=\color{#C00000}{\int_{-t}^0\frac{(x+t)^n}{1+e^x}\mathrm{d}x}+\color{#00A000}{\int_0^\infty\frac{(x+t)^n}{1+e^x}\mathrm{d}x}\tag{1} \end{align} $$ Note that the first integral on the right side of $(1)$ is $$ \begin{align} \color{#C00000}{\int_{-t}^0\frac{(x+t)^n}{1+e^x}\mathrm{d}x} &=t^{n+1}\int_{-1}^0\frac{(1+x)^n}{1+e^{tx}}\mathrm{d}x\\ &=t^{n+1}\int_{-1}^0(1+x)^n\,\mathrm{d}x-\color{#0000FF}{t^{n+1}\int_{-1}^0\frac{(1+x)^n}{1+e^{tx}}e^{tx}\,\mathrm{d}x}\\ &=\frac{t^{n+1}}{n+1}+O\left(t^n\right)\tag{2} \end{align} $$ because $$ \begin{align} \color{#0000FF}{t^{n+1}\int_{-1}^0\frac{(1+x)^n}{1+e^{tx}}e^{tx}\,\mathrm{d}x} &\le t^{n+1}\int_{-1}^0(1+x)^n\,e^{tx}\mathrm{d}x\\ &\le t^{n+1}\int_{-1}^0e^{nx}\,e^{tx}\mathrm{d}x\\ &\le\frac{t^{n+1}}{n+t}\\[6pt] &=O\left(t^n\right)\tag{3} \end{align} $$ Furthermore, by dominated convergence $$ \lim_{t\to\infty}\int_0^\infty\left(1+\frac xt\right)^n\,e^{-x}\,\mathrm{d}x =1\tag{4} $$ therefore, the second integral on the right side of $(1)$ is $$ \begin{align} \color{#00A000}{\int_0^\infty\frac{(x+t)^n}{1+e^x}\mathrm{d}x} &\le t^n\int_0^\infty\left(1+\frac xt\right)^n\,e^{-x}\,\mathrm{d}x\\ &=O\left(t^n\right)\tag{5} \end{align} $$ Combining $(1)$, $(2)$, and $(5)$, we get $$ \int_0^\infty\frac{x^n}{1+e^{x-t}}\mathrm{d}x =\frac{t^{n+1}}{n+1}+O\left(t^n\right)\tag{6} $$