This is a generalization (and answer to) Writting $S = \sum_{k=0}^{\infty} \frac{1}{(r_1+k+1)(r_2+k+1)(r_3+k+1)}$ as a rational function of $r_1,r_2$ and $r_3$.
Let $S =\sum_{n=1}^{\infty} \dfrac1{\prod_{k=1}^m (n+a_k)} $ where $m \ge 2$ and the $a_k$ are positive integers.
Find an expression for $S$ involving only finite sums.
I get $S =-\sum_{k=1}^mc_kH_{a_k} $ where $c_k= \dfrac1{\prod_{j=1, j\ne k}^m (a_j-a_k)} $ and $H_n = \sum_{j=1}^n \dfrac1{j}$.
It is interesting that a quite good approximation to $S$ is $-\sum_{k=1}^mc_k\ln(a_k) $.
Provided that $\text{Re}(a),\text{Re}(b)>0$ we have $$ \sum_{n\geq 0}\frac{1}{(n+a)(n+b)}=\frac{\psi(a)-\psi(b)}{a-b},\tag{1} $$ hence by partial fraction decomposition, assuming $\text{Re}(c)>0$,
$$\begin{eqnarray*} \sum_{n\geq 0}\frac{1}{(n+a)(n+b)(n+c)}&=&\frac{1}{c-a}\left[\frac{\psi(a)-\psi(b)}{a-b}-\frac{\psi(b)-\psi(c)}{b-c}\right]\\&=&\frac{(c-b)\psi(a)+(a-c)\psi(b)+(a-b)\psi(c)}{(a-c)(a-b)(b-c)}\end{eqnarray*}\tag{2} $$ with $\psi(z)=\frac{d}{dz}\log\Gamma(z)=\frac{\Gamma'(z)}{\Gamma(z)}=H_{z-1}-\gamma$.
If more variables are involved the approach is just the same.