Factorization

In mathematics, factorization (or factorisation, see English spelling differences) or factoring consists of writing a number or another mathematical object as a product of several factors, usually smaller or simpler objects of the same kind. For example, 3 × 5 is an integer factorization of 15, and (x − 2)(x + 2) is a polynomial factorization of x² − 4.

Factorization is not usually considered meaningful within number systems possessing division, such as the real or complex numbers, since any x {\displaystyle x} can be trivially written as ( x y ) × ( 1 / y ) {\displaystyle (xy)\times (1/y)} whenever y {\displaystyle y} is not zero. However, a meaningful factorization for a rational number or a rational function can be obtained by writing it in lowest terms and separately factoring its numerator and denominator.

Factorization was first considered by ancient Greek mathematicians in the case of integers. They proved the fundamental theorem of arithmetic, which asserts that every positive integer may be factored into a product of prime numbers, which cannot be further factored into integers greater than 1. Moreover, this factorization is unique up to the order of the factors. Although integer factorization is a sort of inverse to multiplication, it is much more difficult algorithmically, a fact which is exploited in the RSA cryptosystem to implement public-key cryptography.

Polynomial factorization has also been studied for centuries. In elementary algebra, factoring a polynomial reduces the problem of finding its roots to finding the roots of the factors. Polynomials with coefficients in the integers or in a field possess the unique factorization property, a version of the fundamental theorem of arithmetic with prime numbers replaced by irreducible polynomials. In particular, a univariate polynomial with complex coefficients admits a unique (up to ordering) factorization into linear polynomials: this is a version of the fundamental theorem of algebra. In this case, the factorization can be done with root-finding algorithms. The case of polynomials with integer coefficients is fundamental for computer algebra. There are efficient computer algorithms for computing (complete) factorizations within the ring of polynomials with rational number coefficients (see factorization of polynomials).

A commutative ring possessing the unique factorization property is called a unique factorization domain. There are number systems, such as certain rings of algebraic integers, which are not unique factorization domains. However, rings of algebraic integers satisfy the weaker property of Dedekind domains: ideals factor uniquely into prime ideals.

Factorization may also refer to more general decompositions of a mathematical object into the product of smaller or simpler objects. For example, every function may be factored into the composition of a surjective function with an injective function. Matrices possess many kinds of matrix factorizations. For example, every matrix has a unique LUP factorization as a product of a lower triangular matrix L with all diagonal entries equal to one, an upper triangular matrix U, and a permutation matrix P; this is a matrix formulation of Gaussian elimination.

By the fundamental theorem of arithmetic, every integer greater than 1 has a unique (up to the order of the factors) factorization into prime numbers, which are those integers which cannot be further factorized into the product of integers greater than one.

For computing the factorization of an integer n, one needs an algorithm for finding a divisor q of n or deciding that n is prime. When such a divisor is found, the repeated application of this algorithm to the factors q and n / q gives eventually the complete factorization of n.^[1]

For finding a divisor q of n, if any, it suffices to test all values of q such that 1 < q and q² ≤ n. In fact, if r is a divisor of n such that r² > n, then q = n / r is a divisor of n such that q² ≤ n.

If one tests the values of q in increasing order, the first divisor that is found is necessarily a prime number, and the cofactor r = n / q cannot have any divisor smaller than q. For getting the complete factorization, it suffices thus to continue the algorithm by searching a divisor of r that is not smaller than q and not greater than √r.

There is no need to test all values of q for applying the method. In principle, it suffices to test only prime divisors. This needs to have a table of prime numbers that may be generated for example with the sieve of Eratosthenes. As the method of factorization does essentially the same work as the sieve of Eratosthenes, it is generally more efficient to test for a divisor only those numbers for which it is not immediately clear whether they are prime or not. Typically, one may proceed by testing 2, 3, 5, and the numbers > 5, whose last digit is 1, 3, 7, 9 and the sum of digits is not a multiple of 3.

This method works well for factoring small integers, but is inefficient for larger integers. For example, Pierre de Fermat was unable to discover that the 6th Fermat number

1 + 2 2 5 = 1 + 2 32 = 4 294 967 297 {\displaystyle 1+2^{2^{5}}=1+2^{32}=4\,294\,967\,297}

is not a prime number. In fact, applying the above method would require more than 10000 divisions, for a number that has 10 decimal digits.

There are more efficient factoring algorithms. However they remain relatively inefficient, as, with the present state of the art, one cannot factorize, even with the more powerful computers, a number of 500 decimal digits that is the product of two randomly chosen prime numbers. This ensures the security of the RSA cryptosystem, which is widely used for secure internet communication.

For factoring n = 1386 into primes:

• Start with division by 2: the number is even, and n = 2 · 693. Continue with 693, and 2 as a first divisor candidate.
• 693 is odd (2 is not a divisor), but is a multiple of 3: one has 693 = 3 · 231 and n = 2 · 3 · 231. Continue with 231, and 3 as a first divisor candidate.
• 231 is also a multiple of 3: one has 231 = 3 · 77, and thus n = 2 · 3² · 77. Continue with 77, and 3 as a first divisor candidate.
• 77 is not a multiple of 3, since the sum of its digits is 14, not a multiple of 3. It is also not a multiple of 5 because its last digit is 7. The next odd divisor to be tested is 7. One has 77 = 7 · 11, and thus n = 2 · 3² · 7 · 11. This shows that 7 is prime (easy to test directly). Continue with 11, and 7 as a first divisor candidate.
• As 7² > 11, one has finished. Thus 11 is prime, and the prime factorization is

1386 = 2 · 3² · 7 · 11.

Manipulating expressions is the basis of algebra. Factorization is one of the most important methods for expression manipulation for several reasons. If one can put an equation in a factored form E⋅F = 0, then the problem of solving the equation splits into two independent (and generally easier) problems E = 0 and F = 0. When an expression can be factored, the factors are often much simpler, and may thus offer some insight on the problem. For example,

x 3 − a x 2 − b x 2 − c x 2 + a b x + a c x + b c x − a b c {\displaystyle x^{3}-ax^{2}-bx^{2}-cx^{2}+abx+acx+bcx-abc}

having 16 multiplications, 4 subtractions and 3 additions, may be factored into the much simpler expression

( x − a ) ( x − b ) ( x − c ) , {\displaystyle (x-a)(x-b)(x-c),}

with only two multiplications and three subtractions. Moreover, the factored form immediately gives roots x = a,b,c as the roots of the polynomial.

On the other hand, factorization is not always possible, and when it is possible, the factors are not always simpler. For example, x 10 − 1 {\displaystyle x^{10}-1} can be factored into two irreducible factors x − 1 {\displaystyle x-1} and x 9 + x 8 + ⋯ + x 2 + x + 1 {\displaystyle x^{9}+x^{8}+\cdots +x^{2}+x+1} .

Various methods have been developed for finding factorizations; some are described below.

Solving algebraic equations may be viewed as a problem of polynomial factorization. In fact, the fundamental theorem of algebra can be stated as follows: every polynomial in x of degree n with complex coefficients may be factorized into n linear factors x − a i , {\displaystyle x-a_{i},} for i = 1, ..., n, where the a_is are the roots of the polynomial.^[2] Even though the structure of the factorization is known in these cases, the a_is generally cannot be computed in terms of radicals (n^th roots), by the Abel–Ruffini theorem. In most cases, the best that can be done is computing approximate values of the roots with a root-finding algorithm.

The systematic use of algebraic manipulations for simplifying expressions (more specifically equations) may be dated to 9th century, with al-Khwarizmi's book The Compendious Book on Calculation by Completion and Balancing, which is titled with two such types of manipulation.

However, even for solving quadratic equations, the factoring method was not used before Harriot's work published in 1631, ten years after his death.^[3] In his book Artis Analyticae Praxis ad Aequationes Algebraicas Resolvendas, Harriot drew tables for addition, subtraction, multiplication and division of monomials, binomials, and trinomials. Then, in a second section, he set up the equation aa − ba + ca = + bc, and showed that this matches the form of multiplication he had previously provided, giving the factorization (a − b)(a + c).^[4]

The following methods apply to any expression that is a sum, or that may be transformed into a sum. Therefore, they are most often applied to polynomials, though they also may be applied when the terms of the sum are not monomials, that is, the terms of the sum are a product of variables and constants.

It may occur that all terms of a sum are products and that some factors are common to all terms. In this case, the distributive law allows factoring out this common factor. If there are several such common factors, it is preferable to divide out the greatest such common factor. Also, if there are integer coefficients, one may factor out the greatest common divisor of these coefficients.

For example,^[5] 6 x 3 y 2 + 8 x 4 y 3 − 10 x 5 y 3 = 2 x 3 y 2 ( 3 + 4 x y − 5 x 2 y ) , {\displaystyle 6x^{3}y^{2}+8x^{4}y^{3}-10x^{5}y^{3}=2x^{3}y^{2}(3+4xy-5x^{2}y),} since 2 is the greatest common divisor of 6, 8, and 10, and x 3 y 2 {\displaystyle x^{3}y^{2}} divides all terms.

Grouping terms may allow using other methods for getting a factorization.

For example, to factor 4 x 2 + 20 x + 3 x y + 15 y , {\displaystyle 4x^{2}+20x+3xy+15y,} one may remark that the first two terms have a common factor x, and the last two terms have the common factor y. Thus 4 x 2 + 20 x + 3 x y + 15 y = ( 4 x 2 + 20 x ) + ( 3 x y + 15 y ) = 4 x ( x + 5 ) + 3 y ( x + 5 ) . {\displaystyle 4x^{2}+20x+3xy+15y=(4x^{2}+20x)+(3xy+15y)=4x(x+5)+3y(x+5).} Then a simple inspection shows the common factor x + 5, leading to the factorization 4 x 2 + 20 x + 3 x y + 15 y = ( 4 x + 3 y ) ( x + 5 ) . {\displaystyle 4x^{2}+20x+3xy+15y=(4x+3y)(x+5).}

In general, this works for sums of 4 terms that have been obtained as the product of two binomials. Although not frequently, this may work also for more complicated examples.

Sometimes, some term grouping reveals part of a recognizable pattern. It is then useful to add and subtract terms to complete the pattern.

A typical use of this is the completing the square method for getting the quadratic formula.

Another example is the factorization of x 4 + 1. {\displaystyle x^{4}+1.} If one introduces the non-real square root of –1, commonly denoted i, then one has a difference of squares x 4 + 1 = ( x 2 + i ) ( x 2 − i ) . {\displaystyle x^{4}+1=(x^{2}+i)(x^{2}-i).} However, one may also want a factorization with real number coefficients. By adding and subtracting 2 x 2 , {\displaystyle 2x^{2},} and grouping three terms together, one may recognize the square of a binomial: x 4 + 1 = ( x 4 + 2 x 2 + 1 ) − 2 x 2 = ( x 2 + 1 ) 2 − ( x 2 ) 2 = ( x 2 + x 2 + 1 ) ( x 2 − x 2 + 1 ) . {\displaystyle x^{4}+1=(x^{4}+2x^{2}+1)-2x^{2}=(x^{2}+1)^{2}-\left(x{\sqrt {2}}\right)^{2}=\left(x^{2}+x{\sqrt {2}}+1\right)\left(x^{2}-x{\sqrt {2}}+1\right).} Subtracting and adding 2 x 2 {\displaystyle 2x^{2}} also yields the factorization: x 4 + 1 = ( x 4 − 2 x 2 + 1 ) + 2 x 2 = ( x 2 − 1 ) 2 + ( x 2 ) 2 = ( x 2 + x − 2 − 1 ) ( x 2 − x − 2 − 1 ) . {\displaystyle x^{4}+1=(x^{4}-2x^{2}+1)+2x^{2}=(x^{2}-1)^{2}+\left(x{\sqrt {2}}\right)^{2}=\left(x^{2}+x{\sqrt {-2}}-1\right)\left(x^{2}-x{\sqrt {-2}}-1\right).} These factorizations work not only over the complex numbers, but also over any field, where either –1, 2 or –2 is a square. In a finite field, the product of two non-squares is a square; this implies that the polynomial x 4 + 1 , {\displaystyle x^{4}+1,} which is irreducible over the integers, is reducible modulo every prime number. For example, x 4 + 1 ≡ ( x + 1 ) 4 ( mod 2 ) ; {\displaystyle x^{4}+1\equiv (x+1)^{4}{\pmod {2}};} x 4 + 1 ≡ ( x 2 + x − 1 ) ( x 2 − x − 1 ) ( mod 3 ) , {\displaystyle x^{4}+1\equiv (x^{2}+x-1)(x^{2}-x-1){\pmod {3}},} since 1 2 ≡ − 2 ( mod 3 ) ; {\displaystyle 1^{2}\equiv -2{\pmod {3}};} x 4 + 1 ≡ ( x 2 + 2 ) ( x 2 − 2 ) ( mod 5 ) , {\displaystyle x^{4}+1\equiv (x^{2}+2)(x^{2}-2){\pmod {5}},} since 2 2 ≡ − 1 ( mod 5 ) ; {\displaystyle 2^{2}\equiv -1{\pmod {5}};} x 4 + 1 ≡ ( x 2 + 3 x + 1 ) ( x 2 − 3 x + 1 ) ( mod 7 ) , {\displaystyle x^{4}+1\equiv (x^{2}+3x+1)(x^{2}-3x+1){\pmod {7}},} since 3 2 ≡ 2 ( mod 7 ) . {\displaystyle 3^{2}\equiv 2{\pmod {7}}.}

Many identities provide an equality between a sum and a product. The above methods may be used for letting the sum side of some identity appear in an expression, which may therefore be replaced by a product.

Below are identities whose left-hand sides are commonly used as patterns (this means that the variables E and F that appear in these identities may represent any subexpression of the expression that has to be factorized).^[6]

• Difference of two squares

E 2 − F 2 = ( E + F ) ( E − F ) {\displaystyle E^{2}-F^{2}=(E+F)(E-F)}

For example,

a 2 + 2 a b + b 2 − x 2 + 2 x y − y 2 = ( a 2 + 2 a b + b 2 ) − ( x 2 − 2 x y + y 2 ) = ( a + b ) 2 − ( x − y ) 2 = ( a + b + x − y ) ( a + b − x + y ) . {\displaystyle {\begin{aligned}a^{2}+&2ab+b^{2}-x^{2}+2xy-y^{2}\\&=(a^{2}+2ab+b^{2})-(x^{2}-2xy+y^{2})\\&=(a+b)^{2}-(x-y)^{2}\\&=(a+b+x-y)(a+b-x+y).\end{aligned}}}

• Sum/difference of two cubes

E 3 + F 3 = ( E + F ) ( E 2 − E F + F 2 ) {\displaystyle E^{3}+F^{3}=(E+F)(E^{2}-EF+F^{2})}

E 3 − F 3 = ( E − F ) ( E 2 + E F + F 2 ) {\displaystyle E^{3}-F^{3}=(E-F)(E^{2}+EF+F^{2})}

• Cauchy identity

a 3 + b 3 + 3 a b ( a + b ) = ( a + b ) 3 {\displaystyle a^{3}+b^{3}+3ab(a+b)=(a+b)^{3}}

a 3 − b 3 − 3 a b ( a − b ) = ( a − b ) 3 {\displaystyle a^{3}-b^{3}-3ab(a-b)=(a-b)^{3}}

• Difference of two fourth powers

E 4 − F 4 = ( E 2 + F 2 ) ( E 2 − F 2 ) = ( E 2 + F 2 ) ( E + F ) ( E − F ) {\displaystyle {\begin{aligned}E^{4}-F^{4}&=(E^{2}+F^{2})(E^{2}-F^{2})\\&=(E^{2}+F^{2})(E+F)(E-F)\end{aligned}}}

• Sum/difference of two nth powers

In the following identities, the factors may often be further factorized:

• Difference, even exponent

E 2 n − F 2 n = ( E n + F n ) ( E n − F n ) {\displaystyle E^{2n}-F^{2n}=(E^{n}+F^{n})(E^{n}-F^{n})}

• Difference, even or odd exponent

E n − F n = ( E − F ) ( E n − 1 + E n − 2 F + E n − 3 F 2 + ⋯ + E F n − 2 + F n − 1 ) = ( E − F ) ∑ k = 0 n − 1 E n − 1 − k F k {\displaystyle E^{n}-F^{n}=(E-F)(E^{n-1}+E^{n-2}F+E^{n-3}F^{2}+\cdots +EF^{n-2}+F^{n-1})=(E-F)\sum _{k=0}^{n-1}{E^{n-1-k}F^{k}}}

This is an example showing that the factors may be much larger than the sum that is factorized.

• Sum, odd exponent

E n + F n = ( E + F ) ( E n − 1 − E n − 2 F + E n − 3 F 2 − ⋯ − E F n − 2 + F n − 1 ) = ( E + F ) ∑ k = 0 n − 1 ( − 1 ) k E n − 1 − k F k {\displaystyle E^{n}+F^{n}=(E+F)(E^{n-1}-E^{n-2}F+E^{n-3}F^{2}-\cdots -EF^{n-2}+F^{n-1})=(E+F)\sum _{k=0}^{n-1}(-1)^{k}{E^{n-1-k}F^{k}}}

(obtained by changing F by –F in the preceding formula)

• Sum, even exponent

If the exponent is a power of two then the expression cannot, in general, be factorized without introducing complex numbers (if E and F contain complex numbers, this may be not the case). If n has an odd divisor, that is if n = pq with p odd, one may use the preceding formula (in "Sum, odd exponent") applied to ( E q ) p + ( F q ) p . {\displaystyle (E^{q})^{p}+(F^{q})^{p}.}

• Trinomials and cubic formulas

x 2 + y 2 + z 2 + 2 ( x y + y z + x z ) = ( x + y + z ) 2 x 3 + y 3 + z 3 − 3 x y z = ( x + y + z ) ( x 2 + y 2 + z 2 − x y − x z − y z ) x 3 + y 3 + z 3 + 3 x 2 ( y + z ) + 3 y 2 ( x + z ) + 3 z 2 ( x + y ) + 6 x y z = ( x + y + z ) 3 x 3 + y 3 + z 3 + 3 ( x + y ) ( y + z ) ( x + z ) = ( x + y + z ) 3 {\displaystyle {\begin{aligned}&x^{2}+y^{2}+z^{2}+2(xy+yz+xz)=(x+y+z)^{2}\\&x^{3}+y^{3}+z^{3}-3xyz=(x+y+z)(x^{2}+y^{2}+z^{2}-xy-xz-yz)\\&x^{3}+y^{3}+z^{3}+3x^{2}(y+z)+3y^{2}(x+z)+3z^{2}(x+y)+6xyz=(x+y+z)^{3}\\&x^{3}+y^{3}+z^{3}+3(x+y)(y+z)(x+z)=(x+y+z)^{3}\\\end{aligned}}}

• Argand identity

x 4 + x 2 y 2 + y 4 = ( x 2 + x y + y 2 ) ( x 2 − x y + y 2 ) {\displaystyle x^{4}+x^{2}y^{2}+y^{4}=(x^{2}+xy+y^{2})(x^{2}-xy+y^{2})}

x 4 + x 2 + 1 = ( x 2 + x + 1 ) ( x 2 − x + 1 ) {\displaystyle x^{4}+x^{2}+1=(x^{2}+x+1)(x^{2}-x+1)}

• Binomial expansions

The binomial theorem supplies patterns that can easily be recognized from the integers that appear in them

In low degree:

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

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

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

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

More generally, the coefficients of the expanded forms of ( a + b ) n {\displaystyle (a+b)^{n}} and ( a − b ) n {\displaystyle (a-b)^{n}} are the binomial coefficients, that appear in the nth row of Pascal's triangle.

The nth roots of unity are the complex numbers each of which is a root of the polynomial x n − 1. {\displaystyle x^{n}-1.} They are thus the numbers e 2 i k π / n = cos ⁡ 2 π k n + i sin ⁡ 2 π k n {\displaystyle e^{2ik\pi /n}=\cos {\tfrac {2\pi k}{n}}+i\sin {\tfrac {2\pi k}{n}}} for k = 0 , … , n − 1. {\displaystyle k=0,\ldots ,n-1.}

It follows that for any two expressions E and F, one has: E n − F n = ( E − F ) ∏ k = 1 n − 1 ( E − F e 2 i k π / n ) {\displaystyle E^{n}-F^{n}=(E-F)\prod _{k=1}^{n-1}\left(E-Fe^{2ik\pi /n}\right)} E n + F n = ∏ k = 0 n − 1 ( E − F e ( 2 k + 1 ) i π / n ) if  n  is even {\displaystyle E^{n}+F^{n}=\prod _{k=0}^{n-1}\left(E-Fe^{(2k+1)i\pi /n}\right)\qquad {\text{if }}n{\text{ is even}}} E n + F n = ( E + F ) ∏ k = 1 n − 1 ( E + F e 2 i k π / n ) if  n  is odd {\displaystyle E^{n}+F^{n}=(E+F)\prod _{k=1}^{n-1}\left(E+Fe^{2ik\pi /n}\right)\qquad {\text{if }}n{\text{ is odd}}}

If E and F are real expressions, and one wants real factors, one has to replace every pair of complex conjugate factors by its product. As the complex conjugate of e i α {\displaystyle e^{i\alpha }} is e − i α , {\displaystyle e^{-i\alpha },} and ( a − b e i α ) ( a − b e − i α ) = a 2 − a b ( e i α + e − i α ) + b 2 e i α e − i α = a 2 − 2 a b cos α + b 2 , {\displaystyle \left(a-be^{i\alpha }\right)\left(a-be^{-i\alpha }\right)=a^{2}-ab\left(e^{i\alpha }+e^{-i\alpha }\right)+b^{2}e^{i\alpha }e^{-i\alpha }=a^{2}-2ab\cos \,\alpha +b^{2},} one has the following real factorizations (one passes from one to the other by changing k into n − k or n + 1 − k, and applying the usual trigonometric formulas: E 2 n − F 2 n = ( E − F ) ( E + F ) ∏ k = 1 n − 1 ( E 2 − 2 E F cos k π n + F 2 ) = ( E − F ) ( E + F ) ∏ k = 1 n − 1 ( E 2 + 2 E F cos k π n + F 2 ) {\displaystyle {\begin{aligned}E^{2n}-F^{2n}&=(E-F)(E+F)\prod _{k=1}^{n-1}\left(E^{2}-2EF\cos \,{\tfrac {k\pi }{n}}+F^{2}\right)\\&=(E-F)(E+F)\prod _{k=1}^{n-1}\left(E^{2}+2EF\cos \,{\tfrac {k\pi }{n}}+F^{2}\right)\end{aligned}}} E 2 n + F 2 n = ∏ k = 1 n ( E 2 + 2 E F cos ( 2 k − 1 ) π 2 n + F 2 ) = ∏ k = 1 n ( E 2 − 2 E F cos ( 2 k − 1 ) π 2 n + F 2 ) {\displaystyle {\begin{aligned}E^{2n}+F^{2n}&=\prod _{k=1}^{n}\left(E^{2}+2EF\cos \,{\tfrac {(2k-1)\pi }{2n}}+F^{2}\right)\\&=\prod _{k=1}^{n}\left(E^{2}-2EF\cos \,{\tfrac {(2k-1)\pi }{2n}}+F^{2}\right)\end{aligned}}}

The cosines that appear in these factorizations are algebraic numbers, and may be expressed in terms of radicals (this is possible because their Galois group is cyclic); however, these radical expressions are too complicated to be used, except for low values of n. For example, a 4 + b 4 = ( a 2 − 2 a b + b 2 ) ( a 2 + 2 a b + b 2 ) . {\displaystyle a^{4}+b^{4}=(a^{2}-{\sqrt {2}}ab+b^{2})(a^{2}+{\sqrt {2}}ab+b^{2}).} a 5 − b 5 = ( a − b ) ( a 2 + 1 − 5 2 a b + b 2 ) ( a 2 + 1 + 5 2 a b + b 2 ) , {\displaystyle a^{5}-b^{5}=(a-b)\left(a^{2}+{\frac {1-{\sqrt {5}}}{2}}ab+b^{2}\right)\left(a^{2}+{\frac {1+{\sqrt {5}}}{2}}ab+b^{2}\right),} a 5 + b 5 = ( a + b ) ( a 2 − 1 − 5 2 a b + b 2 ) ( a 2 − 1 + 5 2 a b + b 2 ) , {\displaystyle a^{5}+b^{5}=(a+b)\left(a^{2}-{\frac {1-{\sqrt {5}}}{2}}ab+b^{2}\right)\left(a^{2}-{\frac {1+{\sqrt {5}}}{2}}ab+b^{2}\right),}

Often one wants a factorization with rational coefficients. Such a factorization involves cyclotomic polynomials. To express rational factorizations of sums and differences or powers, we need a notation for the homogenization of a polynomial: if P ( x ) = a 0 x n + a i x n − 1 + ⋯ + a n , {\displaystyle P(x)=a_{0}x^{n}+a_{i}x^{n-1}+\cdots +a_{n},} its homogenization is the bivariate polynomial P ¯ ( x , y ) = a 0 x n + a i x n − 1 y + ⋯ + a n y n . {\displaystyle {\overline {P}}(x,y)=a_{0}x^{n}+a_{i}x^{n-1}y+\cdots +a_{n}y^{n}.} Then, one has E n − F n = ∏ k ∣ n Q ¯ n ( E , F ) , {\displaystyle E^{n}-F^{n}=\prod _{k\mid n}{\overline {Q}}_{n}(E,F),} E n + F n = ∏ k ∣ 2 n , k ∤ n Q ¯ n ( E , F ) , {\displaystyle E^{n}+F^{n}=\prod _{k\mid 2n,k\not \mid n}{\overline {Q}}_{n}(E,F),} where the products are taken over all divisors of n, or all divisors of 2n that do not divide n, and Q n ( x ) {\displaystyle Q_{n}(x)} is the nth cyclotomic polynomial.

For example, a 6 − b 6 = Q ¯ 1 ( a , b ) Q ¯ 2 ( a , b ) Q ¯ 3 ( a , b ) Q ¯ 6 ( a , b ) = ( a − b ) ( a + b ) ( a 2 − a b + b 2 ) ( a 2 + a b + b 2 ) , {\displaystyle a^{6}-b^{6}={\overline {Q}}_{1}(a,b){\overline {Q}}_{2}(a,b){\overline {Q}}_{3}(a,b){\overline {Q}}_{6}(a,b)=(a-b)(a+b)(a^{2}-ab+b^{2})(a^{2}+ab+b^{2}),} a 6 + b 6 = Q ¯ 4 ( a , b ) Q ¯ 12 ( a , b ) = ( a 2 + b 2 ) ( a 4 − a 2 b 2 + b 4 ) , {\displaystyle a^{6}+b^{6}={\overline {Q}}_{4}(a,b){\overline {Q}}_{12}(a,b)=(a^{2}+b^{2})(a^{4}-a^{2}b^{2}+b^{4}),} since the divisors of 6 are 1, 2, 3, 6, and the divisors of 12 that do not divide 6 are 4 and 12.

For polynomials, factorization is strongly related with the problem of solving algebraic equations. An algebraic equation has the form

P ( x )   = def   a 0 x n + a 1 x n − 1 + ⋯ + a n = 0 , {\displaystyle P(x)\ \,{\stackrel {\text{def}}{=}}\ \,a_{0}x^{n}+a_{1}x^{n-1}+\cdots +a_{n}=0,}

where P(x) is a polynomial in x with a 0 ≠ 0. {\displaystyle a_{0}\neq 0.} A solution of this equation (also called a root of the polynomial) is a value r of x such that

P ( r ) = 0. {\displaystyle P(r)=0.}

If P ( x ) = Q ( x ) R ( x ) {\displaystyle P(x)=Q(x)R(x)} is a factorization of P(x) = 0 as a product of two polynomials, then the roots of P(x) are the union of the roots of Q(x) and the roots of R(x). Thus solving P(x) = 0 is reduced to the simpler problems of solving Q(x) = 0 and R(x) = 0.

Conversely, the factor theorem asserts that, if r is a root of P(x) = 0, then P(x) may be factored as

P ( x ) = ( x − r ) Q ( x ) , {\displaystyle P(x)=(x-r)Q(x),}

where Q(x) is the quotient of Euclidean division of P(x) = 0 by the linear (degree one) factor x − r.

If the coefficients of P(x) are real or complex numbers, the fundamental theorem of algebra asserts that P(x) has a real or complex root. Using the factor theorem recursively, it results that

P ( x ) = a 0 ( x − r 1 ) ⋯ ( x − r n ) , {\displaystyle P(x)=a_{0}(x-r_{1})\cdots (x-r_{n}),}

where r 1 , … , r n {\displaystyle r_{1},\ldots ,r_{n}} are the real or complex roots of P, with some of them possibly repeated. This complete factorization is unique up to the order of the factors.

If the coefficients of P(x) are real, one generally wants a factorization where factors have real coefficients. In this case, the complete factorization may have some quadratic (degree two) factors. This factorization may easily be deduced from the above complete factorization. In fact, if r = a + ib is a non-real root of P(x), then its complex conjugate s = a − ib is also a root of P(x). So, the product

( x − r ) ( x − s ) = x 2 − ( r + s ) x + r s = x 2 − 2 a x + a 2 + b 2 {\displaystyle (x-r)(x-s)=x^{2}-(r+s)x+rs=x^{2}-2ax+a^{2}+b^{2}}

is a factor of P(x) with real coefficients. Repeating this for all non-real factors gives a factorization with linear or quadratic real factors.

For computing these real or complex factorizations, one needs the roots of the polynomial, which may not be computed exactly, and only approximated using root-finding algorithms.

In practice, most algebraic equations of interest have integer or rational coefficients, and one may want a factorization with factors of the same kind. The fundamental theorem of arithmetic may be generalized to this case, stating that polynomials with integer or rational coefficients have the unique factorization property. More precisely, every polynomial with rational coefficients may be factorized in a product

P ( x ) = q P 1 ( x ) ⋯ P k ( x ) , {\displaystyle P(x)=q\,P_{1}(x)\cdots P_{k}(x),}

where q is a rational number and P 1 ( x ) , … , P k ( x ) {\displaystyle P_{1}(x),\ldots ,P_{k}(x)} are non-constant polynomials with integer coefficients that are irreducible and primitive; this means that none of the P i ( x ) {\displaystyle P_{i}(x)} may be written as the product two polynomials (with integer coefficients) that are neither 1 nor −1 (integers are considered as polynomials of degree zero). Moreover, this factorization is unique up to the order of the factors and the signs of the factors.

There are efficient algorithms for computing this factorization, which are implemented in most computer algebra systems. See Factorization of polynomials. Unfortunately, these algorithms are too complicated to use for paper-and-pencil computations. Besides the heuristics above, only a few methods are suitable for hand computations, which generally work only for polynomials of low degree, with few nonzero coefficients. The main such methods are described in next subsections.

Every polynomial with rational coefficients, may be factorized, in a unique way, as the product of a rational number and a polynomial with integer coefficients, which is primitive (that is, the greatest common divisor of the coefficients is 1), and has a positive leading coefficient (coefficient of the term of the highest degree). For example:

− 10 x 2 + 5 x + 5 = ( − 5 ) ⋅ ( 2 x 2 − x − 1 ) {\displaystyle -10x^{2}+5x+5=(-5)\cdot (2x^{2}-x-1)}

1 3 x 5 + 7 2 x 2 + 2 x + 1 = 1 6 ( 2 x 5 + 21 x 2 + 12 x + 6 ) {\displaystyle {\frac {1}{3}}x^{5}+{\frac {7}{2}}x^{2}+2x+1={\frac {1}{6}}(2x^{5}+21x^{2}+12x+6)}

In this factorization, the rational number is called the content, and the primitive polynomial is the primitive part. The computation of this factorization may be done as follows: firstly, reduce all coefficients to a common denominator, for getting the quotient by an integer q of a polynomial with integer coefficients. Then one divides out the greater common divisor p of the coefficients of this polynomial for getting the primitive part, the content being p / q . {\displaystyle p/q.} Finally, if needed, one changes the signs of p and all coefficients of the primitive part.

This factorization may produce a result that is larger than the original polynomial (typically when there are many coprime denominators), but, even when this is the case, the primitive part is generally easier to manipulate for further factorization.

The factor theorem states that, if r is a root of a polynomial

P ( x ) = a 0 x n + a 1 x n − 1 + ⋯ + a n − 1 x + a n , {\displaystyle P(x)=a_{0}x^{n}+a_{1}x^{n-1}+\cdots +a_{n-1}x+a_{n},}

meaning P(r) = 0, then there is a factorization

P ( x ) = ( x − r ) Q ( x ) , {\displaystyle P(x)=(x-r)Q(x),}

where

Q ( x ) = b 0 x n − 1 + ⋯ + b n − 2 x + b n − 1 , {\displaystyle Q(x)=b_{0}x^{n-1}+\cdots +b_{n-2}x+b_{n-1},}

with a 0 = b 0 {\displaystyle a_{0}=b_{0}} . Then polynomial long division or synthetic division give:

b i = a 0 r i + ⋯ + a i − 1 r + a i    for    i = 1 , … , n − 1. {\displaystyle b_{i}=a_{0}r^{i}+\cdots +a_{i-1}r+a_{i}\ {\text{ for }}\ i=1,\ldots ,n{-}1.}

This may be useful when one knows or can guess a root of the polynomial.

For example, for P ( x ) = x 3 − 3 x + 2 , {\displaystyle P(x)=x^{3}-3x+2,} one may easily see that the sum of its coefficients is 0, so r = 1 is a root. As r + 0 = 1, and r 2 + 0 r − 3 = − 2 , {\displaystyle r^{2}+0r-3=-2,} one has

x 3 − 3 x + 2 = ( x − 1 ) ( x 2 + x − 2 ) . {\displaystyle x^{3}-3x+2=(x-1)(x^{2}+x-2).}

For polynomials with rational number coefficients, one may search for roots which are rational numbers. Primitive part-content factorization (see above) reduces the problem of searching for rational roots to the case of polynomials with integer coefficients having no non-trivial common divisor.

If x = p q {\displaystyle x={\tfrac {p}{q}}} is a rational root of such a polynomial

P ( x ) = a 0 x n + a 1 x n − 1 + ⋯ + a n − 1 x + a n , {\displaystyle P(x)=a_{0}x^{n}+a_{1}x^{n-1}+\cdots +a_{n-1}x+a_{n},}

the factor theorem shows that one has a factorization

P ( x ) = ( q x − p ) Q ( x ) , {\displaystyle P(x)=(qx-p)Q(x),}

where both factors have integer coefficients (the fact that Q has integer coefficients results from the above formula for the quotient of P(x) by x − p / q {\displaystyle x-p/q} ).

Comparing the coefficients of degree n and the constant coefficients in the above equality shows that, if p q {\displaystyle {\tfrac {p}{q}}} is a rational root in reduced form, then q is a divisor of a 0 , {\displaystyle a_{0},} and p is a divisor of a n . {\displaystyle a_{n}.} Therefore, there is a finite number of possibilities for p and q, which can be systematically examined.^[7]

For example, if the polynomial

P ( x ) = 2 x 3 − 7 x 2 + 10 x − 6 {\displaystyle P(x)=2x^{3}-7x^{2}+10x-6}

has a rational root p q {\displaystyle {\tfrac {p}{q}}} with q > 0, then p must divide 6; that is p ∈ { ± 1 , ± 2 , ± 3 , ± 6 } , {\displaystyle p\in \{\pm 1,\pm 2,\pm 3,\pm 6\},} and q must divide 2, that is q ∈ { 1 , 2 } . {\displaystyle q\in \{1,2\}.} Moreover, if x < 0, all terms of the polynomial are negative, and, therefore, a root cannot be negative. That is, one must have

p q ∈ { 1 , 2 , 3 , 6 , 1 2 , 3 2 } . {\displaystyle {\tfrac {p}{q}}\in \{1,2,3,6,{\tfrac {1}{2}},{\tfrac {3}{2}}\}.}

A direct computation shows that only 3 2 {\displaystyle {\tfrac {3}{2}}} is a root, so there can be no other rational root. Applying the factor theorem leads finally to the factorization 2 x 3 − 7 x 2 + 10 x − 6 = ( 2 x − 3 ) ( x 2 − 2 x + 2 ) . {\displaystyle 2x^{3}-7x^{2}+10x-6=(2x-3)(x^{2}-2x+2).}

The above method may be adapted for quadratic polynomials, leading to the ac method of factorization.^[8]

Consider the quadratic polynomial

P ( x ) = a x 2 + b x + c {\displaystyle P(x)=ax^{2}+bx+c}

with integer coefficients. If it has a rational root, its denominator must divide a evenly and it may be written as a possibly reducible fraction r 1 = r a . {\displaystyle r_{1}={\tfrac {r}{a}}.} By Vieta's formulas, the other root r 2 {\displaystyle r_{2}} is

r 2 = − b a − r 1 = − b a − r a = − b + r a = s a , {\displaystyle r_{2}=-{\frac {b}{a}}-r_{1}=-{\frac {b}{a}}-{\frac {r}{a}}=-{\frac {b+r}{a}}={\frac {s}{a}},}

with s = − ( b + r ) . {\displaystyle s=-(b+r).} Thus the second root is also rational, and Vieta's second formula r 1 r 2 = c a {\displaystyle r_{1}r_{2}={\frac {c}{a}}} gives

s a r a = c a , {\displaystyle {\frac {s}{a}}{\frac {r}{a}}={\frac {c}{a}},}

that is

r s = a c and r + s = − b . {\displaystyle rs=ac\quad {\text{and}}\quad r+s=-b.}

Checking all pairs of integers whose product is ac gives the rational roots, if any.

In summary, if a x 2 + b x + c {\displaystyle ax^{2}+bx+c} has rational roots there are integers r and s such r s = a c {\displaystyle rs=ac} and r + s = − b {\displaystyle r+s=-b} (a finite number of cases to test), and the roots are r a {\displaystyle {\tfrac {r}{a}}} and s a . {\displaystyle {\tfrac {s}{a}}.} In other words, one has the factorization

a ( a x 2 + b x + c ) = ( a x − r ) ( a x − s ) . {\displaystyle a(ax^{2}+bx+c)=(ax-r)(ax-s).}

For example, let consider the quadratic polynomial

6 x 2 + 13 x + 6. {\displaystyle 6x^{2}+13x+6.}

Inspection of the factors of ac = 36 leads to 4 + 9 = 13 = b, giving the two roots

r 1 = − 4 6 = − 2 3 and r 2 = − 9 6 = − 3 2 , {\displaystyle r_{1}=-{\frac {4}{6}}=-{\frac {2}{3}}\quad {\text{and}}\quad r_{2}=-{\frac {9}{6}}=-{\frac {3}{2}},}

and the factorization

6 x 2 + 13 x + 6 = 6 ( x + 2 3 ) ( x + 3 2 ) = ( 3 x + 2 ) ( 2 x + 3 ) . {\displaystyle 6x^{2}+13x+6=6(x+{\tfrac {2}{3}})(x+{\tfrac {3}{2}})=(3x+2)(2x+3).}

Any univariate quadratic polynomial a x 2 + b x + c {\displaystyle ax^{2}+bx+c} can be factored using the quadratic formula:

a x 2 + b x + c = a ( x − α ) ( x − β ) = a ( x − − b + b 2 − 4 a c 2 a ) ( x − − b − b 2 − 4 a c 2 a ) , {\displaystyle ax^{2}+bx+c=a(x-\alpha )(x-\beta )=a\left(x-{\frac {-b+{\sqrt {b^{2}-4ac}}}{2a}}\right)\left(x-{\frac {-b-{\sqrt {b^{2}-4ac}}}{2a}}\right),}

where α {\displaystyle \alpha } and β {\displaystyle \beta } are the two roots of the polynomial.

If a, b, c are all real, the factors are real if and only if the discriminant b 2 − 4 a c {\displaystyle b^{2}-4ac} is non-negative. Otherwise, the quadratic polynomial cannot be factorized into non-constant real factors.

The quadratic formula is valid when the coefficients belong to any field of characteristic different from two, and, in particular, for coefficients in a finite field with an odd number of elements.^[9]

There are also formulas for roots of cubic and quartic polynomials, which are, in general, too complicated for practical use. The Abel–Ruffini theorem shows that there are no general root formulas in terms of radicals for polynomials of degree five or higher.

It may occur that one knows some relationship between the roots of a polynomial and its coefficients. Using this knowledge may help factoring the polynomial and finding its roots. Galois theory is based on a systematic study of the relations between roots and coefficients, that include Vieta's formulas.

Here, we consider the simpler case where two roots x 1 {\displaystyle x_{1}} and x 2 {\displaystyle x_{2}} of a polynomial P ( x ) {\displaystyle P(x)} satisfy the relation

x 2 = Q ( x 1 ) , {\displaystyle x_{2}=Q(x_{1}),}

where Q is a polynomial.

This implies that x 1 {\displaystyle x_{1}} is a common root of P ( Q ( x ) ) {\displaystyle P(Q(x))} and P ( x ) . {\displaystyle P(x).} It is therefore a root of the greatest common divisor of these two polynomials. It follows that this greatest common divisor is a non constant factor of P ( x ) . {\displaystyle P(x).} Euclidean algorithm for polynomials allows computing this greatest common factor.

For example,^[10] if one know or guess that: P ( x ) = x 3 − 5 x 2 − 16 x + 80 {\displaystyle P(x)=x^{3}-5x^{2}-16x+80} has two roots that sum to zero, one may apply Euclidean algorithm to P ( x ) {\displaystyle P(x)} and P ( − x ) . {\displaystyle P(-x).} The first division step consists in adding P ( x ) {\displaystyle P(x)} to P ( − x ) , {\displaystyle P(-x),} giving the remainder of

− 10 ( x 2 − 16 ) . {\displaystyle -10(x^{2}-16).}

Then, dividing P ( x ) {\displaystyle P(x)} by x 2 − 16 {\displaystyle x^{2}-16} gives zero as a new remainder, and x − 5 as a quotient, leading to the complete factorization

x 3 − 5 x 2 − 16 x + 80 = ( x − 5 ) ( x − 4 ) ( x + 4 ) . {\displaystyle x^{3}-5x^{2}-16x+80=(x-5)(x-4)(x+4).}

The integers and the polynomials over a field share the property of unique factorization, that is, every nonzero element may be factored into a product of an invertible element (a unit, ±1 in the case of integers) and a product of irreducible elements (prime numbers, in the case of integers), and this factorization is unique up to rearranging the factors and shifting units among the factors. Integral domains which share this property are called unique factorization domains (UFD).

Greatest common divisors exist in UFDs, but not every integral domain in which greatest common divisors exist (known as a GCD domain) is a UFD. Every principal ideal domain is a UFD.

A Euclidean domain is an integral domain on which is defined a Euclidean division similar to that of integers. Every Euclidean domain is a principal ideal domain, and thus a UFD.

In a Euclidean domain, Euclidean division allows defining a Euclidean algorithm for computing greatest common divisors. However this does not imply the existence of a factorization algorithm. There is an explicit example of a field F such that there cannot exist any factorization algorithm in the Euclidean domain F[x] of the univariate polynomials over F.

In algebraic number theory, the study of Diophantine equations led mathematicians, during 19th century, to introduce generalizations of the integers called algebraic integers. The first ring of algebraic integers that have been considered were Gaussian integers and Eisenstein integers, which share with usual integers the property of being principal ideal domains, and have thus the unique factorization property.

Unfortunately, it soon appeared that most rings of algebraic integers are not principal and do not have unique factorization. The simplest example is Z [ − 5 ] , {\displaystyle \mathbb {Z} [{\sqrt {-5}}],} in which

9 = 3 ⋅ 3 = ( 2 + − 5 ) ( 2 − − 5 ) , {\displaystyle 9=3\cdot 3=(2+{\sqrt {-5}})(2-{\sqrt {-5}}),}

and all these factors are irreducible.

This lack of unique factorization is a major difficulty for solving Diophantine equations. For example, many wrong proofs of Fermat's Last Theorem (probably including Fermat's "truly marvelous proof of this, which this margin is too narrow to contain") were based on the implicit supposition of unique factorization.

This difficulty was resolved by Dedekind, who proved that the rings of algebraic integers have unique factorization of ideals: in these rings, every ideal is a product of prime ideals, and this factorization is unique up the order of the factors. The integral domains that have this unique factorization property are now called Dedekind domains. They have many nice properties that make them fundamental in algebraic number theory.

Matrix rings are non-commutative and have no unique factorization: there are, in general, many ways of writing a matrix as a product of matrices. Thus, the factorization problem consists of finding factors of specified types. For example, the LU decomposition gives a matrix as the product of a lower triangular matrix by an upper triangular matrix. As this is not always possible, one generally considers the "LUP decomposition" having a permutation matrix as its third factor.

See Matrix decomposition for the most common types of matrix factorizations.

A logical matrix represents a binary relation, and matrix multiplication corresponds to composition of relations. Decomposition of a relation through factorization serves to profile the nature of the relation, such as a difunctional relation.

• Euler's factorization method for integers
• Fermat's factorization method for integers
• Monoid factorisation
• Multiplicative partition
• Table of Gaussian integer factorizations

• Burnside, William Snow; Panton, Arthur William (1960) [1912], The Theory of Equations with an introduction to the theory of binary algebraic forms (Volume one), Dover
• Dickson, Leonard Eugene (1922), "First Course in the Theory of Equations", Nature, 109 (2746), New York: John Wiley & Sons: 773, Bibcode:1922Natur.109R.773., doi:10.1038/109773c0
• Fite, William Benjamin (1921), College Algebra (Revised), Boston: D. C. Heath & Co.
• Klein, Felix (1925), Elementary Mathematics from an Advanced Standpoint; Arithmetic, Algebra, Analysis, Dover
• Selby, Samuel M. (1970), CRC Standard Mathematical Tables (18th ed.), The Chemical Rubber Co.

• Wolfram Alpha can factorize too.