Since $(e_1, e_2) \cdot (x,y) = (x,y)$ should work for any $x,y \in \mathbb{R}$, we can pick some usefull values for them, like $0$ or $1$.
Let $(x,y) = (1,0)$. Then $(e_1, e_2) \cdot (1,0) = (e_1, 0)$ should be equal to $(1,0)$. So, $e_1 = 1$.
Let $(x,y) = (0,1)$. Then $(e_1, e_2) \cdot (0,1) = (e_2, 0)$ should be equal to $(0,1)$. So, $e_2 = 0$.
Now we can check that the pair $(1,0)$ works as an identity. In fact, it does: $(1, 0) \cdot (x,y) = (x\cdot 1 - 2 \cdot y \cdot 0, y \cdot 1 + x \cdot 0) = (x, y)$. Multiplication on the right can be checked similarily.
To make the whole construction more clear, you can notice that it resembles complex numbers multiplication, but with a "$2$" inside the multiplication formula. This is a consequence of the fact that this formula is exactly the multiplication in factor ring $\mathbb{R}[X]/(X^2 + 2)$ (while ordinary complex numbers are defined as $\mathbb{R}[X]/(X^2 + 1)$), that is, a factorization of the ring of polynomials over $\mathbb{R}$ by the ideal generated by $X^2+2$.