A Pythagorean quadruple is a tuple of integers a, b, c, and d, such that a2 + b2 + c2 = d2. They are solutions of a Diophantine equation and often only positive integer values are considered. However, to provide a more complete geometric interpretation, the integer values can be allowed to be negative and zero (thus allowing Pythagorean triples to be included) with the only condition being that d > 0. In this setting, a Pythagorean quadruple (a, b, c, d) defines a cuboid with integer side lengths |a|, |b|, and |c|, whose space diagonal has integer length d; with this interpretation, Pythagorean quadruples are thus also called Pythagorean boxes. In this article we will assume, unless otherwise stated, that the values of a Pythagorean quadruple are all positive integers.
Parametrization of primitive quadruples A Pythagorean quadruple is called primitive if the greatest common divisor of its entries is 1. Every Pythagorean quadruple is an integer multiple of a primitive quadruple. The set of primitive Pythagorean quadruples for which a is odd can be generated by the formulas
a = m 2 + n 2 − p 2 − q 2 , b = 2 ( m q + n p ) , c = 2 ( n q − m p ) , d = m 2 + n 2 + p 2 + q 2 , {\displaystyle {\begin{aligned}a&=m^{2}+n^{2}-p^{2}-q^{2},\\b&=2(mq+np),\\c&=2(nq-mp),\\d&=m^{2}+n^{2}+p^{2}+q^{2},\end{aligned}}}
where m, n, p, q are non-negative integers with greatest common divisor 1 such that m + n + p + q is odd. Thus, all primitive Pythagorean quadruples are characterized by the identity
( m 2 + n 2 + p 2 + q 2 ) 2 = ( 2 m q + 2 n p ) 2 + ( 2 n q − 2 m p ) 2 + ( m 2 + n 2 − p 2 − q 2 ) 2 . {\displaystyle (m^{2}+n^{2}+p^{2}+q^{2})^{2}=(2mq+2np)^{2}+(2nq-2mp)^{2}+(m^{2}+n^{2}-p^{2}-q^{2})^{2}.}
Alternate parametrization All Pythagorean quadruples (including non-primitives, and with repetition, though a, b, and c do not appear in all possible orders) can be generated from two positive integers a and b as follows: If a and b have different parity, let p be any factor of a2 + b2 such that p2 < a2 + b2. Then c = a2 + b2 − p2/2p and d = a2 + b2 + p2/2p. Note that p = d − c. A similar method exists for generating all Pythagorean quadruples for which a and b are both even. Let l = a/2 and m = b/2 and let n be a factor of l2 + m2 such that n2 < l2 + m2. Then c = l2 + m2 − n2/n and d = l2 + m2 + n2/n. This method generates all Pythagorean quadruples exactly once each when l and m run through all pairs of natural numbers and n runs through all permissible values for each pair. No such method exists if both a and b are odd, in which case no solutions exist as can be seen by the parametrization in the previous section.
… excerpt ends here. Continue reading the full article.


