For odd $n$ the Fermat equation $x^n + y^n = z^n$ factorises as $$(x + y)(x + \zeta y) \cdots (x + \zeta^{n-1}y) = z^n,$$
where $\zeta = e^{2 \pi i/n}$. I tried seeing this was true by multiplying the factors out by hand but got into a mess. Is there a neat way of showing this (e.g. is there a useful identity that can be used involving products of the $\zeta^i$ etc.?
We will make use of the following theorem
Theorem Let $P(x)$ be a polynomial and $P(z)=0$, then $x-z$ divides $P(x)$.
Let $n$ be odd and write $$P(x) = x^n + y^n,$$ then $P(-\zeta^r y)=0$ for all $n$ powers of $r$.
Therefore $(x+y)(x+\zeta y)(x+\zeta^2 y)\cdots$ divides $P(x)$ but it has the same degree so its equal (up to a constant factor which is easily seen to be 1).
Let $n$ be odd. The $n$ distinct roots of the polynomial $X^n-1$ are $1,\zeta,\dots,\zeta^{n-1}$ $$ X^n-1=(X-1)(X-\zeta)\cdots(X-\zeta^{n-1}) $$ Apply the transformation $X\mapsto -\frac{x}{y}$ $$ \left(-\frac{x}{y}\right)^n-1=\left(-\frac{x}{y}-1\right)\left(-\frac{x}{y}-\zeta\right)\cdots\left(-\frac{x}{y}-\zeta^{n-1}\right) $$ Finally, multiply through by $-y^n$ $$ x^n+y^n=(x+y)\cdots(x+\zeta^{n-1}y) $$
Fix $x \neq 0$. Let
$$P(y)=(x + y)(x + \zeta y) \cdots (x + \zeta^{n-1}y) -x^n-y^n$$
Then $P(-x)=P(-\zeta x)=P(-\zeta^2 x)=..=P(-\zeta^{n-1} x)=0$ and $P(0)=x^n-x^n=0$.
Thus $P$ is a polynomial of degree at most $n$ which has at least $n+1$ roots. This proves that
$$(x + y)(x + \zeta y) \cdots (x + \zeta^{n-1}y) =x^n+y^n$$
for all $y$ and all $x \neq 0$.
The case $x=0$ is trivial.
