Deje $P(x) = \prod\limits_{k=1}^n( a_k + b_k x)$ donde $n > 1$, $a_k, b_k > 0$ y las raíces $-\frac{a_k}{b_k}$ son distintos.
Considere la integral
$$\int_0^\infty \frac{dx}{P(x)} = \lim_{\Lambda\to\infty}\int_0^\Lambda \frac{dx}{P(x)}\tag{*1}$$
Puesto que todas las raíces de $P(x)$ son distintos, tenemos las siguientes parcial fracción de descomposición
$$\frac{1}{P(x)} = \sum_{j=1}^n \frac{1}{a_j + b_j x}\prod_{k=1,\ne j}^n \frac{1}{a_k - b_k\frac{a_j}{b_j}}
= \sum_{j=1}^n \frac{1}{a_j + b_j x}\prod_{k=1,\ne j}^n \frac{b_j}{a_k b_j - b_k a_j}\etiqueta{*2}$$
Aviso
$$\int_0^\Lambda \frac{dx}{a_j + b_j x} = \frac{1}{b_j}\left(\log(a_j + b_j \Lambda) - \log(a_j)\right)
= \frac{1}{b_j}\left[\log\Lambda + \log\left(\frac{b_j}{a_j} + \frac{1}{\Lambda}\right)\right]\etiqueta{*3}$$
Cuando se conecta $(*2)$ en RHS de $(*1)$ e integrar, $(*3)$ nos dice
hay un término proporcional a $\log\Lambda$. El coeficiente del término
es igual a
$$
\sum_{j=1}^n \frac{1}{b_j}\prod_{k=1,\ne j}^n \frac{b_j}{a_k b_j - b_k a_j}
= \lim_{x\to\infty}
\sum_{j=1}^n \frac{x}{a_j+b_jx}\prod_{k=1,\ne j}^n \frac{b_j}{a_k b_j - b_k a_j}
= \lim_{x\to\infty} \frac{x}{P(x)} = 0$$
Aviso por el plazo restante $(*3)$ converge a$\frac{1}{b_j}\log\left(\frac{b_j}{a_j}\right)$$\Lambda \to \infty$, obtenemos:
$$\int_0^\infty \frac{dx}{P(x)} =
\sum_{j=1}^n
\log\left(\frac{b_j}{a_j}\right)b_j^{n-2}
\prod_{k=1,\ne j}^n \frac{1}{a_k b_j - b_k a_j}$$
Para $(a_1,b_1,a_2,b_2,a_3,b_3) = (a,b,c,d,1,p)$, la integral se convierte en
$$\frac{p\,\mathrm{log}\left( p\right) }{\left( a\,p-b\right) \,\left( c\,p-d\right) }+\frac{d\,\mathrm{log}\left( \frac{d}{c}\right) }{\left( a\,d-b\,c\right) \,\left( d-c\,p\right) }+\frac{b\,\mathrm{log}\left( \frac{b}{a}\right) }{\left( b\,c-a\,d\right) \,\left( b-a\,p\right) }$$
Con la ayuda de un CAS, he verificado que este coincide con la respuesta que recibe de Mathematica.