Алгоритм для одной суммы обратных величин

Описан алгоритм для приближенного вычисления набора сумм вида
$$ \varphi_k=\sum_{j=1}^nc_j/(\lambda_j+\lambda_k),\;1\leqslant k\leqslant n, \text{ где }0<\alpha\leqslant \lambda_j\leqslant \beta. $$
Время работы алгоритма $O(n(t+\log n)\Psi(t+\log n))$, если $\varphi_k$ вычисляется с точностью $2^{-t}$; здесь функция $\Psi(l)$ обозначает время умножения двух целых чисел двоичной длины $l$.