En primer lugar, se utiliza la definición de Euler-Mascheroni número
$$\gamma =\lim_{n\to \infty}\left(\sum_{k=1}^n\frac1k-\log n\right) \tag 1$$
A continuación, vamos a escribir la suma de los intereses, como
$$\begin{align}
\sum_{k=0}^{n}\frac{1}{pk+1}&=1+\frac1p \left(\sum_{k=1}^n\left(\frac{1}{k+1/p}-\frac1k\right)+\sum_{k=1}^n\frac1k \right)\\\\
&=1+\frac1p \sum_{k=1}^n\left(\frac{1}{k+1/p}-\frac1k\right)+\frac1p\left(\sum_{k=1}^n\frac1k -\log n\right)+\frac1p\log n\\\\
\sum_{k=0}^{n}\frac{1}{pk+1}-\frac1p\log n&=1+\frac1p \sum_{k=1}^n\left(\frac{1}{k+1/p}-\frac1k\right)+\frac1p\left(\sum_{k=1}^n\frac1k -\log n\right)\\\\
&=1+\frac1p \sum_{k=1}^n\left(\int_0^1x^{k+1/p-1}\,dx-\int_0^1x^{k-1}\,dx\right)\\\\
&+\frac1p\left(\sum_{k=1}^n\frac1k -\log n\right)\\\\
&=1+\frac1p \left(\int_0^1\frac{x^{1/p}-1}{1-x}\,dx\right)+\frac1p\left(\sum_{k=1}^n\frac1k -\log n\right) \tag 2\\\\
\end{align}$$
Para $p=3$, podemos evaluar la integral en $(2)$ en forma cerrada. De continuar, hemos de hacer cumplir la sustitución de $x\to x^3$. A continuación,
$$\begin{align}
\int_0^1\frac{x^{1/3}-1}{1-x}\,dx&=3\int_0^1\frac{(x-1)x^2}{1-x^3}\,dx\\\\
&=-3\int_0^1\frac{x^2}{x^2+x+1}\,dx\\\\
&=-3+\frac{\pi\sqrt{3}}{6}+\frac32 \log 3 \tag 3
\end{align}$$
El uso de $(3)$ $(2)$ $p=3$ y tomando el límite cuando $n\to \infty$ rendimientos
$$\bbox[5px,border:2px solid #C0A000]{\lim_{n\to \infty}\left(\sum_{k=0}^{n}\frac{1}{3k+1}-\frac13\log n\right)=\frac 13\gamma +\frac{\pi\sqrt{3}}{18}+\frac12 \log 3}$$
como iba a ser mostrado.
Para los valores de $p\in N$ diferente de la $3$, uno todavía puede llevar a cabo la integral en $(2)$ y llegar al límite
$$\bbox[5px,border:2px solid #C0A000]{\lim_{n\to \infty}\left(\sum_{k=0}^{n}\frac{1}{pk+1}-\frac1p\log n\right)=\frac 1p\gamma +\int_0^1\frac{1-x^{p-1}}{1-x^p}\,dx} \tag 4$$
NOTA:
El uso de Gauss, digamma teorema, podemos evaluar el lado derecho de la $(4)$ en forma cerrada. Obtenemos
$$\begin{align}
\lim_{n\to \infty}\left(\sum_{k=0}^{n}\frac{1}{pk+1}-\frac1p\log n\right)=\frac 1p\gamma +\frac1p\left(\log (2p)+\frac{\pi}{2}\cot\left(\frac{\pi}{p}\right)-2\sum_{k=1}^{\lfloor\frac{p-1}{2}\rfloor}\cos\left(\frac{2\pi k}{p}\right)\log\left(\sin\left(\frac{\pi k}{p}\right)\right)\right)
\end{align}$$