If it were $e^z$ being expanded at $z=1$, it would be relatively easy. We would take advantage of the fact that $$e^z =\sum_{-\infty}^\infty a_n(z-1)^n$$ where $a_n$ is given by the closed-loop integral $$a_n = \frac{1}{2\pi i} \int_\gamma \frac{e^z}{(z-1)^{n+1}}dz$$ If $n+1\leq0$, then the function inside of the integral is completely defined, so $a_n=0$. So we only worry about when $n \geq 0$. For this, we take advantage of Cauchy's integral formula to get $$a_n = \frac{e}{n!}$$ And hence $$e^z = \sum_{n=0}^\infty e(z-1)^n$$ at $z=1$
I wanted to apply the same system to evaluating $e^{1/z}$ at $z=1$. But unfortunately, $\frac{d}{dz}e^{1/z} \neq e^{1/z}$, so the Cauchy integral formula will not be of use. Is there any trick for evaluating this one?
Let $0<r<1$ and $\gamma$ be the counterclockwise contour $|z-1|=r$. Then we have for $n\ge 1$ and $|z-1|<1$
$$\begin{align} a_n&=\frac1{2\pi i}\oint_{\gamma}\frac{e^{1/z}}{(z-1)^{n+1}}\,dz\\\\ &\overbrace{=}^{z=1/w}\frac{(-1)^n}{2\pi i}\oint_{\gamma'}\frac{w^{n-1}e^w}{(w-1)^{n+1}}\,dw\\\\ &=\frac{(-1)^n}{n!} \lim_{w\to 1}\frac{d^n(w^{n-1}e^w)}{dw^n}\\\\ &=\frac{(-1)^n}{n!} \lim_{w\to 1}\sum_{k=0}^n \binom{n}{k}\frac{d^{n-k}e^w}{dw^{n-k}}\frac{d^kw^{n-1}}{dw^k}\\\\ &=\frac{(-1)^n e}{n!}\sum_{k=0}^n \binom{n}{k}\frac{(n-1)!}{(n-1-k)!}\\\\ &=(-1)^n (n-1)! e \sum_{k=0}^{n-1}\frac{1}{k!(n-k)!(n-1-k)!} \end{align}$$