According to p.244 in "Magnus, W., Oberhettinger, F., Soni, R.: Formulas and Theorems for the Special Functions of Mathematical Physics, Grundlehren der mathematischen Wissenschaften 52, Springer, Berlin et al. (1966)", we have, $$\int_{0}^\infty e^{-zx}\, x^\alpha\, L^{\alpha}_m(x)\, dx = \frac{\Gamma(m+\alpha+1)}{\Gamma(m+1)} z^{-m-\alpha-1}\,(z-1)^{m}; \quad \text{Re}\, z> 0,\, \text{Re}\, \alpha> 0.$$ I am looking for the value of the integral $$I:= \int_{0}^\infty (x+1)^p\, e^{-ax}\, x^\alpha\, L^{\alpha}_m(x)\, dx =?; \quad a> 0,\,p\in\mathbb Z.$$
The Laguerre polynomial of type $\alpha$ and degree $m$ is given by
$$L^{\alpha}_m(x)=\sum_{k=0}^{m}(-1)^k \binom{m+\alpha}{m-k} \frac{x^k}{k!}.$$
Then, by this definition, we get
$$I:= \sum_{k=0}^{m} \frac{(-1)^k}{k!} \binom{m+\alpha}{m-k} \int_{0}^\infty (x+1)^p\, e^{-ax}\, x^{\alpha+k}\, dx$$
So to calculate I, it is enough to calculate the following integral:
$$J:=\int_{0}^\infty (x+1)^p\, e^{-ax}\, x^{\alpha+k}\, dx=?$$
According to p.277 in "Magnus, W., Oberhettinger, F., Soni, R.: Formulas and Theorems for the Special Functions of Mathematical Physics, Grundlehren der mathematischen Wissenschaften 52, Springer, Berlin et al. (1966)", we have,
$$U(\gamma,c,a)=\frac{1}{\Gamma(\gamma)}\int_{0}^\infty (1+x)^{c-\gamma-1}\, e^{-ax}\, x^{\gamma-1}\, dx; \quad \text{Re}\, \gamma> 0,\, \text{Re}\,a> 0.$$
Then for $\gamma=\alpha+k+1$ and $c=\alpha+p+k+2$, we get
$$J:=\Gamma(\alpha+k+1) U(\alpha+k+1,\alpha+p+k+2,a).$$
Then,
$$I:= \sum_{k=0}^{m} \frac{(-1)^k}{k!} \binom{m+\alpha}{m-k} \Gamma(\alpha+k+1) U(\alpha+k+1,\alpha+p+k+2,a).$$