Creo que el argumento que me convence es éste:
Supongamos que \(r\equiv s \mod(n)\), \(s\equiv r\mod(n)\).
O sea, son "asociados" según tu terminología (¿esa es la terminología usual?).
Tenemos que existen \(\alpha,\beta \in\mathbb Z_n\) tales que
\(s\equiv \alpha r\),
\(r\equiv\beta s\),
de donde
\(r\equiv \gamma r\), donde \(\gamma = \alpha\beta\).
Esa ecuación, con \(\gamma\) como incógnita, tiene solución, obviamente,
y todas las soluciones \(\gamma\) se pueden encontrar mediante
\(\gamma = \xi + j\, n / \rho\),
donde \(\rho = \gcd(n,r)\),
y \(\xi\) es la única solución a la ecuación:
\[ \dfrac r\rho \xi \equiv \dfrac r\rho \mod(\nu); \qquad n/\gcd(n,r).\]
Notamos que hay solución única porque
\(R=\dfrac r\rho\) es coprimo con
\(N=n/\gcd(n,r)=n/\rho\).
Y entonces esa única solución no puede ser otra que \(\xi=1\).
Luego, las soluciones son
\(\gamma = \gamma_j = 1 + j \, N\),
para \(0\leq j<\rho\).
Si un primo \(p\) fuera divisor común de \(n\) y de \(\gamma\),
entonces no puede ser divisor de \(N\),
porque en tal caso sería divisor de \(1=j\,N-\gamma\).
Así que \(p\mid \rho=\gcd(n,r)\).
En particular, \(p\leq \rho \).
Sea \(j\) el mínimo número \(0\leq j<\rho\) tal que \(p\mid_{\mathbb Z} (1+jN)\).
Si \(j=0\), entonces estamos ante la unidad \(1\), que no es divisible por \(p\).
Por lo tanto, asumimos que \(j>0\).
Sea \(j = qp+c\), dado por el algoritmo de la división, al dividir por \(p\).
Si \(c>0\), se tiene \(c<j\), y además
\(p\mid_{\mathbb Z} (1+jN-pqN) = (1+cN)\), contradiciendo la minimalidad de \(j\)
si es que era \(q>0\), por lo tanto \(q = 0,j<p\).
LO QUE SIGUE YA NO SIRVE...
Luego \(c=0\).
Pero en ese caso \(j\) es múltiplo de \(p\), lo que implica que \(p\mid 1\), absurdo.
Por lo tanto, no hay divisores primos comunes de \(n\) y alguno de los números
\(\gamma=\gamma_j=(1+jN)\).
Es decir, son coprimos, y por lo tanto todos ellos son unidades.