![Smile :)](./images/smilies/icon_smile.gif)
Entrambe le dimostrazioni che conosco (che in effetti sono filosoficamente la stessa) arrivano a quella che hai scritto anche tu
$ \omega^n=a\pm \omega$
con $a$ intero.
Questa, vogliamo abbatterla a colpi di congruenze, che non promette bene perché di soluzioni ce ne sono, e anche non tanto banali (n=13). Quando ci sono di mezzo equazioni esponenziali, la speranza è di prendere congruenze modulo cose piuttosto grandi, fare in modo che le quantità in gioco diventino periodiche ad un certo punto, ma le soluzioni si trovino tutte nell'antiperiodo. (Spero di essermi spiegato. Esempio idiota: $2^n=8x+4$, modulo potenze di 2 si conclude).
Uno può lanciarsi in ardite congruenze già qui, con la complicazione di trovarci in $\mathbb{Z}[\omega]$, e finire per fare più o meno quello che ti era venuto in mente a scuola, come fa Nagell; oppure si può andare "al piano di sotto" con la successione $b_n$ che hai definito.
Per vedere che $b_n$ vale $\pm 1$ solo in quei casi lì, proviamo a costringere n in varie sottosuccessioni, fino a raggiungere un assurdo; a quel punto gli unici casi possibili saranno solo quei pochi casi "bassi" che ci eravamo lasciati da parte da fare a mano nel procedere della dimostrazione.
Primo passo: come si comporta $b_n$ modulo 16 ?
Visto che puntiamo a studiare $b_n$ quando n è in particolari sottosuccessioni, è utile ricavarci qualche formula per cose come $b_{nk}$, $b_{nk+1}$.
Si possono ricavare, per esempio, da $\omega^n=b_{n+1}-b_n \bar \omega$ elevando alla $k$ e sviluppando il binomio sulla destra, e poi liberandoci delle $\omega$ per rimanere con una complicata espressione per ricorrenza sugli interi che esprima $b_{nk+1}$ in termini di $b_n$, $b_{n+1}$ e i $b_j$ fino a $k$.
Buon lavoro.