Is this $f(x) = x+1$ the only solution to this functional equation.

I am considering the problem of finding all functions $f:(0,\infty)\to(0,\infty)$ satisfying the functional equation:

$$f\big(xf(y)+f(x)\big) = 2f(x)+xy\text.$$

I have been able to prove the following three results/properties:

1. $f$ is not surjective.
2. $f$ does not have any fixed points.
3. $f(x)=x+1$ is a solution.

My intuition tells me that $x\mapsto x+1$ is the only solution, but I have not been successful in proving or disproving this claim.

Any ideas on how I can make further progress is appreciated.

As suggested in comments, here is a solution by user pco in https://artofproblemsolving.com/community/c6h611705p3637387:

Let $P(x,y)$ be the assertion $f\big(xf(y)+f(x)\big)=2f(x)+xy$. Let $a=f(1)$.

$P(1,x) \implies f\big(f(x)+a\big)=x+2a$ and so $(2a, +\infty) \subseteq f(\mathbb{R})$

$P\big(1,f(x)+a\big) \implies f(x+3a)=f(x)+3a$ and so $f(x+3ka)=f(x)+3ka$ $\forall x >0, \forall k \in \mathbb{Z}_{\ge 0}$.

Let $x,y>0$ and $n \in \mathbb{N}$ such that $y+3na>2ax+f(x)$. Then there is $t>0$ such that $y+3na=xf(t)+f(x)$. Comparing $P(x,t)$ with $P(x,t+3a)$, we get $f(y+3na+3ax)=f(y+3na)+3ax$ and so $f(y+3ax)=f(y)+3ax$. So $f(x+y)=x+f(y)$ $\forall x,y >0$ and so $f(x)=x+c$ $\forall x$ and for some $c \in \mathbb{R}$.

Plugging this back in original equation, we get $c=1$ and so the unique solution is $f(x)=x+1$ $\forall x$.