7.4  Existence and Uniqueness of the Determinant

In the first section of this chapter, signed area obeyed three rules: it is linear in each side, it vanishes when two sides are equal, and the unit square has area 1. In the second section we named functions obeying the first two rules, alternating multilinear forms, and in the third we built the sign of a permutation. This section puts the pieces together. We prove that on Fn there is exactly one alternating n-linear form that takes the value 1 on the standard basis, and we find a formula for it. That form is the determinant.

7.4.1. The goal and the plan

Here is the statement we are after.

There is exactly one function D:Fn×⋯×Fn→F (with n factors) that is n-linear (Definition 7.2.1), alternating (Definition 7.2.2), and satisfies D(e1,…,en)=1.

It has two halves of a very different nature. Uniqueness says: if such a D exists, its values are forced. Existence says: some function really has all three properties. We prove uniqueness first, because its proof is how we find the formula. Suppose D exists, expand everything using the rules, and see what is left. The formula that comes out is then the only candidate, and existence reduces to checking that the candidate obeys the rules.

Throughout, a list of n columns a1,…,an∈Fn and the matrix A=(a1⋯an)∈Mn(F) carry the same information, and we pass between them freely. We write aij for the i-th entry of aj, so that aj=a1je1+a2je2+⋯+anjen=∑i=1naijei.

7.4.2. Rearranging the arguments

Expanding a form on basis vectors will produce values like f(e3,e1,e2): the standard basis in a scrambled order. Alternating forms handle scrambling by a sign, and the sign is exactly the one built from inversions.

Lemma 7.4.1 (Permuting the arguments of an alternating form)

Let f be an alternating n-linear form on a vector space V over F, let v1,…,vn∈V, and let σ∈Sn. Then f(vσ(1),vσ(2),…,vσ(n))=sgn(σ)f(v1,v2,…,vn).

Idea. A transposition swaps two arguments, and each swap costs a factor −1 by Theorem 7.2.1. Every permutation is a product of transpositions, so the value picks up (−1)k for a product of k of them. Multiplicativity of the sign says that (−1)k is sgn(σ), whichever product we chose.

Proof. If n=1, then S1={id} and sgn(id)=1, since the identity has no inversions (Definition 7.3.2); there is nothing to prove. Let n≥2. For π∈Sn, write fπ=f(vπ(1),…,vπ(n)).

First let π∈Sn and let τ=(p q) be a transposition. The j-th argument of fπτ is vπ(τ(j)). For j≠p,q this is vπ(j); for j=p it is vπ(q); for j=q it is vπ(p). So the arguments of fπτ are those of fπ with positions p and q exchanged. By Theorem 7.2.1, swapping two arguments changes the sign, so (∗)fπτ=−fπ.

Now let σ∈Sn. By Theorem 0.10.4, σ=τ1τ2⋯τk for some transpositions τ1,…,τk with k≥1. Put π0=id and πm=τ1⋯τm, so that πm=πm−1τm. Applying (∗) k times, fσ=fπk=−fπk−1=⋯=(−1)kfπ0=(−1)kf(v1,…,vn). On the other hand, by Theorem 7.3.1 and Corollary 7.3.2, sgn(σ)=sgn(τ1)⋯sgn(τk)=(−1)k. This proves the lemma.

The proof chose a factorization of σ into transpositions, and different choices can have different lengths. The answer does not depend on the choice, because only the parity of k enters, and Corollary 7.3.3 says the parity is determined by σ. In the proof this is absorbed into the equation (−1)k=sgn(σ).

7.4.3. Uniqueness

We now expand an alternating form completely. It helps to see the case n=2 first.

Theorem 7.4.2 (Uniqueness of Alternating Forms)

Let f be an alternating n-linear form on Fn, and let a1,…,an∈Fn have entries aij as above. Then f(a1,…,an)=f(e1,…,en)∑σ∈Snsgn(σ)aσ(1)1aσ(2)2⋯aσ(n)n. In particular, an alternating n-linear form on Fn is determined by the single value f(e1,…,en).

Idea. For n=2, write a1=a11e1+a21e2 and a2=a12e1+a22e2. Linearity in the first argument, then in the second, gives four terms: f(a1,a2)=a11f(e1,a2)+a21f(e2,a2)=a11a12f(e1,e1)+a11a22f(e1,e2)+a21a12f(e2,e1)+a21a22f(e2,e2). The first and last terms vanish because two arguments are equal, and f(e2,e1)=−f(e1,e2). What is left is (a11a22−a21a12)f(e1,e2), the formula for n=2.

In general the four terms become one term for each way of choosing, for every argument j, which basis vector eφ(j) to take from aj. A choice is a function φ:{1,…,n}→{1,…,n}, and there are nn of them. The ones that pick some basis vector twice die. The survivors pick every basis vector once, so they are permutations, and Lemma 7.4.1 converts each into a sign.

Proof. Write [k]={1,…,k} for k≥1. We first expand the arguments one at a time.

Claim. For each k∈{0,1,…,n}, f(a1,…,an)=∑φ:[k]→[n]aφ(1)1⋯aφ(k)kf(eφ(1),…,eφ(k),ak+1,…,an), where the sum runs over all functions φ:[k]→[n]. For k=0 we read the right side as the single term f(a1,…,an).

Proof. We use induction on k. For k=0 both sides are f(a1,…,an). Suppose the claim holds for some k<n. In each term, the argument in position k+1 is ak+1=∑i=1nai,k+1ei. Since f is linear in argument k+1, f(eφ(1),…,eφ(k),ak+1,…,an)=∑i=1nai,k+1f(eφ(1),…,eφ(k),ei,ak+2,…,an). Substituting this into the induction hypothesis gives a sum over all pairs (φ,i) with φ:[k]→[n] and i∈[n]. Such a pair is the same thing as a function ψ:[k+1]→[n], namely ψ(j)=φ(j) for j≤k and ψ(k+1)=i; and each ψ arises from exactly one pair. The term for the pair (φ,i) is aψ(1)1⋯aψ(k+1),k+1f(eψ(1),…,eψ(k+1),ak+2,…,an). This is the claim for k+1.

Taking k=n in the claim, (∗∗)f(a1,…,an)=∑φ:[n]→[n]aφ(1)1⋯aφ(n)nf(eφ(1),…,eφ(n)). We sort the functions φ into two kinds.

Case 1: φ is not injective. Then φ(p)=φ(q) for some p≠q, so the arguments in positions p and q of f(eφ(1),…,eφ(n)) are equal. Since f is alternating, this value is 0, and the whole term is 0.

Case 2: φ is injective. Since [n] is finite, φ is bijective by Theorem 0.4.9, so φ=σ for some σ∈Sn. By Lemma 7.4.1 with vi=ei, f(eσ(1),…,eσ(n))=sgn(σ)f(e1,…,en).

Hence only the terms of Case 2 survive in (∗∗), there is exactly one for each σ∈Sn, and f(a1,…,an)=∑σ∈Snaσ(1)1⋯aσ(n)nsgn(σ)f(e1,…,en), which is the stated formula after moving the scalar f(e1,…,en) to the front. This proves the theorem.

So the three rules leave no freedom at all. If a normalized alternating form exists, it must be a1,…,an ⟼ ∑σ∈Snsgn(σ)aσ(1)1⋯aσ(n)n. Out of nn terms, the rules killed all but n!. For n=3 that is 6 out of 27.

7.4.4. Existence

The candidate has been found. Now we check it against the rules, forgetting where it came from.

Theorem 7.4.3 (Existence: the Leibniz Formula)

Define D:Fn×⋯×Fn→F by D(a1,…,an)=∑σ∈Snsgn(σ)aσ(1)1aσ(2)2⋯aσ(n)n, where aij is the i-th entry of aj. Then D is an alternating n-linear form on Fn, and D(e1,…,en)=1.

Idea. Linearity is visible: each product takes exactly one entry from each column. For alternation, suppose columns p and q are equal. Pair each σ with σ(p q). The two permutations have opposite signs. Their products use the same entries, except that the entries taken from columns p and q trade places, and those two columns are equal. So each pair contributes x−x=0.

Why pair terms instead of arguing “swapping the equal columns changes D to −D, so D=−D and D=0”? Because D=−D gives 2D=0, and over a field with 1+1=0, such as F2, that says nothing. The pairing argument never divides by 2.

Proof. Multilinearity. Fix k∈[n] and fix all columns except the k-th. For a fixed σ, the product aσ(1)1⋯aσ(n)n contains exactly one entry of ak, namely aσ(k)k, the σ(k)-th coordinate of ak. So the term for σ equals cσaσ(k)k, where cσ=sgn(σ)∏j≠kaσ(j)j does not depend on ak. The map ak↦aσ(k)k is linear, since coordinates of a sum are sums of coordinates and coordinates of cak are c times the coordinates. So each term is linear in ak, and hence so is their sum D.

Alternation. Suppose ap=aq with p<q, that is, aip=aiq for all i. Let τ=(p q), let An={σ∈Sn:sgn(σ)=1} be the even permutations (Definition 7.3.3), so that Sn∖An={σ∈Sn:sgn(σ)=−1} consists of the odd ones. Since sgn(στ)=sgn(σ)sgn(τ)=−sgn(σ) by Theorem 7.3.1 and Corollary 7.3.2, the map σ↦στ sends An to Sn∖An. It is a bijection An→Sn∖An, because ρ↦ρτ sends Sn∖An back to An and σττ=σ. Every σ∈Sn has sign 1 or −1, so D(a1,…,an)=∑σ∈An(∏j=1naσ(j)j−∏j=1naστ(j)j). Fix σ∈An and compare the two products factor by factor. For j≠p,q, στ(j)=σ(j), so the factors agree. For j=p, the second product has aσ(q)p=aσ(q)q, which is the factor of the first product at j=q. For j=q, the second product has aσ(p)q=aσ(p)p, the factor of the first product at j=p. So the two products consist of the same n factors in a different order, and they are equal because multiplication in F is commutative. Every bracket is 0, and D(a1,…,an)=0.

Normalization. For aj=ej, aij is 1 if i=j and 0 otherwise. The product aσ(1)1⋯aσ(n)n is therefore 1 if σ(j)=j for all j, and 0 otherwise. Only σ=id contributes, and sgn(id)=1. Hence D(e1,…,en)=1. This proves the theorem.

Uniqueness and existence together answer the question of the first section: the rules of signed area single out exactly one function.

7.4.5. The determinant

We now have a function of square matrices that is forced on us by three reasonable rules, and it deserves a name. It is the central notion of this chapter.

The determinant is the only way to attach a number to n vectors in Fn that is linear in each vector, vanishes when two vectors coincide, and gives 1 on the standard basis.

Definition 7.4.1 (Determinant)

Let n≥1 and A=(aij)∈Mn(F). The determinant of A is detA:=∑σ∈Snsgn(σ)aσ(1)1aσ(2)2⋯aσ(n)n ∈F. For columns a1,…,an∈Fn we also write det(a1,…,an) for the determinant of the matrix (a1⋯an).

In words: run over all n! permutations σ. For each, pick from column j the entry in row σ(j); this picks one entry in every column and, because σ is a bijection, one entry in every row. Multiply the n picked entries, attach the sign of σ, and add up. The determinant is defined only for square matrices, and its value lies in the field F of the entries.

Well-definedness. The formula is a finite sum, so it always gives a number. What the name “the determinant” really claims is that nothing else satisfies the three rules, and that is Theorem 7.4.2; Theorem 7.4.3 shows det does satisfy them. The sign sgn(σ) is fixed by counting inversions (Definition 7.3.2), and by Corollary 7.3.3 any factorization of σ into transpositions gives the same sign, so either method may be used to evaluate the formula.

Examples.

  • n=1. S1={id}, so det(a)=a. A degenerate case, but it is the base of every induction on size.
  • n=2. S2={id,(1 2)}, with signs 1 and −1. The identity picks a11a22, and (1 2) picks a21a12. So det(abcd)=ad−cb, the value derived from the rules in Example (ad−bc from the Rules), and the number whose non-vanishing decides invertibility in Theorem 0.9.8.
  • n=3. The six permutations of S3 and their signs (a 3-cycle has sign 1, a transposition −1) give detA=a11a22a33−a21a12a33−a31a22a13−a11a32a23+a21a32a13+a31a12a23, the terms coming from id,(1 2),(1 3),(2 3),(1 2 3),(1 3 2) in this order. For instance σ=(1 2 3) has σ(1)=2, σ(2)=3, σ(3)=1, and picks a21a32a13.
  • A numerical 3×3. For A=(213041521), the six terms in the same order are 8, −0, −60, −4, 0 and 5, so detA=8−60−4+5=−51.
  • Identity and zero. detIn=1 by the normalization in Theorem 7.4.3, and det0=0, since every product contains a zero factor.

Non-example by minimal change. Delete the sign. The permanent perA=∑σaσ(1)1⋯aσ(n)n is still n-linear in the columns, by the same argument as in Theorem 7.4.3, and still gives 1 at In. But it is not alternating: per(1111)=1+1=2, which is non-zero over R although the two columns are equal. The sign is exactly what makes the pairing argument cancel.

Why this definition. Each rule has a job. Without normalization, the uniqueness theorem only pins the form down up to the scalar f(e1,…,en). Without alternation, the permanent and many other functions qualify. Multilinearity is what let us expand at all. The choice of columns rather than rows is a convention: we show below that the transpose has the same determinant, so rows would give the same function. We have avoided “Sarrus’s rule” on purpose. It is a picture of the six-term formula for n=3 and nothing more, and for n=4 its diagonal pattern would produce 8 products while the Leibniz formula has 24.

Warning. The determinant is not linear. It is linear in each column separately, not in the matrix. For n=2, detI2=1 but det(I2+I2)=det(2I2)=4≠1+1. In general each product in the Leibniz formula has n factors, and scaling A scales every factor, so det(cA)=cndetA,in particulardet(−A)=(−1)ndetA. So det(−A)=detA for 2×2 matrices, not −detA.

Quick check. Which permutations give a non-zero term in the determinant of (001020300), and what is the determinant?

Solution. Column 1 has its only non-zero entry in row 3, column 2 in row 2, and column 3 in row 1. A non-zero term needs aσ(j)j≠0 for every j, so σ(1)=3, σ(2)=2, σ(3)=1, that is, σ=(1 3), of sign −1. The determinant is −3⋅2⋅1=−6.

The payoff for the forms of the second section is immediate. There we saw that the alternating n-linear forms on Fn form a vector space (Theorem 7.2.3) and promised its dimension.

Corollary 7.4.4 (Alternating n-Forms on Fn Form a Line)

The vector space of alternating n-linear forms on Fn has dimension 1, with basis (det). Explicitly, every alternating n-linear form f on Fn satisfies f(a1,…,an)=f(e1,…,en)det(a1,…,an).

Proof. By Theorem 7.4.3, det is an element of this space, and it is not the zero form since det(e1,…,en)=1. So the list (det) is linearly independent. By Theorem 7.4.2 and Definition 7.4.1, every f in the space satisfies the displayed formula, so f=f(e1,…,en)det lies in the span of det. Hence (det) is a basis, and the dimension is 1.

This corollary is the most useful form of uniqueness. To prove that some function g of n columns is a multiple of the determinant, check that g is multilinear and alternating, and compute g(e1,…,en). The multiplicativity of the determinant in the next section is proved exactly this way.

7.4.6. Transposes and triangular matrices

The definition read the matrix column by column. Reading it row by row gives nothing new.

Theorem 7.4.5 (Determinant of the Transpose)

For every A∈Mn(F), detA⊤=detA. Consequently, detA is also an alternating n-linear function of the rows of A, with value 1 at In.

Idea. A term of detA⊤ picks the entries ajσ(j), one in each row. The same n positions are also “one in each column”: position (j,σ(j)) is (σ−1(i),i) with i=σ(j). So the term is the term of detA for σ−1, and σ−1 has the same sign as σ.

Proof. The (i,j)-entry of A⊤ is aji (Definition 0.9.10), so by Definition 7.4.1 detA⊤=∑σ∈Snsgn(σ)a1σ(1)a2σ(2)⋯anσ(n). Fix σ. As j runs through [n], i=σ(j) runs through [n] exactly once, and j=σ−1(i). Since multiplication in F is commutative, we may reorder the factors by i: ∏j=1najσ(j)=∏i=1naσ−1(i)i. Next, sgn(σ)sgn(σ−1)=sgn(σσ−1)=sgn(id)=1 by Theorem 7.3.1, and both signs lie in {1,−1}, so sgn(σ−1)=sgn(σ). Therefore detA⊤=∑σ∈Snsgn(σ−1)∏i=1naσ−1(i)i=∑ρ∈Snsgn(ρ)∏i=1naρ(i)i=detA, where the second equality substitutes ρ=σ−1; this is allowed because σ↦σ−1 is a bijection of Sn onto itself, being its own inverse.

For the second statement, let r1,…,rn be the rows of A. The columns of A⊤ are r1⊤,…,rn⊤, so detA=detA⊤=det(r1⊤,…,rn⊤). The map r↦r⊤ is linear (Theorem 0.9.5), so by Theorem 7.4.3 the function (r1,…,rn)↦det(r1⊤,…,rn⊤) is linear in each ri. If rp=rq for some p≠q, then rp⊤=rq⊤, and the value is 0. Finally In⊤=In. This proves the theorem.

So every statement about columns has a twin about rows. In particular Theorem 7.2.1 applies to rows: swapping two rows changes the sign of det, and adding a multiple of one row to another does not change it. The next section turns this into a method of computation.

Before that, one class of matrices whose determinant can be read off directly. Recall that A is upper triangular if aij=0 whenever i>j (Definition 0.9.5).

Theorem 7.4.6 (Determinant of a Triangular Matrix)

If A∈Mn(F) is upper triangular or lower triangular, then detA=a11a22⋯ann.

Idea. Which permutations can pick only entries on or above the diagonal? Picking aσ(j)j with σ(j)≤j for every j forces σ=id: the numbers σ(1),…,σ(n) are 1,…,n in some order, so they have the same sum, and none of them is allowed to exceed its partner.

Proof. Suppose A is upper triangular, and let σ∈Sn with σ≠id. We claim σ(j)>j for some j. Otherwise σ(j)≤j for all j. Since σ is a bijection of [n], ∑j=1nσ(j)=∑j=1nj, so ∑j=1n(j−σ(j))=0 is a sum of non-negative integers equal to 0. Hence σ(j)=j for every j, contradicting σ≠id. For such j, aσ(j)j=0 by upper triangularity, so the term of σ is 0. Only the identity remains, and detA=sgn(id)a11⋯ann=a11⋯ann.

If A is lower triangular, then A⊤ is upper triangular with the same diagonal, and detA=detA⊤=a11⋯ann by Theorem 7.4.5. This proves the theorem.

For example, det(27−103500−4)=2⋅3⋅(−4)=−24: six terms, five of which vanish. Diagonal matrices are a special case, detdiag(d1,…,dn)=d1⋯dn. Together with Chapter 3’s result that an upper triangular matrix is invertible exactly when its diagonal entries are non-zero (Lemma 3.6.2), this already shows that for upper triangular matrices, A is invertible if and only if detA≠0. The next section proves this for all square matrices.

7.4.7. Determinants over a commutative ring

The characteristic polynomial at the end of this chapter is the determinant of a matrix whose entries are polynomials, such as (x−1−2−3x−4). Polynomials form no field, since x has no inverse in F[x]. So we record exactly what this section used.

Remark. Commutative rings. A commutative ring is a set R with an addition and a multiplication satisfying all the axioms of a field (Definition 0.7.1) except, possibly, the existence of multiplicative inverses (F8). The integers Z and the polynomials F[x] are examples (Theorem 0.8.1). For A∈Mn(R), the Leibniz formula of Definition 7.4.1 still defines detA∈R, reading sgn(σ)r as r or −r.

Every proof in this section uses only addition, subtraction and multiplication, together with commutativity, and never a multiplicative inverse. So the following hold verbatim over any commutative ring R, with “linear in each argument” meaning additive and compatible with multiplication by elements of R: Lemma 7.4.1, Theorem 7.4.2, Theorem 7.4.3, Theorem 7.4.5 and Theorem 7.4.6. The only outside input is the swap rule of Theorem 7.2.1, and it too is division-free: expanding 0=f(…,u+v,…,u+v,…) by additivity leaves f(…,u,…,v,…)+f(…,v,…,u,…)=0. The parts of Theorem 7.2.1 about dependent arguments are not used; they involve solving for one vector in terms of others, which needs division.

In the characteristic-polynomial section we will apply these results with R=F[x]. The next two sections flag, in the same way, which of their proofs survive the passage to a ring.

7.4.8. Exercises

A. Check your understanding

Exercise (A1)
  • (a) State the three properties that characterize det:Mn(F)→F, and write down the Leibniz formula.
  • (b) How many terms does the Leibniz formula have for n=4? How many of the corresponding permutations have sign 1?
  • (c) True or false: det(A+B)=detA+detB for all A,B∈M2(R). Justify your answer.
  • (d) True or false: det(−A)=−detA for all A∈M3(R). Justify your answer.
  • (e) In the proof of Theorem 7.4.2, why does a term with a non-injective φ vanish?
  • (f) Explain why the existence proof pairs σ with σ(p q) instead of arguing “D=−D, so D=0”.

Solution.

  • (a) det is n-linear in the columns, alternating (zero whenever two columns are equal), and detIn=1. By Theorem 7.4.2 and Theorem 7.4.3 it is the only such function, and detA=∑σ∈Snsgn(σ)aσ(1)1⋯aσ(n)n.
  • (b) One term per permutation: 4!=24. In the proof of Theorem 7.4.3, σ↦σ(1 2) is a bijection from the permutations of sign 1 to those of sign −1, so each set has 12 elements.
  • (c) False. For A=B=I2, det(A+B)=det(2I2)=4 but detA+detB=2.
  • (d) True. Each Leibniz product has 3 factors, each multiplied by −1, so det(−A)=(−1)3detA=−detA. (For n=2 it is false.)
  • (e) If φ(p)=φ(q) with p≠q, the form is evaluated with the equal arguments eφ(p) in positions p and q, and an alternating form vanishes there.
  • (f) D=−D only gives 2D=0. Over a field in which 2=0, such as F2, this does not imply D=0. The pairing shows that the terms cancel in pairs, which works over every field.

B. Practice

Exercise (B1: Six-term determinants)

Using the Leibniz formula, listing all six terms, compute the determinants of A=(312104251),B=(1−1203142−2)∈M3(R). Then compute detA⊤ in the same way, and hence verify Theorem 7.4.5 for this A.

Solution. We use the order id,(1 2),(1 3),(2 3),(1 2 3),(1 3 2) and the six-term formula from the examples after Definition 7.4.1: a11a22a33, −a21a12a33, −a31a22a13, −a11a32a23, a21a32a13, a31a12a23. For A: 3⋅0⋅1=0; −1⋅1⋅1=−1; −2⋅0⋅2=0; −3⋅5⋅4=−60; 1⋅5⋅2=10; 2⋅1⋅4=8. Hence detA=0−1+0−60+10+8=−43.

For B: 1⋅3⋅(−2)=−6; −0⋅(−1)⋅(−2)=0; −4⋅3⋅2=−24; −1⋅2⋅1=−2; 0⋅2⋅2=0; 4⋅(−1)⋅1=−4. Hence detB=−6+0−24−2+0−4=−36.

For A⊤=(312105241): 3⋅0⋅1=0; −1⋅1⋅1=−1; −2⋅0⋅2=0; −3⋅4⋅5=−60; 1⋅4⋅2=8; 2⋅1⋅5=10. The sum is −43=detA, as Theorem 7.4.5 predicts; the six terms are the same numbers, with the two 3-cycles exchanged, because (1 2 3)−1=(1 3 2).

Exercise (B2: The determinant of Pσ)

Let σ∈Sn and let Pσ=(eσ(1)⋯eσ(n)) be its permutation matrix (Definition 3.4.2).

  • (a) Prove that detPσ=sgn(σ).
  • (b) Hence compute det(001100010) and det(0100100000010010).

Solution.

  • (a) By Theorem 7.4.3, det is an alternating n-linear form on Fn. By Lemma 7.4.1 with vi=ei, detPσ=det(eσ(1),…,eσ(n))=sgn(σ)det(e1,…,en)=sgn(σ)detIn=sgn(σ).
  • (b) The first matrix has columns e2,e3,e1, so it is Pσ with σ(1)=2, σ(2)=3, σ(3)=1, that is, σ=(1 2 3). A 3-cycle is (1 3)(1 2) (Example (Computing with cycles)), a product of two transpositions, so its sign is 1 and the determinant is 1. The second matrix has columns e2,e1,e4,e3, so it is Pσ with σ=(1 2)(3 4), of sign (−1)2=1 by Theorem 7.3.1 and Corollary 7.3.2. Its determinant is 1.
Exercise (B3: A zero row)

Let A∈Mn(F) have a zero row. Prove directly from Definition 7.4.1, without using Theorem 7.4.5, that detA=0. Then prove the same for a zero column.

Solution. Suppose row i of A is zero, and let σ∈Sn. Since σ is surjective, there is j with σ(j)=i. The factor aσ(j)j=aij of the term for σ is 0, so the term is 0. This holds for every σ, so detA=0.

Suppose column j is zero. For every σ, the term for σ contains the factor aσ(j)j=0, so again every term vanishes and detA=0.

C. Going deeper

Exercise (C1: Too many zeros)

Let A∈Mn(F).

  • (a) Prove that if more than n2−n entries of A are zero, then detA=0.
  • (b) Show that the bound is sharp: give, for each n, a matrix with exactly n2−n zero entries and non-zero determinant.

Solution.

  • (a) If more than n2−n entries are zero, then fewer than n entries are non-zero. For σ∈Sn, the term for σ is a product of the n entries in positions (σ(1),1),…,(σ(n),n). These positions are distinct, since their column indices are distinct. They cannot all hold non-zero entries, because there are fewer than n such entries. So at least one factor is 0, and the term is 0. Hence detA=0.
  • (b) In has exactly n2−n zero entries, and detIn=1≠0.
Exercise (C2: Counting the cost)

A computer performs about 109 multiplications per second.

  • (a) Show that evaluating the Leibniz formula term by term for an n×n matrix takes n!(n−1) multiplications. Estimate the time this takes for n=20.
  • (b) Gaussian elimination (Chapter 3) brings A to upper triangular form. Suppose no row swaps are needed. At step k (for k=1,…,n−1), each of the n−k rows below the pivot requires one division to find its multiplier and n−k multiplications to update its entries in columns k+1,…,n. Show that the total number of multiplications and divisions is (n−1)n(n+1)3, and evaluate it for n=20.

In the next section we will see that the determinant can be read off from this triangular form with n−1 further multiplications.

Solution.

  • (a) There is one term for each σ∈Sn, that is, n! terms, and each is a product of n entries, which takes n−1 multiplications. For n=20, 20!⋅19=46225138155356160000≈4.6×1019. At 109 per second this takes about 4.6×1010 seconds, roughly 1460 years.
  • (b) Step k costs (n−k)(1+(n−k)) operations. With m=n−k running from 1 to n−1, the total is ∑m=1n−1m(m+1)=∑m=1n−1m2+∑m=1n−1m=(n−1)n(2n−1)6+(n−1)n2=(n−1)n((2n−1)+3)6=(n−1)n(n+1)3. For n=20 this is 19⋅20⋅213=2660, a few microseconds. The Leibniz formula is a definition and a proof tool, not an algorithm.
Exercise (C3: The permanent, and a change of field)

For A∈Mn(F) define perA=∑σ∈Snaσ(1)1⋯aσ(n)n.

  • (a) Prove that per is n-linear in the columns and perIn=1.
  • (b) Let n≥2. Prove that over R, per is not alternating.
  • (c) Prove that over F2, perA=detA for every A∈Mn(F2). Explain why this does not contradict (b) and Theorem 7.4.2.

Solution.

  • (a) The proof of multilinearity in Theorem 7.4.3 never used the signs: for fixed k, the term for σ is a constant times aσ(k)k, which is linear in the k-th column. The normalization argument there also applies: only σ=id contributes at In, with product 1.
  • (b) Let J∈Mn(R) have every entry equal to 1. Its first two columns are equal, but every term of perJ is 1, so perJ=n!≠0 in R.
  • (c) In F2, −1=1, so sgn(σ)⋅r=r for every r, and the two formulas agree term by term. There is no contradiction: (b) is a statement over R, where n!≠0. Over F2, perJ=n!⋅1=0 for n≥2, and per=det is alternating, as Theorem 7.4.3 guarantees. Uniqueness is about a fixed field, and over each field there is exactly one normalized alternating form.

Found a typo or a mistake? Report it, or read this page’s source

The text of this book is licensed under CC BY-NC-SA 4.0.