Seeing the answer by Andre Nicolas to this question, it got me thinking: is there a way one can derive the formula for $a^n-b^n$ in a more general setting than an euclidian domain, or at least the fact that $a^n-b^n$ is divisible by $a-b$? It seems like all there's needed for it to work is that the ring be commutative, because then: $$ (a-b)(\sum\limits_{k=1}^n a^{n-k}b^{k-1}) = \sum\limits_{k=1}^n a^{n-k+1}b^{k-1} - \sum\limits_{k=1}^n a^{n-k} b^k = a^n - b^n, $$ commutativity being used for the second equality.
Of course, one can do the above check, but it's not nearly as elegant as using the existence of a division algorithm for example, as in the answer I mentioned. So I'm asking if there's a way to at least deduce the divisibility by $a-b$ while also refusing to guess what the formula is.