In algebra, Pfister's sixteen-square identity is a non-bilinear identity of form
(
x
1
2
+
x
2
2
+
x
3
2
+
⋯
+
x
16
2
)
(
y
1
2
+
y
2
2
+
y
3
2
+
⋯
+
y
16
2
)
=
z
1
2
+
z
2
2
+
z
3
2
+
⋯
+
z
16
2
It was first proven to exist by H. Zassenhaus and W. Eichhorn in the 1960s, and independently by Pfister around the same time. There are several versions, a concise one of which is
z
1
=
x
1
y
1
−
x
2
y
2
−
x
3
y
3
−
x
4
y
4
−
x
5
y
5
−
x
6
y
6
−
x
7
y
7
−
x
8
y
8
+
u
1
y
9
−
u
2
y
10
−
u
3
y
11
−
u
4
y
12
−
u
5
y
13
−
u
6
y
14
−
u
7
y
15
−
u
8
y
16
z
2
=
x
2
y
1
+
x
1
y
2
+
x
4
y
3
−
x
3
y
4
+
x
6
y
5
−
x
5
y
6
−
x
8
y
7
+
x
7
y
8
+
u
2
y
9
+
u
1
y
10
+
u
4
y
11
−
u
3
y
12
+
u
6
y
13
−
u
5
y
14
−
u
8
y
15
+
u
7
y
16
z
3
=
x
3
y
1
−
x
4
y
2
+
x
1
y
3
+
x
2
y
4
+
x
7
y
5
+
x
8
y
6
−
x
5
y
7
−
x
6
y
8
+
u
3
y
9
−
u
4
y
10
+
u
1
y
11
+
u
2
y
12
+
u
7
y
13
+
u
8
y
14
−
u
5
y
15
−
u
6
y
16
z
4
=
x
4
y
1
+
x
3
y
2
−
x
2
y
3
+
x
1
y
4
+
x
8
y
5
−
x
7
y
6
+
x
6
y
7
−
x
5
y
8
+
u
4
y
9
+
u
3
y
10
−
u
2
y
11
+
u
1
y
12
+
u
8
y
13
−
u
7
y
14
+
u
6
y
15
−
u
5
y
16
z
5
=
x
5
y
1
−
x
6
y
2
−
x
7
y
3
−
x
8
y
4
+
x
1
y
5
+
x
2
y
6
+
x
3
y
7
+
x
4
y
8
+
u
5
y
9
−
u
6
y
10
−
u
7
y
11
−
u
8
y
12
+
u
1
y
13
+
u
2
y
14
+
u
3
y
15
+
u
4
y
16
z
6
=
x
6
y
1
+
x
5
y
2
−
x
8
y
3
+
x
7
y
4
−
x
2
y
5
+
x
1
y
6
−
x
4
y
7
+
x
3
y
8
+
u
6
y
9
+
u
5
y
10
−
u
8
y
11
+
u
7
y
12
−
u
2
y
13
+
u
1
y
14
−
u
4
y
15
+
u
3
y
16
z
7
=
x
7
y
1
+
x
8
y
2
+
x
5
y
3
−
x
6
y
4
−
x
3
y
5
+
x
4
y
6
+
x
1
y
7
−
x
2
y
8
+
u
7
y
9
+
u
8
y
10
+
u
5
y
11
−
u
6
y
12
−
u
3
y
13
+
u
4
y
14
+
u
1
y
15
−
u
2
y
16
z
8
=
x
8
y
1
−
x
7
y
2
+
x
6
y
3
+
x
5
y
4
−
x
4
y
5
−
x
3
y
6
+
x
2
y
7
+
x
1
y
8
+
u
8
y
9
−
u
7
y
10
+
u
6
y
11
+
u
5
y
12
−
u
4
y
13
−
u
3
y
14
+
u
2
y
15
+
u
1
y
16
z
9
=
x
9
y
1
−
x
10
y
2
−
x
11
y
3
−
x
12
y
4
−
x
13
y
5
−
x
14
y
6
−
x
15
y
7
−
x
16
y
8
+
x
1
y
9
−
x
2
y
10
−
x
3
y
11
−
x
4
y
12
−
x
5
y
13
−
x
6
y
14
−
x
7
y
15
−
x
8
y
16
z
10
=
x
10
y
1
+
x
9
y
2
+
x
12
y
3
−
x
11
y
4
+
x
14
y
5
−
x
13
y
6
−
x
16
y
7
+
x
15
y
8
+
x
2
y
9
+
x
1
y
10
+
x
4
y
11
−
x
3
y
12
+
x
6
y
13
−
x
5
y
14
−
x
8
y
15
+
x
7
y
16
z
11
=
x
11
y
1
−
x
12
y
2
+
x
9
y
3
+
x
10
y
4
+
x
15
y
5
+
x
16
y
6
−
x
13
y
7
−
x
14
y
8
+
x
3
y
9
−
x
4
y
10
+
x
1
y
11
+
x
2
y
12
+
x
7
y
13
+
x
8
y
14
−
x
5
y
15
−
x
6
y
16
z
12
=
x
12
y
1
+
x
11
y
2
−
x
10
y
3
+
x
9
y
4
+
x
16
y
5
−
x
15
y
6
+
x
14
y
7
−
x
13
y
8
+
x
4
y
9
+
x
3
y
10
−
x
2
y
11
+
x
1
y
12
+
x
8
y
13
−
x
7
y
14
+
x
6
y
15
−
x
5
y
16
z
13
=
x
13
y
1
−
x
14
y
2
−
x
15
y
3
−
x
16
y
4
+
x
9
y
5
+
x
10
y
6
+
x
11
y
7
+
x
12
y
8
+
x
5
y
9
−
x
6
y
10
−
x
7
y
11
−
x
8
y
12
+
x
1
y
13
+
x
2
y
14
+
x
3
y
15
+
x
4
y
16
z
14
=
x
14
y
1
+
x
13
y
2
−
x
16
y
3
+
x
15
y
4
−
x
10
y
5
+
x
9
y
6
−
x
12
y
7
+
x
11
y
8
+
x
6
y
9
+
x
5
y
10
−
x
8
y
11
+
x
7
y
12
−
x
2
y
13
+
x
1
y
14
−
x
4
y
15
+
x
3
y
16
z
15
=
x
15
y
1
+
x
16
y
2
+
x
13
y
3
−
x
14
y
4
−
x
11
y
5
+
x
12
y
6
+
x
9
y
7
−
x
10
y
8
+
x
7
y
9
+
x
8
y
10
+
x
5
y
11
−
x
6
y
12
−
x
3
y
13
+
x
4
y
14
+
x
1
y
15
−
x
2
y
16
z
16
=
x
16
y
1
−
x
15
y
2
+
x
14
y
3
+
x
13
y
4
−
x
12
y
5
−
x
11
y
6
+
x
10
y
7
+
x
9
y
8
+
x
8
y
9
−
x
7
y
10
+
x
6
y
11
+
x
5
y
12
−
x
4
y
13
−
x
3
y
14
+
x
2
y
15
+
x
1
y
16
If all
x
i
and
y
i
with
i
>
8
are set equal to zero, then it reduces to Degen's eight-square identity (in blue). The
u
i
are
u
1
=
(
a
x
1
2
+
x
2
2
+
x
3
2
+
x
4
2
+
x
5
2
+
x
6
2
+
x
7
2
+
x
8
2
)
x
9
−
2
x
1
(
b
x
1
x
9
+
x
2
x
10
+
x
3
x
11
+
x
4
x
12
+
x
5
x
13
+
x
6
x
14
+
x
7
x
15
+
x
8
x
16
)
c
u
2
=
(
x
1
2
+
a
x
2
2
+
x
3
2
+
x
4
2
+
x
5
2
+
x
6
2
+
x
7
2
+
x
8
2
)
x
10
−
2
x
2
(
x
1
x
9
+
b
x
2
x
10
+
x
3
x
11
+
x
4
x
12
+
x
5
x
13
+
x
6
x
14
+
x
7
x
15
+
x
8
x
16
)
c
u
3
=
(
x
1
2
+
x
2
2
+
a
x
3
2
+
x
4
2
+
x
5
2
+
x
6
2
+
x
7
2
+
x
8
2
)
x
11
−
2
x
3
(
x
1
x
9
+
x
2
x
10
+
b
x
3
x
11
+
x
4
x
12
+
x
5
x
13
+
x
6
x
14
+
x
7
x
15
+
x
8
x
16
)
c
u
4
=
(
x
1
2
+
x
2
2
+
x
3
2
+
a
x
4
2
+
x
5
2
+
x
6
2
+
x
7
2
+
x
8
2
)
x
12
−
2
x
4
(
x
1
x
9
+
x
2
x
10
+
x
3
x
11
+
b
x
4
x
12
+
x
5
x
13
+
x
6
x
14
+
x
7
x
15
+
x
8
x
16
)
c
u
5
=
(
x
1
2
+
x
2
2
+
x
3
2
+
x
4
2
+
a
x
5
2
+
x
6
2
+
x
7
2
+
x
8
2
)
x
13
−
2
x
5
(
x
1
x
9
+
x
2
x
10
+
x
3
x
11
+
x
4
x
12
+
b
x
5
x
13
+
x
6
x
14
+
x
7
x
15
+
x
8
x
16
)
c
u
6
=
(
x
1
2
+
x
2
2
+
x
3
2
+
x
4
2
+
x
5
2
+
a
x
6
2
+
x
7
2
+
x
8
2
)
x
14
−
2
x
6
(
x
1
x
9
+
x
2
x
10
+
x
3
x
11
+
x
4
x
12
+
x
5
x
13
+
b
x
6
x
14
+
x
7
x
15
+
x
8
x
16
)
c
u
7
=
(
x
1
2
+
x
2
2
+
x
3
2
+
x
4
2
+
x
5
2
+
x
6
2
+
a
x
7
2
+
x
8
2
)
x
15
−
2
x
7
(
x
1
x
9
+
x
2
x
10
+
x
3
x
11
+
x
4
x
12
+
x
5
x
13
+
x
6
x
14
+
b
x
7
x
15
+
x
8
x
16
)
c
u
8
=
(
x
1
2
+
x
2
2
+
x
3
2
+
x
4
2
+
x
5
2
+
x
6
2
+
x
7
2
+
a
x
8
2
)
x
16
−
2
x
8
(
x
1
x
9
+
x
2
x
10
+
x
3
x
11
+
x
4
x
12
+
x
5
x
13
+
x
6
x
14
+
x
7
x
15
+
b
x
8
x
16
)
c
and,
a
=
−
1
,
b
=
0
,
c
=
x
1
2
+
x
2
2
+
x
3
2
+
x
4
2
+
x
5
2
+
x
6
2
+
x
7
2
+
x
8
2
.
The identity shows that, in general, the product of two sums of sixteen squares is the sum of sixteen rational squares. Incidentally, the
u
i
also obey,
u
1
2
+
u
2
2
+
u
3
2
+
u
4
2
+
u
5
2
+
u
6
2
+
u
7
2
+
u
8
2
=
x
9
2
+
x
10
2
+
x
11
2
+
x
12
2
+
x
13
2
+
x
14
2
+
x
15
2
+
x
16
2
No sixteen-square identity exists involving only bilinear functions since Hurwitz's theorem states an identity of the form
(
x
1
2
+
x
2
2
+
x
3
2
+
⋯
+
x
n
2
)
(
y
1
2
+
y
2
2
+
y
3
2
+
⋯
+
y
n
2
)
=
z
1
2
+
z
2
2
+
z
3
2
+
⋯
+
z
n
2
with the
z
i
bilinear functions of the
x
i
and
y
i
is possible only for n ∈ {1, 2, 4, 8} . However, the more general Pfister's theorem (1965) shows that if the
z
i
are rational functions of one set of variables, hence has a denominator, then it is possible for all
n
=
2
m
. There are also non-bilinear versions of Euler's four-square and Degen's eight-square identities.