Calculation of a strange series Is it possible to find an expression for:
$$S(N)=\sum_{k=0}^{+\infty}\frac{1}{\sum_{n=0}^{N}k^n}?$$
For $N=1$ we have
$$S(1) = \displaystyle\sum_{k=0}^{+\infty}\frac{1}{1 + k} = \displaystyle\sum_{k=1}^{+\infty}\frac{1}{k}$$
which is the (divergent) harmonic series. Thus, $S (1) = \infty$.
For $N=2$ this sum is:
$$S(2)=\sum_{k=0}^{+\infty}\frac{1}{1+k+k^2}$$
which can be expressed as:
$$S(2)=-1+\frac{1}{3}\sqrt 3 \pi \tanh(\frac{1}{2}\pi\sqrt 3)\approx 0.798$$
For $N=3$ we have:
$$S(3)=\frac{1}{4}\Psi(I)+\frac{1}{4I}\Psi(I)-\frac{1}{4I}\pi\coth(\pi)+\frac{1}{4}\pi\coth(\pi)+\frac{1/}{4}\Psi(1+I)-\frac{1}{4I}\Psi(1+I)-\frac{1}{2}+\frac{1}{2}\gamma \approx 0.374$$

Let $T(N) = S(N-1)$. Then
$$ \begin{align*}T(n) 
&= 1 + \frac{1}{n} + \sum_{k=2}^{\infty} \frac{1}{k^{n-1}+k^{n-2}+\cdots+k+1} \\
&= 1 + \frac{1}{n} + \sum_{k=2}^{\infty} \frac{k - 1}{k^n - 1} \\
&= 1 + \frac{1}{n} + \sum_{k=2}^{\infty} \frac{1}{n} \sum_{l=1}^{n-1} \frac{\omega_l (\omega_l - 1)}{k - \omega_l}  \\
&= 1 + \frac{1}{n} + \sum_{k=0}^{\infty} \frac{1}{n} \sum_{l=1}^{n-1} \frac{\omega_l (\omega_l - 1)}{k + 2 - \omega_l},
\end{align*}$$
where $\omega_l = \exp\left(\tfrac{2\pi l i}{n}\right)$. Since 
$$ \frac{1}{n} \sum_{l=0}^{n-1} \omega_l (\omega_l - 1) = 0, $$
we may write
$$ \begin{align*}T(n) 
&= 1 + \frac{1}{n} + \sum_{k=0}^{\infty} \frac{1}{n} \sum_{l=1}^{n-1} \omega_l (\omega_l - 1) \left( \frac{1}{k + 2 - \omega_l} - \frac{1}{k+1} \right) \\
&= 1 + \frac{1}{n} + \frac{1}{n} \sum_{l=1}^{n-1} \omega_l (\omega_l - 1) \sum_{k=0}^{\infty} \left( \frac{1}{k + 2 - \omega_l} - \frac{1}{k+1} \right) \\
&= 1 + \frac{1}{n} - \frac{1}{n} \sum_{l=1}^{n-1} \omega_l (\omega_l - 1) \left( \gamma + \psi_0 (2 - \omega_l) \right) \\
&= 1 + \frac{1}{n} - \frac{1}{n} \sum_{l=1}^{n-1} \omega_l (\omega_l - 1) \psi_0 (2 - \omega_l).
\end{align*}$$
