Let $n \in \mathbb{N}$ and $a_1, \dots, a_n, b_1, \dots, b_n > 0$. Define the function $$ P(x) = \Bigg(\sum_{i=1}^n \frac{a_i}{x-b_i}\Bigg) - 1, \qquad x \in \mathbb{R}_+ \setminus \{b_1, \dots, b_n\}. $$
How can I check that $P$ has $n$ distinct positive roots? I have rearranged $P(x) = 0$ as $$ \sum_{i=1} a_i \prod_{j \neq i} (x-b_j) = \prod_{i=1}^n (x-b_i), $$ but did not know how to proceed. I have also looked at Sturm chains without much success.