Consider a $X_1, ... X_n \sim~ Poisson(\lambda)$, I want to obtain the UMVUE of $P_{\lambda} (X=r)$.
This is my approach: $\operatorname{\mathbb{E}}_{\theta}[h(t)] = P_{\lambda} (X=r)$. The probability is dependent on $\lambda$ so I am inclined to find the conditional probability and proceed as $\operatorname{\mathbb{E}}_{\theta}[h(t)] = P(X=r, \lambda)/ P(\lambda)$. I am not sure if this is necessary but would appreciate some advice.
Incorporating provided suggestions I get: $$\sum_{t=0}^{\infty} h(t) \lambda^t\exp(-t)/ t! = \lambda^r\exp(-r)/ r!$$ I will expand the RHS: $\lambda^r\exp(-r)/ r! = \lambda^r\sum_{n=0}^{\infty} (-r)^{n}/ n!r!$
But I am stuck.