Pythagorean triple

Generating a triple

Euclid's formula is a fundamental formula for generating Pythagorean triples given an arbitrary pair of integers m and n with m > n > 0. The formula states that the integers

$${\displaystyle a=m^2-n^2,b=2mn,c=m^2+n^2}$$

form a Pythagorean triple. For example, the integers

$${\displaystyle m=2,n=1}$$

generate the primitive triple (3, 4, 5):

$${\displaystyle a=2^2-1^2=3,b=2\times 2\times 1=4,c=2^2+1^2=5.}$$

The triple generated by Euclid's formula is primitive if and only if m and n are coprime and exactly one of them is even. When both m and n are odd, then a, b, and c will be even, and the triple will not be primitive; however, dividing a, b, and c by 2 will yield a primitive triple when m and n are coprime.

Every primitive triple arises (after the exchange of a and b, if a is even) from a unique pair of coprime numbers m, n, one of which is even. It follows that there are infinitely many primitive Pythagorean triples. This relationship of a, b and c to m and n from Euclid's formula is referenced throughout the rest of this article.

Despite generating all primitive triples, Euclid's formula does not produce all triples—for example, (9, 12, 15) cannot be generated using integer m and n. This can be remedied by inserting an additional parameter k to the formula: every Pythagorean triple is generated uniquely by

$${\displaystyle a=k\cdot (m^2-n^2),b=k\cdot (2mn),c=k\cdot (m^2+n^2)}$$

where m, n, and k are positive integers with m > n, and with m and n coprime and not both odd.

That these formulas generate Pythagorean triples can be verified by expanding a² + b² using elementary algebra and verifying that the result equals c². Since every Pythagorean triple can be divided through by some integer k to obtain a primitive triple, every triple can be generated uniquely by using the formula with m and n to generate its primitive counterpart and then multiplying through by k as in the last equation.

Choosing m and n from certain integer sequences gives interesting results. For example, if m and n are consecutive Pell numbers, a and b will differ by 1.

Many formulas for generating triples with particular properties have been developed since the time of Euclid.

Proof of Euclid's formula

That satisfaction of Euclid's formula by a, b, c is sufficient for the triangle to be Pythagorean is apparent from the fact that for positive integers m and n, m > n, the a, b, and c given by the formula are all positive integers, and from the fact that

$${\displaystyle a^2+b^2=(m^2-n^2{)}^2+(2mn{)}^2=(m^2+n^2{)}^2=c^2.}$$

A proof of the necessity that a, b, c be expressed by Euclid's formula for any primitive Pythagorean triple is as follows. All such primitive triples can be written as (a, b, c) where a² + b² = c² and a, b, c are coprime. Thus a, b, c are pairwise coprime (if a prime number divided two of them, it would be forced also to divide the third one). As a and b are coprime, at least one of them is odd. If we suppose that a is odd, then b is even and c is odd (if both a and b were odd, c would be even, and c² would be a multiple of 4, while a² + b² would be congruent to 2 modulo 4, as an odd square is congruent to 1 modulo 4).

From ${\displaystyle a^2+b^2=c^2,}$ assume a is odd. We obtain ${\displaystyle c^2-a^2=b^2}$ and hence ${\displaystyle (c-a)(c+a)=b^2.}$ Then ${\displaystyle \frac{(c+a)}{b}=\frac{b}{(c-a)}.}$ Since ${\displaystyle \frac{(c+a)}{b}}$ is rational, we set it equal to ${\displaystyle \frac{m}{n}}$ in lowest terms. Thus ${\displaystyle \frac{(c-a)}{b}=\frac{n}{m},}$ being the reciprocal of ${\displaystyle \frac{(c+a)}{b}.}$ Then solving

$${\displaystyle \frac{c}{b}+\frac{a}{b}=\frac{m}{n},\frac{c}{b}-\frac{a}{b}=\frac{n}{m}}$$

for ${\displaystyle \frac{c}{b}}$ and ${\displaystyle \frac{a}{b}}$ gives

$${\displaystyle \frac{c}{b}=\frac{1}{2}\left(\frac{m}{n}+\frac{n}{m}\right)=\frac{m^2+n^2}{2mn},\frac{a}{b}=\frac{1}{2}\left(\frac{m}{n}-\frac{n}{m}\right)=\frac{m^2-n^2}{2mn}.}$$

As ${\displaystyle \frac{m}{n}}$ is fully reduced, m and n are coprime, and they cannot both be even. If they were both odd, the numerator of ${\displaystyle \frac{m^2-n^2}{2mn}}$ would be a multiple of 4 (because an odd square is congruent to 1 modulo 4), and the denominator 2mn would not be a multiple of 4. Since 4 would be the minimum possible even factor in the numerator and 2 would be the maximum possible even factor in the denominator, this would imply a to be even despite defining it as odd. Thus one of m and n is odd and the other is even, and the numerators of the two fractions with denominator 2mn are odd. Thus these fractions are fully reduced (an odd prime dividing this denominator divides one of m and n but not the other; thus it does not divide m² ± n²). One may thus equate numerators with numerators and denominators with denominators, giving Euclid's formula $${\displaystyle a=m^2-n^2,b=2mn,c=m^2+n^2}$$ with m and n coprime and of opposite parities.

A longer but more commonplace proof is given in Maor (2007) and Sierpiński (2003). Another proof is given in Diophantine equation § Example of Pythagorean triples, as an instance of a general method that applies to every homogeneous Diophantine equation of degree two.

Interpretation of parameters in Euclid's formula

Suppose the sides of a Pythagorean triangle have lengths m² − n², 2mn, and m² + n², and suppose the angle between the leg of length m² − n² and the hypotenuse of length m² + n² is denoted as β. Then ${\displaystyle \tan \frac{\beta }{2}=\frac{n}{m}}$ and the full-angle trigonometric values are ${\displaystyle \sin \beta =\frac{2mn}{m^2+n^2}}$, ${\displaystyle \cos \beta =\frac{m^2-n^2}{m^2+n^2}}$, and ⁠${\displaystyle \tan \beta =\frac{2mn}{m^2-n^2}}$⁠.

A variant

The following variant of Euclid's formula is sometimes more convenient, as being more symmetric in m and n (same parity condition on m and n).

If m and n are two odd integers such that m > n, then $${\displaystyle a=mn,b=\frac{m^2-n^2}{2},c=\frac{m^2+n^2}{2}}$$ are three integers that form a Pythagorean triple, which is primitive if and only if m and n are coprime. Conversely, every primitive Pythagorean triple arises (after the exchange of a and b, if a is even) from a unique pair m > n > 0 of coprime odd integers.

Not exchanging a and b

In the presentation above, it is said that all Pythagorean triples are uniquely obtained from Euclid's formula "after the exchange of a and b, if a is even". Euclid's formula and the variant above can be merged as follows to avoid this exchange, leading to the following result.

Every primitive Pythagorean triple can be uniquely written $${\displaystyle a=2\epsilon mn,b=\epsilon (m^2-n^2),c=\epsilon (m^2+n^2),}$$ where m and n are positive coprime integers, and ${\displaystyle \epsilon =\frac{1}{2}}$ if m and n are both odd, and ${\displaystyle \epsilon =1}$ otherwise. Equivalently, ${\displaystyle \epsilon =\frac{1}{2}}$ if a is odd, and ${\displaystyle \epsilon =1}$ if a is even.

Elementary properties of primitive Pythagorean triples

General properties

The properties of a primitive Pythagorean triple (a, b, c) with a < b < c (without specifying which of a or b is even and which is odd) include:

• At most one of a, b, c is a square.
• The area of a Pythagorean triangle cannot be the square or twice the square of an integer.
• Exactly one of a, b is divisible by 2 (is even), and the hypotenuse c is always odd.
• Exactly one of a, b is divisible by 3, but never c.
• Exactly one of a, b is divisible by 4, but never c (because c is never even).
• Exactly one of a, b, c is divisible by 5.
• The largest number that always divides abc is 60.
• Any odd number of the form 2m+1, where m is an integer and m>1, can be the odd leg of a primitive Pythagorean triple. See almost-isosceles primitive Pythagorean triples section below. However, only even numbers divisible by 4 can be the even leg of a primitive Pythagorean triple. This is because Euclid's formula for the even leg given above is 2mn and one of m or n must be even.
• The hypotenuse c (which is always odd) is the sum of two squares. This requires all of its prime factors to be primes of the form 4n + 1. Therefore, c is of the form 4n + 1. A sequence of possible hypotenuse numbers for a primitive Pythagorean triple can be found at (sequence A008846 in the OEIS).
• The area (K = ab/2) is a congruent number divisible by 6.
• Only two sides of a primitive Pythagorean triple can be simultaneously prime because by Euclid's formula for generating a primitive Pythagorean triple, one of the legs must be composite and even. However, only one side can be an integer of perfect power ${\displaystyle p\ge 2}$ because if two sides were integers of perfect powers with equal exponent ${\displaystyle p}$ it would contradict the fact that there are no integer solutions to the Diophantine equation ${\displaystyle x^{2p}\pm y^{2p}=z^2}$, with ${\displaystyle x}$, ${\displaystyle y}$ and ${\displaystyle z}$ being pairwise coprime.
• There are no Pythagorean triangles in which the hypotenuse and one leg are the legs of another Pythagorean triangle; this is one of the equivalent forms of Fermat's right triangle theorem.
• No primitive Pythagorean triangle has an integer altitude from the hypotenuse; that is, every primitive Pythagorean triangle is indecomposable.
• The set of all primitive Pythagorean triples forms a rooted ternary tree in a natural way; see Tree of primitive Pythagorean triples.
• Neither of the acute angles of a Pythagorean triangle can be a rational number of degrees. (This follows from Niven's theorem.)

Special cases

In addition, special Pythagorean triples with certain additional properties can be guaranteed to exist:

• Every integer greater than 2 that is not congruent to 2 mod 4 (in other words, every integer greater than 2 which is not of the form 4k + 2) is part of a primitive Pythagorean triple. (If the integer has the form 4k, one may take n = 1 and m = 2k in Euclid's formula; if the integer is 2k + 1, one may take n = k and m = k + 1.)
• Every integer greater than 2 is part of a primitive or non-primitive Pythagorean triple. For example, the integers 6, 10, 14, and 18 are not part of primitive triples, but are part of the non-primitive triples (6, 8, 10), (14, 48, 50) and (18, 80, 82).
• There exist infinitely many Pythagorean triples in which the hypotenuse and the longest leg differ by exactly one. Such triples are necessarily primitive and have the form (2n + 1, 2n² + 2n, 2n² + 2n +1). This results from Euclid's formula by remarking that the condition implies that the triple is primitive and must verify (m² + n²) - 2mn = 1. This implies (m – n)² = 1, and thus m = n + 1. The above form of the triples results thus of substituting m for n + 1 in Euclid's formula.
• There exist infinitely many primitive Pythagorean triples in which the hypotenuse and the longest leg differ by exactly two. They are all primitive, and are obtained by putting n = 1 in Euclid's formula. More generally, for every integer k > 0, there exist infinitely many primitive Pythagorean triples in which the hypotenuse and the odd leg differ by 2k². They are obtained by putting n = k in Euclid's formula.
• There exist infinitely many Pythagorean triples in which the two legs differ by exactly one. For example, 20² + 21² = 29²; these are generated by Euclid's formula when ${\displaystyle \frac{m-n}{n}}$ is a convergent to ${\displaystyle \sqrt{2}.}$
• For each positive integer k, there exist k Pythagorean triples with different hypotenuses and the same area.
• For each positive integer k, there exist at least k different primitive Pythagorean triples with the same leg a, where a is some positive integer (the length of the even leg is 2mn, and it suffices to choose a with many factorizations, for example a = 4b, where b is a product of k different odd primes; this produces at least 2^k different primitive triples).
• For each positive integer k, there exist at least k different Pythagorean triples with the same hypotenuse.
• If c = p^e is a prime power, there exists a primitive Pythagorean triple a² + b² = c² if and only if the prime p has the form 4n + 1; this triple is unique up to the exchange of a and b.
• More generally, a positive integer c is the hypotenuse of a primitive Pythagorean triple if and only if each prime factor of c is congruent to 1 modulo 4; that is, each prime factor has the form 4n + 1. In this case, the number of primitive Pythagorean triples (a, b, c) with a < b is 2^k−1, where k is the number of distinct prime factors of c.
• There exist infinitely many Pythagorean triples with square numbers for both the hypotenuse c and the sum of the legs a + b. According to Fermat, the smallest such triple has sides a = 4,565,486,027,761; b = 1,061,652,293,520; and c = 4,687,298,610,289. Here a + b = 2,372,159² and c = 2,165,017². This is generated by Euclid's formula with parameter values m = 2,150,905 and n = 246,792.
• There exist non-primitive Pythagorean triangles with integer altitude from the hypotenuse. Such Pythagorean triangles are known as decomposable since they can be split along this altitude into two separate and smaller Pythagorean triangles.

Geometry of Euclid's formula

Rational points on a unit circle

Euclid's formula for a Pythagorean triple

$${\displaystyle a=m^2-n^2,b=2mn,c=m^2+n^2}$$

can be understood in terms of the geometry of rational points on the unit circle (Trautman 1998).

In fact, a point in the Cartesian plane with coordinates (x, y) belongs to the unit circle if x² + y² = 1. The point is rational if x and y are rational numbers, that is, if there are coprime integers a, b, c such that $${\displaystyle (\frac{a}{c}{)}^2+(\frac{b}{c}{)}^2=1.}$$

By multiplying both members by c², one can see that the rational points on the circle are in one-to-one correspondence with the primitive Pythagorean triples.

The unit circle may also be defined by a parametric equation $${\displaystyle x=\frac{1-t^2}{1+t^2},y=\frac{2t}{1+t^2}.}$$ Euclid's formula for Pythagorean triples and the inverse relationship t = y / (1 + x) mean that, except for (−1, 0), a point (x, y) on the circle is rational if and only if the corresponding value of t is a rational number. Note that t = y / (1 + x) = b / (c + a) = n / m is also the tangent of half of the angle that is opposite the triangle side of length b, and rational values of this half-angle tangent can be used to generate all primitive Pythagorean triples; see Generating all primitive Pythagorean triples using half-angle tangents.

Stereographic approach

There is a correspondence between points on the unit circle with rational coordinates and primitive Pythagorean triples. At this point, Euclid's formulae can be derived either by methods of trigonometry or equivalently by using the stereographic projection.

For the stereographic approach, suppose that P′ is a point on the x-axis with rational coordinates

$${\displaystyle P^{\prime }=\left(\frac{m}{n},0\right).}$$

Then, it can be shown by basic algebra that the point P has coordinates

$${\displaystyle P=\left(\frac{2\left(\frac{m}{n}\right)}{{\left(\frac{m}{n}\right)}^2+1},\frac{{\left(\frac{m}{n}\right)}^2-1}{{\left(\frac{m}{n}\right)}^2+1}\right)=\left(\frac{2mn}{m^2+n^2},\frac{m^2-n^2}{m^2+n^2}\right).}$$

This establishes that each rational point of the x-axis goes over to a rational point of the unit circle. The converse, that every rational point of the unit circle comes from such a point of the x-axis, follows by applying the inverse stereographic projection. Suppose that P(x, y) is a point of the unit circle with x and y rational numbers. Then the point P′ obtained by stereographic projection onto the x-axis has coordinates

$${\displaystyle \left(\frac{x}{1-y},0\right)}$$

which is rational.

In terms of algebraic geometry, the algebraic variety of rational points on the unit circle is birational to the affine line over the rational numbers. The unit circle is thus called a rational curve, and it is this fact which enables an explicit parameterization of the (rational number) points on it by means of rational functions.

Related triangle geometry

When the lengths of the sides of a right triangle are positive integers, a number of other lengths, including the diameter of the circumcircle, the radius of the incircle, and the radii of the three excircles, are also positive integers. Moreover this geometry gives rise to integer solutions to Descartes' circle equation.

As for any right triangle, the converse of Thales' theorem says that the diameter of the circumcircle of a right triangle equals the hypotenuse; hence for Pythagorean triples, the circumdiameter is a positive integer. The circumradius may, however, be a half-integer.

For a general triangle with sides a, b, c and area K, the inradius r and exradii r_a, r_b, r_c, associated with sides a, b, c, satisfy K = rs = r_a(s − a) = r_b(s − b) = r_c(s − c), where ${\displaystyle s=\frac{1}{2}(a+b+c)}$ is the semiperimeter. When the angle opposite c is right, so that ${\displaystyle K=\frac{1}{2}ab}$ and ${\displaystyle a^2+b^2=c^2}$, it is immediate that ${\displaystyle K=s(s-c)=(s-a)(s-b),}$ which gives

${\displaystyle \begin{array}{cccc}r& =s-c=\frac{1}{2}(a+b-c),& r_a& =s-b=\frac{1}{2}(a+c-b),\\ r_b& =s-a=\frac{1}{2}(b+c-a),& r_c& =s=\frac{1}{2}(a+b+c).\end{array}}$

If (a,b,c) is a Pythagorean triple, then s, and hence all four radii, are positive integers. For a primitive triple, Euclid's formula implies that r = n(m − n), and that the radii of the excircles opposite the sides m² − n², 2mn, and the hypotenuse m² + n² are respectively m(m − n), n(m + n), and m(m + n).

Euclid's formula also implies that for primitive triples one of ${\displaystyle c\pm a}$ and ${\displaystyle c\pm b}$ is a perfect square (for both choices of sign), and the other is twice a perfect square (for both choices of sign). Hence ${\displaystyle \frac{(c\pm a)(c\pm b)}{2}}$ is a perfect square (for all Pythagorean triples and all choices of sign). The sign choices (−,−), (−,+), (+,−), (+,+) give the squares of the inradius and exradii ${\displaystyle r_b,}$ ${\displaystyle r_a,}$ and ${\displaystyle r_c,}$ respectively.

From ${\displaystyle K=s(s-c)=(s-a)(s-b)}$ it follows that that ${\displaystyle K}$ and any one of ${\displaystyle s,}$ ${\displaystyle s-c,}$ ${\displaystyle s-a,}$ or ${\displaystyle s-b,}$ determine all the side lengths of the right triangle. By rescaling, either of the ratios

${\displaystyle \begin{array}{cc}\frac{K}{s^2}& =\frac{s-c}{s}=\frac{(s-c{)}^2}{K}=\frac{n(m-n)}{m(m+n)}\\ \frac{K}{(s-a{)}^2}& =\frac{s-b}{s-a}=\frac{(s-b{)}^2}{K}=\frac{m(m-n)}{n(m+n)}\end{array}}$

uniquely determines a primitive Pythagorean triple.

The incircle and excircles of a right triangle may be translated to produce four mutually tangent circles, which is the configuration treated by Descartes' theorem. To obtain this configuration, the centers of the incircle and of the excircles tangent to the sides of lengths ${\displaystyle a}$ and ${\displaystyle b}$ move to the triangle vertices, ${\displaystyle (0,0),}$ ${\displaystyle (a,0),}$ and ${\displaystyle (0,b),}$ while the excircle of the hypotenuse, of radius ${\displaystyle s,}$ becomes the outer Soddy circle, with center at ${\displaystyle (a,b).}$ Descartes' theorem states that the bends of the four circles,

${\displaystyle \begin{array}{cccc}k& =\frac{1}{r}=\frac{1}{s-c},& k_a& =\frac{1}{r_a}=\frac{1}{s-b},\\ k_b& =\frac{1}{r_b}=\frac{1}{s-a},& k_c& =-\frac{1}{r_c}=-\frac{1}{s},\end{array}}$

which are the curvatures of the circles, apart from a minus sign factor for the outer circle, satisfy the equation

${\displaystyle {\left(k+k_a+k_b+k_c\right)}^2=2\left(k^2+k_a^2+k_b^2+k_c^2\right).}$

Multiplying each of the bends by the area ${\displaystyle K}$ gives an integer solution to this equation:

${\displaystyle K\left(k,k_a,k_b,k_c\right)=(s,s-a,s-b,-(s-c)).}$

Pythagorean triangles in a 2D lattice

A 2D lattice is a regular array of isolated points where if any one point is chosen as the Cartesian origin (0, 0), then all the other points are at (x, y) where x and y range over all positive and negative integers. Any Pythagorean triangle with triple (a, b, c) can be drawn within a 2D lattice with vertices at coordinates (0, 0), (a, 0) and (0, b). The count of lattice points lying strictly within the bounds of the triangle is given by ${\displaystyle \frac{(a-1)(b-1)-gcd(a,b)+1}{2};}$ for primitive Pythagorean triples this interior lattice count is ${\displaystyle \frac{(a-1)(b-1)}{2}.}$ The area (by Pick's theorem equal to one less than the interior lattice count plus half the boundary lattice count) equals ${\displaystyle \frac{ab}{2}}$ .

The first occurrence of two primitive Pythagorean triples sharing the same area occurs with triangles with sides (20, 21, 29), (12, 35, 37) and common area 210 (sequence A093536 in the OEIS). The first occurrence of two primitive Pythagorean triples sharing the same interior lattice count occurs with (18108, 252685, 253333), (28077, 162964, 165365) and interior lattice count 2287674594 (sequence A225760 in the OEIS). Three primitive Pythagorean triples have been found sharing the same area: (4485, 5852, 7373), (3059, 8580, 9109), (1380, 19019, 19069) with area 13123110. As yet, no set of three primitive Pythagorean triples have been found sharing the same interior lattice count.

Enumeration of primitive Pythagorean triples

By Euclid's formula all primitive Pythagorean triples can be generated from integers ${\displaystyle m}$ and ${\displaystyle n}$ with ${\displaystyle m>n>0}$, ${\displaystyle m+n}$ odd and ${\displaystyle gcd(m,n)=1.}$ Hence there is a 1 to 1 mapping of rationals (in lowest terms) to primitive Pythagorean triples where ${\displaystyle \frac{n}{m}}$ is in the interval ${\displaystyle (0,1)}$ and ${\displaystyle m+n}$ odd.

The reverse mapping from a primitive triple ${\displaystyle (a,b,c)}$ where ${\displaystyle c>b>a>0}$ to a rational ${\displaystyle \frac{n}{m}}$ is achieved by studying the two sums ${\displaystyle a+c}$ and ${\displaystyle b+c.}$ One of these sums will be a square that can be equated to ${\displaystyle (m+n{)}^2}$ and the other will be twice a square that can be equated to ${\displaystyle 2m^2.}$ It is then possible to determine the rational ${\displaystyle \frac{n}{m}.}$

In order to enumerate primitive Pythagorean triples the rational can be expressed as an ordered pair ${\displaystyle (n,m)}$ and mapped to an integer using a pairing function such as Cantor's pairing function. An example can be seen at (sequence A277557 in the OEIS). It begins

${\displaystyle 8,18,19,32,33,34,\dots }$ and gives rationals

${\displaystyle \frac{1}{2},\frac{2}{3},\frac{1}{4},\frac{3}{4},\frac{2}{5},\frac{1}{6},\dots }$ these, in turn, generate primitive triples

${\displaystyle (3,4,5),(5,12,13),(8,15,17),(7,24,25),(20,21,29),(12,35,37),\dots }$

Spinors and the modular group

Pythagorean triples can likewise be encoded into a square matrix of the form $${\displaystyle X=\left[\begin{array}{cc}c+b& a\\ a& c-b\end{array}\right].}$$ A matrix of this form is symmetric. Furthermore, the determinant of X is $${\displaystyle detX=c^2-a^2-b^2}$$ which is zero when (a,b,c) satisfies the Pythagorean theorem. If X corresponds to a Pythagorean triple, then as a matrix it must have rank 1.

Since X is symmetric and of rank 1, it follows from a result in linear algebra that there is a column vector ξ = [m n]^T such that the outer product

holds, where the T denotes the matrix transpose. Since ξ and -ξ produce the same Pythagorean triple, the vector ξ can be considered a spinor (for the Lorentz group SO(1, 2)). In abstract terms, the Euclid formula means that each primitive Pythagorean triple can be written as the outer product with itself of a spinor with integer entries, as in (1).

The modular group Γ is the set of 2×2 matrices with integer entries

$${\displaystyle A=\left[\begin{array}{cc}\alpha & \beta \\ \gamma & \delta \end{array}\right]}$$

with determinant equal to one: αδ − βγ = 1. This set forms a group, since the inverse of a matrix in Γ is again in Γ, as is the product of two matrices in Γ. The modular group acts on the collection of all integer spinors. Furthermore, the group is transitive on the collection of integer spinors with relatively prime entries. For if [m n]^T has relatively prime entries, then

$${\displaystyle \left[\begin{array}{cc}m& -v\\ n& u\end{array}\right]\left[\begin{array}{c}1\\ 0\end{array}\right]=\left[\begin{array}{c}m\\ n\end{array}\right]}$$

where u and v are selected (by the Euclidean algorithm) so that mu + nv = 1.

By acting on the spinor ξ in (1), the action of Γ goes over to an action on Pythagorean triples, provided one allows for triples with possibly negative components. Thus if A is a matrix in Γ, then

gives rise to an action on the matrix X in (1). This does not give a well-defined action on primitive triples, since it may take a primitive triple to an imprimitive one. It is convenient at this point (per Trautman 1998) to call a triple (a,b,c) standard if c > 0 and either (a,b,c) are relatively prime or (a/2,b/2,c/2) are relatively prime with a/2 odd. If the spinor [m n]^T has relatively prime entries, then the associated triple (a,b,c) determined by (1) is a standard triple. It follows that the action of the modular group is transitive on the set of standard triples.

Alternatively, restrict attention to those values of m and n for which m is odd and n is even. Let the subgroup Γ(2) of Γ be the kernel of the group homomorphism

$${\displaystyle \mathrm{\Gamma }=\mathrm{S}\mathrm{L}(2,\mathbf{Z})\to \mathrm{S}\mathrm{L}(2,{\mathbf{Z}}_2)}$$

where SL(2,Z₂) is the special linear group over the finite field Z₂ of integers modulo 2. Then Γ(2) is the group of unimodular transformations which preserve the parity of each entry. Thus if the first entry of ξ is odd and the second entry is even, then the same is true of Aξ for all A ∈ Γ(2). In fact, under the action (2), the group Γ(2) acts transitively on the collection of primitive Pythagorean triples (Alperin 2005).

The group Γ(2) is the free group whose generators are the matrices $${\displaystyle U=\left[\begin{array}{cc}1& 2\\ 0& 1\end{array}\right],L=\left[\begin{array}{cc}1& 0\\ 2& 1\end{array}\right].}$$ Consequently, every primitive Pythagorean triple can be obtained in a unique way as a product of copies of the matrices U and L.

Parent/child relationships

By a result of Berggren (1934), all primitive Pythagorean triples can be generated from the (3, 4, 5) triangle by using the three linear transformations T₁, T₂, T₃ below, where a, b, c are sides of a triple:

In other words, every primitive triple will be a "parent" to three additional primitive triples. Starting from the initial node with a = 3, b = 4, and c = 5, the operation T₁ produces the new triple

(3 − (2×4) + (2×5), (2×3) − 4 + (2×5), (2×3) − (2×4) + (3×5)) = (5, 12, 13),

and similarly T₂ and T₃ produce the triples (21, 20, 29) and (15, 8, 17).

The linear transformations T₁, T₂, and T₃ have a geometric interpretation in the language of quadratic forms. They are closely related to (but are not equal to) reflections generating the orthogonal group of x² + y² − z² over the integers.

Relation to Gaussian integers

Alternatively, Euclid's formulae can be analyzed and proved using the Gaussian integers. Gaussian integers are complex numbers of the form α = u + vi, where u and v are ordinary integers and i is the square root of negative one. The units of Gaussian integers are ±1 and ±i. The ordinary integers are called the rational integers and denoted as 'Z'. The Gaussian integers are denoted as Z[i]. The right-hand side of the Pythagorean theorem may be factored in Gaussian integers:

$${\displaystyle c^2=a^2+b^2=(a+bi)\overline{(a+bi)}=(a+bi)(a-bi).}$$

A primitive Pythagorean triple is one in which a and b are coprime, i.e., they share no prime factors in the integers. For such a triple, either a or b is even, and the other is odd; from this, it follows that c is also odd.

The two factors z := a + bi and z* := a − bi of a primitive Pythagorean triple each equal the square of a Gaussian integer. This can be proved using the property that every Gaussian integer can be factored uniquely into Gaussian primes up to units. (This unique factorization follows from the fact that, roughly speaking, a version of the Euclidean algorithm can be defined on them.) The proof has three steps. First, if a and b share no prime factors in the integers, then they also share no prime factors in the Gaussian integers. (Assume a = gu and b = gv with Gaussian integers g, u and v and g not a unit. Then u and v lie on the same line through the origin. All Gaussian integers on such a line are integer multiples of some Gaussian integer h. But then the integer gh ≠ ±1 divides both a and b.) Second, it follows that z and z* likewise share no prime factors in the Gaussian integers. For if they did, then their common divisor δ would also divide z + z* = 2a and z − z* = 2ib. Since a and b are coprime, that implies that δ divides 2 = (1 + i)(1 − i) = i(1 − i)². From the formula c² = zz*, that in turn would imply that c is even, contrary to the hypothesis of a primitive Pythagorean triple. Third, since c² is a square, every Gaussian prime in its factorization is doubled, i.e., appears an even number of times. Since z and z* share no prime factors, this doubling is also true for them. Hence, z and z* are squares.

Thus, the first factor can be written

$${\displaystyle a+bi=\epsilon {\left(m+ni\right)}^2,\epsilon \in \{\pm 1,\pm i\}.}$$

The real and imaginary parts of this equation give the two formulas:

$${\displaystyle \{\begin{array}{ll}\epsilon =+1,& a=+\left(m^2-n^2\right),b=+2mn;\\ \epsilon =-1,& a=-\left(m^2-n^2\right),b=-2mn;\\ \epsilon =+i,& a=-2mn,b=+\left(m^2-n^2\right);\\ \epsilon =-i,& a=+2mn,b=-\left(m^2-n^2\right).\end{array}}$$

For any primitive Pythagorean triple, there must be integers m and n such that these two equations are satisfied. Hence, every Pythagorean triple can be generated from some choice of these integers.

As perfect square Gaussian integers

If we consider the square of a Gaussian integer we get the following direct interpretation of Euclid's formula as representing the perfect square of a Gaussian integer.

$${\displaystyle (m+ni{)}^2=(m^2-n^2)+2mni.}$$

Using the facts that the Gaussian integers are a Euclidean domain and that for a Gaussian integer p ${\displaystyle |p{|}^2}$ is always a square it is possible to show that a Pythagorean triple corresponds to the square of a prime Gaussian integer if the hypotenuse is prime.

If the Gaussian integer is not prime then it is the product of two Gaussian integers p and q with ${\displaystyle |p{|}^2}$ and ${\displaystyle |q{|}^2}$ integers. Since magnitudes multiply in the Gaussian integers, the product must be ${\displaystyle |p||q|}$, which when squared to find a Pythagorean triple must be composite. The contrapositive completes the proof.

Distribution of triples

There are a number of results on the distribution of Pythagorean triples. In the scatter plot, a number of obvious patterns are already apparent. Whenever the legs (a,b) of a primitive triple appear in the plot, all integer multiples of (a,b) must also appear in the plot, and this property produces the appearance of lines radiating from the origin in the diagram.

Within the scatter, there are sets of parabolic patterns with a high density of points and all their foci at the origin, opening up in all four directions. Different parabolas intersect at the axes and appear to reflect off the axis with an incidence angle of 45 degrees, with a third parabola entering in a perpendicular fashion. Within this quadrant, each arc centered on the origin shows that section of the parabola that lies between its tip and its intersection with its semi-latus rectum.

These patterns can be explained as follows. If ${\displaystyle a^2/4n}$ is an integer, then (a, ${\displaystyle |n-a^2/4n|}$, ${\displaystyle n+a^2/4n}$) is a Pythagorean triple. (In fact every Pythagorean triple (a, b, c) can be written in this way with integer n, possibly after exchanging a and b, since ${\displaystyle n=(b+c)/2}$ and a and b cannot both be odd.) The Pythagorean triples thus lie on curves given by ${\displaystyle b=|n-a^2/4n|}$, that is, parabolas reflected at the a-axis, and the corresponding curves with a and b interchanged. If a is varied for a given n (i.e. on a given parabola), integer values of b occur relatively frequently if n is a square or a small multiple of a square. If several such values happen to lie close together, the corresponding parabolas approximately coincide, and the triples cluster in a narrow parabolic strip. For instance, 38² = 1444, 2 × 27² = 1458, 3 × 22² = 1452, 5 × 17² = 1445 and 10 × 12² = 1440; the corresponding parabolic strip around n ≈ 1450 is clearly visible in the scatter plot.

The angular properties described above follow immediately from the functional form of the parabolas. The parabolas are reflected at the a-axis at a = 2n, and the derivative of b with respect to a at this point is –1; hence the incidence angle is 45°. Since the clusters, like all triples, are repeated at integer multiples, the value 2n also corresponds to a cluster. The corresponding parabola intersects the b-axis at right angles at b = 2n, and hence its reflection upon interchange of a and b intersects the a-axis at right angles at a = 2n, precisely where the parabola for n is reflected at the a-axis. (The same is of course true for a and b interchanged.)

Albert Fässler and others provide insights into the significance of these parabolas in the context of conformal mappings.

Special cases and related equations

The Platonic sequence

The case n = 1 of the more general construction of Pythagorean triples has been known for a long time. Proclus, in his commentary to the 47th Proposition of the first book of Euclid's Elements, describes it as follows:

Certain methods for the discovery of triangles of this kind are handed down, one which they refer to Plato, and another to Pythagoras. (The latter) starts from odd numbers. For it makes the odd number the smaller of the sides about the right angle; then it takes the square of it, subtracts unity and makes half the difference the greater of the sides about the right angle; lastly it adds unity to this and so forms the remaining side, the hypotenuse.
...For the method of Plato argues from even numbers. It takes the given even number and makes it one of the sides about the right angle; then, bisecting this number and squaring the half, it adds unity to the square to form the hypotenuse, and subtracts unity from the square to form the other side about the right angle. ... Thus it has formed the same triangle that which was obtained by the other method.

In equation form, this becomes:

a is odd (Pythagoras, c. 540 BC):

$${\displaystyle \text{side }a:\text{side }b=\frac{a^2-1}{2}:\text{side }c=\frac{a^2+1}{2}.}$$

a is even (Plato, c. 380 BC):

$${\displaystyle \text{side }a:\text{side }b={\left(\frac{a}{2}\right)}^2-1:\text{side }c={\left(\frac{a}{2}\right)}^2+1}$$

It can be shown that all Pythagorean triples can be obtained, with appropriate rescaling, from the basic Platonic sequence (a, (a² − 1)/2 and (a² + 1)/2) by allowing a to take non-integer rational values. If a is replaced with the fraction m/n in the sequence, the result is equal to the 'standard' triple generator (2mn, m² − n²,m² + n²) after rescaling. It follows that every triple has a corresponding rational a value which can be used to generate a similar triangle (one with the same three angles and with sides in the same proportions as the original). For example, the Platonic equivalent of (56, 33, 65) is generated by a = m/n = 7/4 as (a, (a² –1)/2, (a²+1)/2) = (56/32, 33/32, 65/32). The Platonic sequence itself can be derived by following the steps for 'splitting the square' described in Diophantus II.VIII.

Arithmetic progressions of three squares

Given three positive integers ⁠${\displaystyle x<y<z}$⁠, their squares are in arithmetic progression if $${\displaystyle z^2-y^2=y^2-x^2,}$$ or equivalently if $${\displaystyle x^2+z^2=2y^2.}$$

The arithmetic progressions of three squares are in one-to-one correspondence with Pythgorean triples.

Indeed, if ⁠${\displaystyle (a,b,c)}$⁠ is a Pythagorean triple such that ⁠${\displaystyle a^2+b^2=c^2}$⁠, then $${\displaystyle (a-b{)}^2,c^2,(a+b{)}^2}$$ are three squares in arithmetic progression, since ⁠${\displaystyle (a-b{)}^2+(a+b{)}^2=2(a^2+b^2)=2c^2}$⁠.

Conversely, if ⁠${\displaystyle x^2<y^2<z^2}$⁠ are three squares in arithmetic progression, such that ⁠${\displaystyle y^2=\frac{x^2+z^2}{2}}$⁠, then $${\displaystyle \left(\frac{z-x}{2},y,\frac{z+x}{2}\right)}$$ is a Pythagorean triple. As these transformations are inverse to each other, they define the announced one-to-one correspondence.

This one-to-one correspondence implies that almost all properties of arithmetic progressions of squares are easily deduced from a corresponding property of Pythagorean triples.

In particular, Euclid's formula implies that all arithmetic progressions of three squares can be obtained as $${\displaystyle (m^2-n^2-2mn{)}^2,(m^2+m^2{)}^2,(m^2-n^2+2mn{)}^2,}$$ where ⁠${\displaystyle m}$⁠ and ⁠${\displaystyle n}$⁠ are two arbitrary positive integers of different parities.

Also, the difference ${\displaystyle z^2-y^2=y^2-x^2}$ between two consecutive squares in an arithmetic progression of squares is always a multiple of 24. This results from the divisibility relations given in § General properties, since this difference is ⁠${\displaystyle 2ab}$⁠, where ⁠${\displaystyle a}$⁠ and ⁠${\displaystyle b}$⁠ are the first elements of the associated Euclidean triple.

It results that the smallest arithmetic progression of three squares is $${\displaystyle (1,25,49).}$$

The Jacobi–Madden equation

The equation,

$${\displaystyle a^4+b^4+c^4+d^4=(a+b+c+d{)}^4}$$

is equivalent to the special Pythagorean triple,

$${\displaystyle (a^2+ab+b^2{)}^2+(c^2+cd+d^2{)}^2=((a+b{)}^2+(a+b)(c+d)+(c+d{)}^2{)}^2}$$

There is an infinite number of solutions to this equation as solving for the variables involves an elliptic curve. Small ones are,

$${\displaystyle a,b,c,d=-2634,955,1770,5400}$$

$${\displaystyle a,b,c,d=-31764,7590,27385,48150}$$

Equal sums of two squares

One way to generate solutions to ${\displaystyle a^2+b^2=c^2+d^2}$ is to parametrize a, b, c, d in terms of integers m, n, p, q as follows:

$${\displaystyle (m^2+n^2)(p^2+q^2)=(mp-nq{)}^2+(np+mq{)}^2=(mp+nq{)}^2+(np-mq{)}^2.}$$

Equal sums of two fourth powers

Given two sets of Pythagorean triples,

$${\displaystyle (a^2-b^2{)}^2+(2ab{)}^2=(a^2+b^2{)}^2}$$

$${\displaystyle (c^2-d^2{)}^2+(2cd{)}^2=(c^2+d^2{)}^2}$$

the problem of finding equal products of a non-hypotenuse side and the hypotenuse,

$${\displaystyle (a^2-b^2)(a^2+b^2)=(c^2-d^2)(c^2+d^2)}$$

is easily seen to be equivalent to the equation,

$${\displaystyle a^4-b^4=c^4-d^4}$$

and was first solved by Euler as ${\displaystyle a,b,c,d=133,59,158,134.}$ Since he showed this is a rational point in an elliptic curve, then there is an infinite number of solutions. In fact, he also found a 7th degree polynomial parameterization.

Descartes' Circle Theorem

For the case of Descartes' circle theorem where all variables are squares,

$${\displaystyle 2(a^4+b^4+c^4+d^4)=(a^2+b^2+c^2+d^2{)}^2}$$

Euler showed this is equivalent to three simultaneous Pythagorean triples,

$${\displaystyle (2ab{)}^2+(2cd{)}^2=(a^2+b^2-c^2-d^2{)}^2}$$

$${\displaystyle (2ac{)}^2+(2bd{)}^2=(a^2-b^2+c^2-d^2{)}^2}$$

$${\displaystyle (2ad{)}^2+(2bc{)}^2=(a^2-b^2-c^2+d^2{)}^2}$$

There is also an infinite number of solutions, and for the special case when ${\displaystyle a+b=c}$, then the equation simplifies to,

$${\displaystyle 4(a^2+ab+b^2)=d^2}$$

with small solutions as ${\displaystyle a,b,c,d=3,5,8,14}$ and can be solved as binary quadratic forms.

Almost-isosceles Pythagorean triples

No Pythagorean triples are isosceles, because the ratio of the hypotenuse to either other side is √2, but √2 cannot be expressed as the ratio of 2 integers.

There are, however, right-angled triangles with integral sides for which the lengths of the non-hypotenuse sides differ by one, such as,

$${\displaystyle 3^2+4^2=5^2}$$

$${\displaystyle {20}^2+{21}^2={29}^2}$$

and an infinite number of others. They can be completely parameterized as,

$${\displaystyle {\left(\frac{x-1}{2}\right)}^2+{\left(\frac{x+1}{2}\right)}^2=y^2}$$

where {x, y} are the solutions to the Pell equation ${\displaystyle x^2-2y^2=-1.}$

If a, b, c are the sides of this type of primitive Pythagorean triple then the solution to the Pell equation is given by the recursive formula

$${\displaystyle a_n=6a_{n-1}-a_{n-2}+2}$$ with ${\displaystyle a_1=3}$ and ${\displaystyle a_2=20}$ $${\displaystyle b_n=6b_{n-1}-b_{n-2}-2}$$ with ${\displaystyle b_1=4}$ and ${\displaystyle b_2=21}$ $${\displaystyle c_n=6c_{n-1}-c_{n-2}}$$ with ${\displaystyle c_1=5}$ and ⁠${\displaystyle c_2=29}$⁠.

This sequence of primitive Pythagorean triples forms the central stem (trunk) of the rooted ternary tree of primitive Pythagorean triples.

When it is the longer non-hypotenuse side and hypotenuse that differ by one, such as in

$${\displaystyle 5^2+{12}^2={13}^2}$$

$${\displaystyle 7^2+{24}^2={25}^2}$$

then the complete solution for the primitive Pythagorean triple a, b, c is

$${\displaystyle a=2m+1,b=2m^2+2m,c=2m^2+2m+1}$$

and

$${\displaystyle (2m+1{)}^2+(2m^2+2m{)}^2=(2m^2+2m+1{)}^2}$$

where integer ${\displaystyle m>0}$ is the generating parameter.

It shows that all odd numbers (greater than 1) appear in this type of almost-isosceles primitive Pythagorean triple. This sequence of primitive Pythagorean triples forms the right hand side outer stem of the rooted ternary tree of primitive Pythagorean triples.

Another property of this type of almost-isosceles primitive Pythagorean triple is that the sides are related such that $${\displaystyle a^b+b^a=Kc}$$ for some integer ${\displaystyle K}$. Or in other words ${\displaystyle a^b+b^a}$ is divisible by ${\displaystyle c}$ such as in $${\displaystyle (5^{12}+{12}^5)/13=18799189}$$.

Fibonacci numbers in Pythagorean triples

Starting with 5, every second Fibonacci number is the length of the hypotenuse of a right triangle with integer sides, or in other words, the largest number in a Pythagorean triple, obtained from the formula $${\displaystyle (F_n{F}_{n+3}{)}^2+(2F_{n+1}{F}_{n+2}{)}^2=F_{2n+3}^2.}$$ The sequence of Pythagorean triangles obtained from this formula has sides of lengths

(3,4,5), (5,12,13), (16,30,34), (39,80,89), ...

The middle side of each of these triangles is the sum of the three sides of the preceding triangle.

Generalizations

There are several ways to generalize the concept of Pythagorean triples.

Pythagorean n-tuple

The expression $${\displaystyle {\left(m_1^2-m_2^2-\dots -m_n^2\right)}^2+\sum _{k=2}^n(2m_1{m}_k{)}^2={\left(m_1^2+\dots +m_n^2\right)}^2}$$ is a Pythagorean n-tuple for any tuple of positive integers (m₁, ..., m_n) with m²
₁ > m²
₂ + ... + m²
_n. The Pythagorean n-tuple can be made primitive by dividing out by the largest common divisor of its values.

Furthermore, any primitive Pythagorean n-tuple a²
₁ + ... + a²
_n = c² can be found by this approach. Use (m₁, ..., m_n) = (c + a₁, a₂, ..., a_n) to get a Pythagorean n-tuple by the above formula and divide out by the largest common integer divisor, which is 2m₁ = 2(c + a₁). Dividing out by the largest common divisor of these (m₁, ..., m_n) values gives the same primitive Pythagorean n-tuple; and there is a one-to-one correspondence between tuples of setwise coprime positive integers (m₁, ..., m_n) satisfying m²
₁ > m²
₂ + ... + m²
_n and primitive Pythagorean n-tuples.

Examples of the relationship between setwise coprime values ${\displaystyle \overrightarrow{m}}$ and primitive Pythagorean n-tuples include:

$${\displaystyle \begin{array}{cc}\overrightarrow{m}=(1)& \leftrightarrow 1^2=1^2\\ \overrightarrow{m}=(2,1)& \leftrightarrow 3^2+4^2=5^2\\ \overrightarrow{m}=(2,1,1)& \leftrightarrow 1^2+2^2+2^2=3^2\\ \overrightarrow{m}=(3,1,1,1)& \leftrightarrow 1^2+1^2+1^2+1^2=2^2\\ \overrightarrow{m}=(5,1,1,2,3)& \leftrightarrow 1^2+1^2+1^2+2^2+3^2=4^2\\ \overrightarrow{m}=(4,1,1,1,1,2)& \leftrightarrow 1^2+1^2+1^2+1^2+1^2+2^2=3^2\\ \overrightarrow{m}=(5,1,1,1,2,2,2)& \leftrightarrow 1^2+1^2+1^2+1^2+2^2+2^2+2^2=4^2\end{array}}$$

Consecutive squares

Since the sum F(k,m) of k consecutive squares beginning with m² is given by the formula,

$${\displaystyle F(k,m)=km(k-1+m)+\frac{k(k-1)(2k-1)}{6}}$$

one may find values (k, m) so that F(k,m) is a square, such as one by Hirschhorn where the number of terms is itself a square,

$${\displaystyle m=\frac{v^4-24v^2-25}{48},k=v^2,F(m,k)=\frac{v^5+47v}{48}}$$

and v ≥ 5 is any integer not divisible by 2 or 3. For the smallest case v = 5, hence k = 25, this yields the well-known cannonball-stacking problem of Lucas,

$${\displaystyle 0^2+1^2+2^2+\cdots +{24}^2={70}^2}$$

a fact which is connected to the Leech lattice.

In addition, if in a Pythagorean n-tuple (n ≥ 4) all addends are consecutive except one, one can use the equation,

$${\displaystyle F(k,m)+p^2=(p+1{)}^2}$$

Since the second power of p cancels out, this is only linear and easily solved for as ${\displaystyle p=\frac{F(k,m)-1}{2}}$ though k, m should be chosen so that p is an integer, with a small example being k = 5, m = 1 yielding,

$${\displaystyle 1^2+2^2+3^2+4^2+5^2+{27}^2={28}^2}$$

Thus, one way of generating Pythagorean n-tuples is by using, for various x,

$${\displaystyle x^2+(x+1{)}^2+\cdots +(x+q{)}^2+p^2=(p+1{)}^2,}$$

where q = n–2 and where

$${\displaystyle p=\frac{(q+1)x^2+q(q+1)x+\frac{q(q+1)(2q+1)}{6}-1}{2}.}$$

Fermat's Last Theorem

A generalization of the concept of Pythagorean triples is the search for triples of positive integers a, b, and c, such that aⁿ + bⁿ = cⁿ, for some n strictly greater than 2. Pierre de Fermat in 1637 claimed that no such triple exists, a claim that came to be known as Fermat's Last Theorem because it took longer than any other conjecture by Fermat to be proved or disproved. The first proof was given by Andrew Wiles in 1994.

n − 1 or n nth powers summing to an nth power

Another generalization is searching for sequences of n + 1 positive integers for which the nth power of the last is the sum of the nth powers of the previous terms. The smallest sequences for known values of n are:

• n = 3: {3, 4, 5; 6}.
• n = 4: {30, 120, 272, 315; 353}
• n = 5: {19, 43, 46, 47, 67; 72}
• n = 7: {127, 258, 266, 413, 430, 439, 525; 568}
• n = 8: {90, 223, 478, 524, 748, 1088, 1190, 1324; 1409}

For the n = 3 case, in which ${\displaystyle x^3+y^3+z^3=w^3,}$ called the Fermat cubic, a general formula exists giving all solutions.

A slightly different generalization allows the sum of (k + 1) nth powers to equal the sum of (n − k) nth powers. For example:

• (n = 3): 1³ + 12³ = 9³ + 10³, made famous by Hardy's recollection of a conversation with Ramanujan about the number 1729 being the smallest number that can be expressed as a sum of two cubes in two distinct ways.

There can also exist n − 1 positive integers whose nth powers sum to an nth power (though, by Fermat's Last Theorem, not for n = 3); these are counterexamples to Euler's sum of powers conjecture. The smallest known counterexamples are

• n = 4: (95800, 217519, 414560; 422481)
• n = 5: (27, 84, 110, 133; 144)

Heronian triangle triples

A Heronian triangle is commonly defined as one with integer sides whose area is also an integer. The lengths of the sides of such a triangle form a Heronian triple (a, b, c). Every Pythagorean triple is a Heronian triple, because at least one of the legs a, b must be even in a Pythagorean triple, so the area ab/2 is an integer. Not every Heronian triple is a Pythagorean triple, however, as the example (4, 13, 15) with area 24 shows.

If (a, b, c) is a Heronian triple, so is (ka, kb, kc) where k is any positive integer; its area will be the integer that is k² times the integer area of the (a, b, c) triangle. The Heronian triple (a, b, c) is primitive provided a, b, c are setwise coprime. (With primitive Pythagorean triples the stronger statement that they are pairwise coprime also applies, but with primitive Heronian triangles the stronger statement does not always hold true, such as with (7, 15, 20).) Here are a few of the simplest primitive Heronian triples that are not Pythagorean triples:

(4, 13, 15) with area 24

(3, 25, 26) with area 36

(7, 15, 20) with area 42

(6, 25, 29) with area 60

(11, 13, 20) with area 66

(13, 14, 15) with area 84

(13, 20, 21) with area 126

By Heron's formula, the extra condition for a triple of positive integers (a, b, c) with a < b < c to be Heronian is that

(a² + b² + c²)² − 2(a⁴ + b⁴ + c⁴)

or equivalently

2(a²b² + a²c² + b²c²) − (a⁴ + b⁴ + c⁴)

be a nonzero perfect square divisible by 16.

Application to cryptography

Primitive Pythagorean triples have been used in cryptography as random sequences and for the generation of keys.