# Sufficiency for Truncated Geometric

Here is a deviant of a question I feel like I have seen several times on truncated exponentials and similar distributions for finding sufficient statistics:

Let $$\mathbb{P}(Y=y)=\theta^y(1-\theta)^{\mathbb{I}\{y\leq z\}}$$ for some $$z\in\mathbb{N}$$ and $$y\in\mathbb{N}$$. Is the MLE for $$\hat{\theta}$$, $$\frac{n\bar{Y}-n}{n-\sum\limits_{i=1}^n \mathbb{I}\{y=z+1\}}$$, sufficient? If not, what else aside from the entire data $$Y=(y_1,\dots,y_n)$$, would constitute a sufficient statistic?

My Attempt: I believe that $$T(Y)=\left(\sum y_i,\sum \mathbb{I}\{y=z+1\})\right)$$ would easily work as a sufficient statistic. The MLE is a function (not 1-1 though) of this statistic which checks the box. As for if the MLE itself is sufficient (which intuitively I don't think is), we have

$$\mathcal{L}(\theta)=\theta^{\sum y_i}(1-\theta)^{n-\sum \mathbb{I}\{y_i=z+1\}},$$

but by inspection, I really see no way of (successfuly) using the factorization theorem here nor anything dealing with the exponential family as this does not reside in that family of distributions.

Question: Is there a more rigorous or correct way of going about this? I am being very hand-wavy here and was wondering if there's an alternative?