A **Pythagorean prime** is a prime number of the form 4*n* + 1. Pythagorean primes are exactly the odd prime numbers that are the sum of two squares; this characterization is Fermat's theorem on sums of two squares.

Equivalently, by the Pythagorean theorem, they are the odd prime numbers *p* for which √*p* is the length of the hypotenuse of a right triangle with integer legs, and they are also the prime numbers *p* for which *p* itself is the hypotenuse of a primitive Pythagorean triangle. For instance, the number 5 is a Pythagorean prime; √5 is the hypotenuse of a right triangle with legs 1 and 2, and 5 itself is the hypotenuse of a right triangle with legs 3 and 4.

The first few Pythagorean primes are

- 5, 13, 17, 29, 37, 41, 53, 61, 73, 89, 97, 101, 109, 113, ... (sequence A002144 in the OEIS).

By Dirichlet's theorem on arithmetic progressions, this sequence is infinite. More strongly, for each *n*, the numbers of Pythagorean and non-Pythagorean primes up to *n* are approximately equal. However, the number of Pythagorean primes up to *n* is frequently somewhat smaller than the number of non-Pythagorean primes; this phenomenon is known as Chebyshev's bias.^{[1]}
For example, the only values of *n* up to 600000 for which there are more Pythagorean than non-Pythagorean odd primes less than or equal to n are 26861 and 26862.^{[2]}

The sum of one odd square and one even square is congruent to 1 mod 4, but there exist composite numbers such as 21 that are 1 mod 4 and yet cannot be represented as sums of two squares.
Fermat's theorem on sums of two squares states that the prime numbers that can be represented as sums of two squares are exactly 2 and the odd primes congruent to 1 mod 4.^{[3]} The representation of each such number is unique, up to the ordering of the two squares.^{[4]}

By using the Pythagorean theorem, this representation can be interpreted geometrically: the Pythagorean primes are exactly the odd prime numbers *p* such that there exists a right triangle, with integer legs, whose hypotenuse has length √*p*. They are also exactly the prime numbers *p* such that there exists a right triangle with integer sides whose hypotenuse has length *p*. For, if the triangle with legs *x* and *y* has hypotenuse length √*p* (with *x* > *y*), then the triangle with legs *x*^{2} − *y*^{2} and 2*xy* has hypotenuse length *p*.^{[5]}

Another way to understand this representation as a sum of two squares involves Gaussian integers, the complex numbers whose real part and imaginary part are both integers.^{[6]}
The norm of a Gaussian integer *x* + *yi* is the number *x*^{2} + *y*^{2}.
Thus, the Pythagorean primes (and 2) occur as norms of Gaussian integers, while other primes do not.
Within the Gaussian integers, the Pythagorean primes are not considered to be prime numbers, because they can be factored as

*p*= (*x*+*yi*)(*x*−*yi*).

Similarly, their squares can be factored in a different way than their integer factorization, as

*p*^{2}= (*x*+*yi*)^{2}(*x*−*yi*)^{2}= (*x*^{2}−*y*^{2}+ 2*xyi*)(*x*^{2}−*y*^{2}− 2*xyi*).

The real and imaginary parts of the factors in these factorizations are the leg lengths of the right triangles having the given hypotenuses.

The law of quadratic reciprocity says that if *p* and *q* are distinct odd primes, at least one of which is Pythagorean, then *p* is a quadratic residue mod *q* if and only if *q* is a quadratic residue mod *p*; by contrast, if neither *p* nor *q* is Pythagorean, then *p* is a quadratic residue mod *q* if and only if *q* is **not** a quadratic residue mod *p*.^{[7]}

In the finite field **Z**/*p* with *p* a Pythagorean prime, the polynomial equation *x*^{2} = −1 has two solutions. This may be expressed by saying that −1 is a quadratic residue mod *p*. In contrast, this equation has no solution in the finite fields **Z**/*p* where *p* is an odd prime but is not Pythagorean.^{[8]}

For every Pythagorean prime *p*, there exists a Paley graph with *p* vertices, representing the numbers modulo *p*, with two numbers adjacent in the graph if and only if their difference is a quadratic residue. This definition produces the same adjacency relation regardless of the order in which the two numbers are subtracted to compute their difference, because of the property of Pythagorean primes that −1 is a quadratic residue.^{[9]}