Giải hệ phương trình \[ \begin{cases} 2z(x+y)+1=x^2-y^2 \\ y^2+z^2=1+2xy+2xz-2yz \\ y(3x^2-1)=-2x(x^2+1) \end{cases} \]