How to solve the given integral of the type: $\int_{0}^{\infty}\tfrac{1}{(a + b x)(c + d x) (1 + p x)}dx$

Let $P(x) = \prod\limits_{k=1}^n( a_k + b_k x)$ where $n > 1$, $a_k, b_k > 0$ and the roots $-\frac{a_k}{b_k}$ are distinct.
Consider the integral

$$\int_0^\infty \frac{dx}{P(x)} = \lim_{\Lambda\to\infty}\int_0^\Lambda \frac{dx}{P(x)}\tag{*1}$$

Since all the roots of $P(x)$ are distinct, we have following partial fraction decomposition

$$\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}\tag{*2}$$

Notice $$\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]\tag{*3}$$ When you plug $(*2)$ into RHS of $(*1)$ and integrate, $(*3)$ tells us there is a term proportional to $\log\Lambda$. The coefficient of that term equals to $$ \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$$ Notice the remaining term in $(*3)$ converges to $\frac{1}{b_j}\log\left(\frac{b_j}{a_j}\right)$ as $\Lambda \to \infty$, we obtain:

$$\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}$$

For $(a_1,b_1,a_2,b_2,a_3,b_3) = (a,b,c,d,1,p)$, the integral becomes $$\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) }$$

With help of a CAS, I have verified this matches the answer you get from Mathematica.