Avevo deciso di postare la mia soluzione dopo il 15 marzo ma ormai...
Indico con $ \displaystyle \epsilon_k $ le radici n-esime di 1,con k variabile da 0 ad n-1( mi trovo meglio

).
Sl ha:
$ \displaystyle x^n-1=\prod_{k=0}^{n-1}(x-\epsilon_k) $
Cambiando x in -x :
$ \displaystyle (-x)^n-1=\prod_{k=0}^{n-1}(-x-\epsilon_k)=(-1)^n \prod_{k=0}^{n-1}(x+\epsilon_k ) $
Moltiplicando risulta:
(*) $ \displaystyle (x^n-1)((-x)^n-1)=(-1)^n\prod_{k=0}^{n-1}(x^2-\epsilon_k ^2)} $
Dintinguiamo ora i casi n pari (n=2p) ed n dispari ( n=2p+1)
(A) n=2p
La
(*) diventa ora:
$ {\displaystyle (x^{2p}-1)^2=\prod_{k=0}^{2p-1}(x^2-\epsilon_k ^2)} $
Poniamo in questa $ \displaystyle x=i\sqrt{\frac{a}{b}} $ dove
i è l'unità immaginaria:
$ \displaystyle [(-1)^p \frac{a^p}{b^p}-1]^2=\prod_{k=0}^{2p-1}(-\frac{a}{b}-\epsilon_k^2)=\frac{(-1)^{2p}}{b^{2p}}\prod_{k=0}^{2p-1}(a+b\epsilon_k^2) $
Ovvero:
$ \displaystyle [(-1)^p a^p-b^p]^2=\prod_{k=0}^{2p-1}(a+b\epsilon_k^2) $
E reintroducendo n si ha la formula per n pari.
FORMULA PER n PARI :
$ \fbox{\begin{equation} \displaystyle \prod_{k=0}^{n-1}(a+b\epsilon_k^2)= ((-1)^{\frac{n}{2}} a^{\frac{n}{2}}-b^{\frac{n}{2}})^2 \end{equation }} $
(B) n=2p+1
In questo caso la
(*) diventa ora :
$ \displaystyle (x^{2p+1)}-1)(-x^{2p+1}-1)=-\prod_{k=0}^{2p}(x^2-\epsilon_k ^2)} $
Ovvero:
$ \displaystyle x^{4p+2}-1=\prod_{k=0}^{2p}(x^2-\epsilon_k ^2)} $
Come prima ,ponendo $ \displaystyle x=i\sqrt{\frac{a}{b}} $,abbiamo:
$ \displaystyle (-1)^{2p+1}(\frac{a}{b})^{2p+1}-1=\prod_{k=0}^{2p}(-\frac{a}{b}-\epsilon_k^2) $
Oppure:
$ \displaystyle -\frac{a^{2p+1}+b^{2p+1}}{b^{2p+1}}=\frac{(-1)^{2p+1}}{b^{2p+1}}\prod_{k=0}^{2p}(a+b\epsilon_k^2) $
Semplificando e ritornando ad n si ha la formula anche per n dispari.
FORMULA PER n DISPARI:
$ \fbox{\begin{equation} \displaystyle \prod_{k=0}^{n-1}(a+b\epsilon_k^2)=a^n+b^n \end{equation } } $