Find all functions $f$ from the positive integers to the positive integers such that such that for all integers $x, y$ we have $$2yf(f(x^2)+x)=f(x+1)f(2xy).$$