Assume $f$ is increasing and $\sum_n 1/f(n)<\infty$. I am trying to prove that $$f(n)\sum_{k=n}^{\infty}\frac{1}{f(k)}\le C n^2$$ for some positive constant $C$ (I believe that $C$ could be arbitrarily small and instead of $n^2$ in RHS we can replaced with $n^{1+\epsilon}$ ). This is equivalent to $$\limsup_{n\to\infty}n^{-2} f(n)\sum_{k=n}^{\infty}\frac{1}{f(k)}<\infty$$
One could bring this inequality into the inegral form $$ \limsup_{x\to\infty}x^{-2} f(x)\int_{x}^{\infty}\frac{dy}{f(y)}<\infty$$ I tried with L'Hopital rule but it does not really works.
Surprisingly, not only is this conjecture false, but it remains false if $Cn^2$ is replaced by any function $B(n)$ of $n$, no matter how fast it grows!
The idea is to use a sequence of values $(f(n))_{n=1}^\infty$ that is constant for very long stretches before switching to a new value. (Indeed, I imagine the OP tested the conjecture against standard smooth functions like powers of $n$, exponential functions, and things like $f(n) = n(\log n)^2$; in general, constant-for-long-stretches functions are a really good complement to "normal" functions for testing convergence-related conjectures.)
For simplicity let's assume that $B(n)$ is positive-integer-valued. Define a sequence of cutoff points $x_j$ and a sequence of values $v_j$ recursively as follows:
In other words, the sequence $\biggl(\dfrac1{f(n)}\biggr)_{n=1}^\infty$ consists of $B(x_1)$ copies of $v_1$, then $B(x_2)$ copies of $v_2$, and so on. In particular, $$ \sum_{n=1}^\infty \frac1{f(n)} = \sum_{j=1}^\infty v_j B(x_j) = \sum_{j=1}^\infty \frac1{2^j} = 1. $$ On the other hand, $$ f(x_j) \sum_{k=x_j}^\infty \frac1{f(k)} > f(x_j) \sum_{k=x_j}^{x_{j+1}-1} \frac1{f(k)} = \sum_{k=x_j}^{x_{j+1}-1} 1 = x_{j+1}-x_j = B(x_j); $$ in particular, $\displaystyle f(n) \sum_{k=n}^\infty \frac1{f(k)} > B(n)$ for infinitely many $n$.