Integrate[Gamma[z + 1] z^-z, {z, 0, \[Infinity]}]