$x+iy$ forma parte del ideal generado por a+ib si se puede escribir como $(ac-bd)+i(cb+ad) = (a+ib)(c+id)$ .
Así que $\begin{align} x&= ac-bd\\ y&= bc+ad \end{align} $
Resolver para $c$ y $d$ encontramos
$\begin{align} c&= (x+bd)/a\\ d&=(ay-bx)/(a^2 + b^2) \end{align}$
Desde $d$ tiene que ser un número entero tenemos que
$a^2 + b^2 | ay-bx$
Por la identidad de Bézout sabemos que $ay-bx$ puede ser cualquier número entero si $a$ y $b$ son coprimas, por lo que tenemos $a^2 + b^2$ clases de equivalencia.
Si $a$ y $b$ no son coprimas (sea $d$ sea el GCD) sólo tenemos $(a^2 + b^2)/d$ clases de equivalencia.
Nunca he utilizado un editor de texto matemático, así que lo siento por la notación. Me siento raro por no usar el hecho de que $a|(x+bd)$ pero el resto debería estar bien. Me gustaría que alguien me explicara cómo se escriben las matemáticas aquí. Gracias.