I'm stuck on this question. My thoughts are to try and first make each of these random variables have mean 0 by subtracting out $\lambda_i$. Then I think you can use some kind of borel cantelli argument by trying to show $P(|S_n-\sum_{i=1}^n\lambda_i|>\epsilon \text{ i.o })=0$ Is this the right track? Any hints or solutions would be appreciated
Source: Old Qual problem
https://dornsife.usc.edu/assets/sites/990/docs/Spring_2012/20121_507a.pdf

Yes, you want to do a Borel-Cantelli type argument. A hint for how to check the decay property that you need:
$$\sum_{n=N}^\infty \frac{\lambda^n}{n!} \leq \frac{e^\lambda \lambda^N}{N!}$$
which follows from the Lagrange remainder formula.