La forma "madura" de pensar en esto es como una función en la esfera de Riemann - este es el plano complejo estándar con un punto $\infty$ también llamado "plano complejo extendido". La función $1/z$ es un caso especial de lo que se denomina transformación fraccionaria lineal - todas ellas son funciones de la forma $\frac{az+b}{cz+d}$ donde los coeficientes están en $\mathbb{C}$ . Una propiedad de las transformaciones lineales fraccionarias es que envían "círculos generalizados" a "círculos generalizados", donde un círculo generalizado es un círculo normal o una recta (una recta se llama círculo generalizado porque se puede ver como un círculo que pasa por el punto en el infinito).
En este caso, vemos que nuestro círculo a través de 0 será enviado a través de un "círculo a través del infinito" porque definimos $1/0$ ser $\infty$ aquí. Así que sus cuatro puntos serán enviados a una línea, y en particular, $1/z_1$ , $1/z_2$ y $1/z_3$ serán colineales.
Las transformaciones lineales fraccionarias y sus propiedades son un tema estándar en el análisis complejo, y ni siquiera hace falta saber análisis para demostrar la propiedad que acabo de mencionar. Cualquier texto estándar sobre análisis complejo debería tener la demostración. Te sugiero que busques una copia del libro de Gamelin, si te interesa.