Take for instance $$ \frac{1}{\sqrt{z^2-6z+1}}=\sum_{\ell=0}^\infty c_\ell z^\ell\ . $$ The coefficient $c_\ell$ of the Maclaurin series are all positive integers $$ \vec{c}=\{1,3,13,63,321,\ldots\} $$ given by the formula $c_\ell=p_\ell(3)$, where $p_\ell(t)$ is a Legendre polynomial.
I am wondering if there exist necessary/sufficient conditions for a function $f(z)$ to ensure that its Maclaurin expansion has (positive) integer coefficients. I have never heard of anything like that, but it seems like a very curious property for a function.
Many thanks for your help.