この方法は前にどこかのサイトで見かけたものです。
$\beginend{align}{ n &\ge 1 \\ S_n \acoloneqq \sum_{k=1}^n \frac{1}{p_k} \\ T_n \acoloneqq \prod_{k=1}^n \lr({1+\frac{1}{p_k}}) }$
$\displaystyle T_n < \prod_{k=1}^n \lr({\frac1{0!}+\frac{1}{1!p_k}+\frac1{2!p_k^2}+\cdots}) = \prod_{k=1}^n e^{\frac1{p_k}} = e^{S_n}$
$\displaystyle \prod_p \lr({1+\frac1p}) =
\prod_p \frac{1-\frac1{p^2}}{1-\frac1p} =
\prod_p \frac{\frac11+\frac1p+\frac1{p^2}+\cdots}{\frac11+\frac1{p^2}+\frac1{p^4}+\cdots} =
\frac{\frac11+\frac12+\frac13+\cdots}{\frac1{1^2}+\frac1{2^2}+\frac1{3^2}+\cdots} = \infty$
である為$T_n$は正の無限大に発散する。
よって$S_n$は正の無限大に発散する。