Let $X_1,..., X_n$ be iid sample from the Poisson distribution with parameter $\lambda$. Find the UMVUE of $\lambda + \lambda^2$.
I know $T := \sum\limits_{i=1}^n X_i$ is complete and sufficient for $\lambda$. Also, $T/n$ is the MLE of $\lambda$. By the invariance principle of the MLE, $T/n + T^2 / n^2$ is the MLE for $\lambda + \lambda^2$.
I tried following the approach for computing UMVUE of $\lambda^3$ from an earlier post, but could not compute the expectation that gives $\lambda + \lambda^2$.