Eine etwas andere Idee
Mache mal mit n+1 eine Polynomdivision
(x^{n+1} - y^{n+1}) : (x - y) = x^n + y * (x^n - y^n) / (x - y)
Jetzt sieht man auf der Rechten seite den Quotienten (x^n - y^n) / (x - y). Der ist ganzzahlig wenn x^n - y^n durch (x - y) zu teilen ist. Das soll aber nach der Induktionsannahme der Fall sein. Daher ist dann auch der Term auf der Linken Seite ganzzahlig und somit ist auch x^{n+1} - y^{n+1} durch (x - y) teilbar.