title: Euler's four-square identity
text: # Euler's four-square identity
  
  In mathematics, Euler's four-square identity says that the product of two numbers, each of which is a sum of four squares, is itself a sum of four squares.
  
  ## Algebraic identity
  
  For any pair of quadruples from a commutative ring, the following expressions are equal:
  
  ${\begin{aligned}&\left(a_{1}^{2}+a_{2}^{2}+a_{3}^{2}+a_{4}^{2}\right)\left(b_{1}^{2}+b_{2}^{2}+b_{3}^{2}+b_{4}^{2}\right)\\[3mu]&\qquad =\left(a_{1}b_{1}-a_{2}b_{2}-a_{3}b_{3}-a_{4}b_{4}\right)^{2}+\left(a_{1}b_{2}+a_{2}b_{1}+a_{3}b_{4}-a_{4}b_{3}\right)^{2}\\[3mu]&\qquad \qquad +\left(a_{1}b_{3}-a_{2}b_{4}+a_{3}b_{1}+a_{4}b_{2}\right)^{2}+\left(a_{1}b_{4}+a_{2}b_{3}-a_{3}b_{2}+a_{4}b_{1}\right)^{2}.\end{aligned}}$
  
  Euler wrote about this identity in a letter dated May 4, 1748 to Goldbach (but he used a different sign convention from the above). It can be verified with elementary algebra.
  
  The identity was used by Lagrange to prove his four square theorem. More specifically, it implies that it is sufficient to prove the theorem for prime numbers, after which the more general theorem follows. The sign convention used above corresponds to the signs obtained by multiplying two quaternions. Other sign conventions can be obtained by changing any $a_{k}$ to $-a_{k}$, and/or any $b_{k}$ to $-b_{k}$.
  
  If the $a_{k}$ and $b_{k}$ are real numbers, the identity expresses the fact that the absolute value of the product of two quaternions is equal to the product of their absolute values, in the same way that the Brahmagupta–Fibonacci two-square identity does for complex numbers. This property is the definitive feature of composition algebras.
  
  Hurwitz's theorem states that an identity of form,
  
  $\left(a_{1}^{2}+a_{2}^{2}+a_{3}^{2}+\dots +a_{n}^{2}\right)\left(b_{1}^{2}+b_{2}^{2}+b_{3}^{2}+\dots +b_{n}^{2}\right)=c_{1}^{2}+c_{2}^{2}+c_{3}^{2}+\dots +c_{n}^{2}$
  
  where the $c_{i}$ are bilinear functions of the $a_{i}$ and $b_{i}$ is possible only for n = 1, 2, 4, or 8.
  
  ### Proof of the identity using quaternions
  
  Comment: The proof of Euler's four-square identity is by simple algebraic evaluation. Quaternions derive from the four-square identity, which can be written as the product of two inner products of 4-dimensional vectors, yielding again an inner product of 4-dimensional vectors: (a·a)(b·b) = (a×b)·(a×b). This defines the quaternion multiplication rule a×b, which simply reflects Euler's identity, and some mathematics of quaternions. Quaternions are, so to say, the "square root" of the four-square identity. But let the proof go on:
  
  Let $\alpha =a_{1}+a_{2}i+a_{3}j+a_{4}k$ and $\beta =b_{1}+b_{2}i+b_{3}j+b_{4}k$ be a pair of quaternions. Their quaternion conjugates are $\alpha ^{*}=a_{1}-a_{2}i-a_{3}j-a_{4}k$ and $\beta ^{*}=b_{1}-b_{2}i-b_{3}j-b_{4}k$. Then
  
  $A:=\alpha \alpha ^{*}=a_{1}^{2}+a_{2}^{2}+a_{3}^{2}+a_{4}^{2}$
  
  and
  
  $B:=\beta \beta ^{*}=b_{1}^{2}+b_{2}^{2}+b_{3}^{2}+b_{4}^{2}.$
  
  The product of these two is $AB=\alpha \alpha ^{*}\beta \beta ^{*}$, where $\beta \beta ^{*}$ is a real number, so it can commute with the quaternion $\alpha ^{*}$, yielding
  
  $AB=\alpha \beta \beta ^{*}\alpha ^{*}.$
  
  No parentheses are necessary above, because quaternions associate. The conjugate of a product is equal to the commuted product of the conjugates of the product's factors, so
  
  $AB=\alpha \beta (\alpha \beta )^{*}=\gamma \gamma ^{*}$
  
  where $\gamma$ is the Hamilton product of $\alpha$ and $\beta$:
  
  ${\begin{aligned}\gamma &=\left(a_{1}+\langle a_{2},a_{3},a_{4}\rangle \right)\left(b_{1}+\langle b_{2},b_{3},b_{4}\rangle \right)\\[3mu]&=a_{1}b_{1}+a_{1}\langle b_{2},\ b_{3},\ b_{4}\rangle +\langle a_{2},\ a_{3},\ a_{4}\rangle b_{1}+\langle a_{2},\ a_{3},\ a_{4}\rangle \langle b_{2},\ b_{3},\ b_{4}\rangle \\[3mu]&=a_{1}b_{1}+\langle a_{1}b_{2},\ a_{1}b_{3},\ a_{1}b_{4}\rangle +\langle a_{2}b_{1},\ a_{3}b_{1},\ a_{4}b_{1}\rangle \\&\qquad -\langle a_{2},\ a_{3},\ a_{4}\rangle \cdot \langle b_{2},\ b_{3},\ b_{4}\rangle +\langle a_{2},\ a_{3},\ a_{4}\rangle \times \langle b_{2},\ b_{3},\ b_{4}\rangle \\[3mu]&=a_{1}b_{1}+\langle a_{1}b_{2}+a_{2}b_{1},\ a_{1}b_{3}+a_{3}b_{1},\ a_{1}b_{4}+a_{4}b_{1}\rangle \\&\qquad -a_{2}b_{2}-a_{3}b_{3}-a_{4}b_{4}+\langle a_{3}b_{4}-a_{4}b_{3},\ a_{4}b_{2}-a_{2}b_{4},\ a_{2}b_{3}-a_{3}b_{2}\rangle \\[3mu]&=(a_{1}b_{1}-a_{2}b_{2}-a_{3}b_{3}-a_{4}b_{4})\\&\qquad +\langle a_{1}b_{2}+a_{2}b_{1}+a_{3}b_{4}-a_{4}b_{3},\ a_{1}b_{3}+a_{3}b_{1}+a_{4}b_{2}-a_{2}b_{4},\ a_{1}b_{4}+a_{4}b_{1}+a_{2}b_{3}-a_{3}b_{2}\rangle \\[3mu]\gamma &=(a_{1}b_{1}-a_{2}b_{2}-a_{3}b_{3}-a_{4}b_{4})+(a_{1}b_{2}+a_{2}b_{1}+a_{3}b_{4}-a_{4}b_{3})i\\&\qquad +(a_{1}b_{3}+a_{3}b_{1}+a_{4}b_{2}-a_{2}b_{4})j+(a_{1}b_{4}+a_{4}b_{1}+a_{2}b_{3}-a_{3}b_{2})k.\end{aligned}}$
  
  Then
  
  ${\begin{aligned}\gamma ^{*}&=(a_{1}b_{1}-a_{2}b_{2}-a_{3}b_{3}-a_{4}b_{4})-(a_{1}b_{2}+a_{2}b_{1}+a_{3}b_{4}-a_{4}b_{3})i\\&\qquad -(a_{1}b_{3}+a_{3}b_{1}+a_{4}b_{2}-a_{2}b_{4})j-(a_{1}b_{4}+a_{4}b_{1}+a_{2}b_{3}-a_{3}b_{2})k.\end{aligned}}$
  
  If $\gamma =r+{\vec {u}}$ where $r$ is the scalar part and ${\vec {u}}=\langle u_{1},u_{2},u_{3}\rangle$ is the vector part, then $\gamma ^{*}=r-{\vec {u}}$ so
  
  ${\begin{aligned}\gamma \gamma ^{*}&=(r+{\vec {u}})(r-{\vec {u}})=r^{2}-r{\vec {u}}+r{\vec {u}}-{\vec {u}}{\vec {u}}=r^{2}+{\vec {u}}\cdot {\vec {u}}-{\vec {u}}\times {\vec {u}}\\&=r^{2}+{\vec {u}}\cdot {\vec {u}}=r^{2}+u_{1}^{2}+u_{2}^{2}+u_{3}^{2}.\end{aligned}}$
  
  So,
  
  ${\begin{aligned}AB=\gamma \gamma ^{*}&=(a_{1}b_{1}-a_{2}b_{2}-a_{3}b_{3}-a_{4}b_{4})^{2}+(a_{1}b_{2}+a_{2}b_{1}+a_{3}b_{4}-a_{4}b_{3})^{2}\\&\qquad +(a_{1}b_{3}+a_{3}b_{1}+a_{4}b_{2}-a_{2}b_{4})^{2}+(a_{1}b_{4}+a_{4}b_{1}+a_{2}b_{3}-a_{3}b_{2})^{2}.\end{aligned}}$
  
  ## Pfister's identity
  
  Pfister found another square identity for any even power:
  
  If the $c_{i}$ are just rational functions of one set of variables, so that each $c_{i}$ has a denominator, then it is possible for all $n=2^{m}$.
  
  Thus, another four-square identity is as follows:
  
  ${\begin{aligned}&\left(a_{1}^{2}+a_{2}^{2}+a_{3}^{2}+a_{4}^{2}\right)\left(b_{1}^{2}+b_{2}^{2}+b_{3}^{2}+b_{4}^{2}\right)\\[5mu]&\quad =\left(a_{1}b_{4}+a_{2}b_{3}+a_{3}b_{2}+a_{4}b_{1}\right)^{2}+\left(a_{1}b_{3}-a_{2}b_{4}+a_{3}b_{1}-a_{4}b_{2}\right)^{2}\\&\quad \qquad +\left(a_{1}b_{2}+a_{2}b_{1}+{\frac {a_{3}u_{1}}{b_{1}^{2}+b_{2}^{2}}}-{\frac {a_{4}u_{2}}{b_{1}^{2}+b_{2}^{2}}}\right)^{2}+\left(a_{1}b_{1}-a_{2}b_{2}-{\frac {a_{4}u_{1}}{b_{1}^{2}+b_{2}^{2}}}-{\frac {a_{3}u_{2}}{b_{1}^{2}+b_{2}^{2}}}\right)^{2}\end{aligned}}$
  
  where $u_{1}$ and $u_{2}$ are given by
  
  ${\begin{aligned}u_{1}&=b_{1}^{2}b_{4}-2b_{1}b_{2}b_{3}-b_{2}^{2}b_{4}\\u_{2}&=b_{1}^{2}b_{3}+2b_{1}b_{2}b_{4}-b_{2}^{2}b_{3}\end{aligned}}$
  
  Incidentally, the following identity is also true:
  
  $u_{1}^{2}+u_{2}^{2}=\left(b_{1}^{2}+b_{2}^{2}\right)^{2}\left(b_{3}^{2}+b_{4}^{2}\right)$
  
  