Applying additivity over equal sums twice, f(x+y)+f(1)=f(x)+f(y+1) and f(y+1)+f(1)=f(y)+f(2); combining gives f(x+y)=f(x)+f(y)+b for all x,y>0, where b=f(2)−2f(1)=f(2)−2. Using yf(x)=f(xf(y)), the involution f(f(y))=y, and this additive relation at x=y=2, one computes 4+2b=2f(2)=f(2f(2))=f(f(2)+f(2))=f(f(2))+f(f(2))+b=4+b, forcing b=0.