$f(x_1,y_1)=f(x_2,y_2) \implies x_1^2-y_1^2=x_2^2-y_2^2$ and $2x_1y_1=2x_2y_2$. Multiplying the first equation by $y_1^2y_2^2$,
$$
x_1^2y_1^2y_2^2-y_1^4y_2^2=x_2^2y_1^2y_2^2-y_2^4y_1^2
$$
Rearranging and using the second equation,
$$
x_1^2y_1^2(y_2^2-y_1^2)+y_1^2y_2^2(y_2^2-y_1^2)=0$$
So
$$
(y_2^2-y_1^2)(x_1^2y_1^2+y_1^2y_2^2)=(y_2^2-y_1^2)y_1^2(x_1^2+y_2^2)=0
$$
There are three cases; $y_1=\pm y_2$ leads to $x_1=\pm x_2$
$y_1=0$ leads to $x_2=x_1=y_2=0$ or $y_2=0 ,x_1=\pm x_2$
$x_1^2+y_1^2=0$ leads to $x_1=y_1=x_2=y_2=0$
Overall, summing them up we get $(x_1,y_1)=\pm(x_2,y_2)$
For surjectivity, if $f(x,y)=(a,b)$ then $x^2-y^2=a$ and $2xy=b$.
You can proceed as is suggested by the other answer; but you can sort of see it immediately via graphs
if $b=0$ then one of $x$ or $y$ is $0$, so that case is simple (depending on sign of $a$ you can choose either $x$ or $y$ to be $0$ and the other one would be $\sqrt{|a|}$)
If $b\neq 0$ the solution to system of equations is decribed by the intersection of the graphs of $y=\frac{b}{2x}$ and $x^2-y^2=a$; But it is quite immediate that the graphs must have an intersection.
Also note that the intersection points of the graphs are symmetric about the origin, and that's another way of seeing that $f$ is 2 to 1 function except at origin.