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.