On Some Euclidean Properties of Matrix Algebras Pierre Lezowski
On some Euclidean properties of matrix algebras Pierre Lezowski
To cite this version:
Pierre Lezowski. On some Euclidean properties of matrix algebras. Journal of Algebra, Elsevier, 2017, 486, pp.157–203. 10.1016/j.jalgebra.2017.05.018. hal-01135202v3
HAL Id: hal-01135202 https://hal.archives-ouvertes.fr/hal-01135202v3 Submitted on 30 Jun 2017
HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. ON SOME EUCLIDEAN PROPERTIES OF MATRIX ALGEBRAS
PIERRE LEZOWSKI
Abstract. Let R be a commutative ring and n ∈ Z>1. We study some Euclidean properties of the algebra Mn (R) of n by n matrices with coefficients in R. In particular, we prove that Mn (R) is a left and right Euclidean ring if and only if R is a principal ideal ring. We also study the Euclidean order type of Mn (R). If R is a K-Hermite ring, then Mn (R) is (4n − 3)-stage left and right Euclidean. We obtain shorter division chains when R is an elementary divisor ring, and even shorter ones when R is a principal ideal ring. If we assume that R is an integral domain, R is a Bézout ring if and only if Mn (R) is ω-stage left and right Euclidean.
1. Introduction In this paper, all rings are nonzero, with unity, but not necessarily com- mutative. An integral domain is a commutative ring with no nontrivial zero divisor. A principal ideal ring (or PIR for short) is a commutative ring in which every ideal is principal. A principal ideal domain (or PID for short) is an integral domain which is a PIR. Given a ring A, we denote by A• the set A 0 and by A the units of A. \ { } × Given a ring A and integers n,m > 0, Mm,n (A) is the set of matrices of elements of R with m rows and n columns; Mn (A) = Mn,n (A); GLn (A) is the subset of Mn (A) of units of Mn (A). Whenever R is commutative, Mn (R) is an algebra, and we are especially interested in its Euclidean properties. In the classical sense, we say that a ring A is right Euclidean if there exists some function ϕ : A Z>0 such that for all a, b A, b = 0, there exists q A such that −→ ∈ 6 ∈ a = bq or ϕ(a bq) < ϕ(b). − However, with this definition, A = Mn (Z) cannot be right Euclidean when n Z>1(see [Kal85, Theorem 2]), so instead, we will use a broader defini- tion,∈ following Samuel [Sam71]. Let us denote by the class of all ordinal numbers. O Definition 1.1. Let A be a ring. We say that A is right Euclidean if there exists a function ϕ : A• such that for all a, b A, b = 0, there exists q A such that −→ O ∈ 6 ∈ (1) a = bq or ϕ(a bq) < ϕ(b). − Date: September 19, 2016. 2010 Mathematics Subject Classification. Primary: 13F07; Secondary: 11A05. Key words and phrases. Euclidean rings, 2-stage Euclidean rings, Euclidean algorithm, Principal Ideal Rings, division chains, elementary divisor ring, K-Hermite ring. 1 2 PIERRE LEZOWSKI
Such a ϕ is then called a right Euclidean stathm (or a right Euclidean func- tion). Obviously, we may define similarly left Euclidean rings and left Euclidean stathms by replacing bq with qb in (1). With this definition, Brungs proved the following property. Proposition 1.2 ([Bru73, Theorem 1]). If A is a (not necessarily commu- tative) left Euclidean ring without nontrivial zero divisors, then Mn (A) is a left Euclidean ring for any n Z 1. ∈ ≥ We will establish the following result.
Theorem 4.1. Let R be a commutative ring and n Z>1. Then Mn (R) is right and left Euclidean if and only if R is a principal∈ ideal ring. To prove it, we will use some technical tools and notations, introduced in Section 2, and proceed in two steps. We will rely on the fact that a PIR is a K-Hermite ring, that is to say that every matrix admits triangular reduction, and even an elementary divisor ring, that is to say that every matrix admits diagonal reduction1. First, we will prove Theorem 4.1 over PIDs in Section 3, which will allow us to extend it to PIRs in Section 4. We will see in Section 5 under which conditions we can compute a quotient of the right Euclidean division (1) for the stathm that we build. As an application, we will see that we can compute continued fractions in a matrix algebra over a PID. In Definition 1.1, the range of the Euclidean stathm may be arbitrary, but for a given right Euclidean ring A, we can try to find a right Euclidean stathm whose range is as “small” as possible. This is formalized by the notion of Euclidean order type of A. Section 6 will be devoted to the study of the Euclidean order type of Mn (R) when R is a PIR. Finally, we will study another generalization of the Euclidean property in Section 7. Instead of allowing ordinals in the range of the stathm, we still consider ϕ : A• Z>0, but we allow several divisions on the right: starting from the pair (a,−→ b), we continue with a pair (b, a bq) for some q A, and 2 − ∈ so forth . After k divisions, we want the remainder rk to satisfy
(2) rk = 0 or ϕ(rk) < ϕ(b). If for all pair of elements of A, we can obtain a k-stage division chain with the k-th remainder rk satisfying (2), we say that A is k-stage right Euclidean. If rk = 0, we say that the division chain is terminating. If for all a, b A, b = 0, there exists a terminating division chain starting from (a, b), we∈ say that6 A is ω-stage right Euclidean. A right Euclidean ring is necessarily ω- stage right Euclidean, but the converse is false in general since such a ring may have non-principal ideals. Alahmadi, Jain, Lam, and Leroy proved the following result about the ω-stage right Euclidean properties of matrix rings. Proposition 1.3 ([AJLL14, Theorem 14]). If A is a (not necessarily com- mutative) ω-stage right Euclidean ring, then so is Mn (A) for any n Z 1. ∈ ≥ 1 See Section 2 for precise definitions. 2 See Section 7 for precise definitions. ON SOME EUCLIDEAN PROPERTIES OF MATRIX ALGEBRAS 3
An immediate consequence of Theorem 4.1 is that Mn (R) is ω-stage right Euclidean if R is a PIR and n Z>1. But we will show the following more precise result in Section 7. ∈
Theorem 7.3. Let R be a commutative ring and n Z>1. Then we have the following properties. ∈ (1) If R is a K-Hermite ring, then for every pair (A, B) M (R) ∈ n × M (R)•, there exists a (4n 3)-stage terminating division chain in n − Mn (R) starting from (A, B). In particular, Mn (R) is ω-stage left and right Euclidean. (2) If R is an elementary divisor ring (e.g. if R is a PIR), then for every pair (A, B) M (R) M (R)•, there exists a (2n 1)-stage ∈ n × n − terminating division chain in Mn (R) starting from (A, B). (3) If R is a PIR, then Mn (R) is 2-stage right and left Euclidean.
Therefore, when R is an integral domain, we can characterize when Mn (R) is ω-stage right and left Euclidean.
Corollary 7.4. Let R be an integral domain and n Z>1. Then Mn (R) is ω-stage left and right Euclidean if and only if R is a∈ Bézout ring, that is to say for all a, b R, there exists d R such that aR + bR = dR. ∈ ∈ 2. Generalities and first remarks 2.1. Notation and terminology. Consider a ring A, a commutative ring R, and m, n, p Z 1. In R, we say that a divides b, denoted by a b, if bR aR. Most of∈ the≥ definitions and the results in this paragraph are| due to Kaplansky⊆ [Kap49]3. We denote by diag(b1, . . . , bn)m,p the matrix in Mm,p (A) with diagonal coefficients b1, . . . , bn. For short, diag(b1, . . . , bn) = diag(b1, . . . , bn)n,n. We write 1n for the identity matrix of size n, 0m,n for the zero matrix with m rows and n columns, 0n = 0n,n. A ring A is a right Bézout ring if for all a, b A, aA + bA = dA, for some d A. Such a d is called a greatest common∈ left divisor of a and b, or gcld for∈ short. We define similarly left Bézout rings and gcrds (i.e. greatest common right divisors). The stable rank of a ring A is the infimum of the positive integers n such that for all a , . . . , a A, 0 n ∈ a0A + + anA = A = b1, . . . , bn A, (3) · · · ⇒∃ ∈ (a + b a )A + + (a + b a )A = A. 1 1 0 · · · n n 0 We denote it by sr A. If sr A = l, then (3) holds for any n l [Vas71, Theorem 1]. The stable rank can be defined as in (3) using left ideals≥ instead of right ideals; the value sr A coincides [Vas71, Theorem 2]. The stable rank of A and matrix rings over A are connected [Vas71, Theorem 3]: sr A 1 sr M (A)=1+ − . n n In the formula above, x is the least integer exceeding x. In particular, ⌈ ⌉ Mn (A) has stable rank 1 if and only if A has stable rank 1. 3 We call “K-Hermite” what he calls “Hermite” to follow the current trend. 4 PIERRE LEZOWSKI
We say that A is a right K-Hermite ring if for all a, b A, there exists ∈ Q GL2 (A) and d A such that a b Q = d 0 . Any right K-Hermite ring∈ is a right Bézout∈ ring. Besides, if A is a right K-Hermite ring, then for all M Mn (A), there exists T GLn (A) such that AT is lower triangular. A left∈ K-Hermite ring satisfies∈ the following condition: for all a, b A, a d ∈ there exists Q GL (A) and d A such that Q = . A K-Hermite ∈ 2 ∈ b 0 ring is a right and left K-Hermite ring. If A is a K-Hermite ring, then Mn (A) is a right K-Hermite ring ([Kap49, Theorem 3.6]). If A is a K-Hermite ring with no nontrivial zero divisors, then for A, B M (A), the gcld of A and ∈ n B is unique up to multiplication by an element of GLn (A) on the right (see [Kap49, Theorem 3.8]). In commutative K-Hermite rings, we can simplify by a (good choice of a) gcd, and it is actually a characterization of these rings. Lemma 2.1 ([GH56, Theorem 3]). Let R be a commutative ring. Then R is a K-Hermite ring if and only if for all a, b R, there exist a , b , d R ∈ ′ ′ ∈ such that a = da′, b = db′, and
a′R + b′R = R. We can clearly extend it to three elements, and we will especially use it in this form: Let R be a commutative K-Hermite ring, a, b, c R. Then there exist a , b , c , d R such that a = da , b = db , c = dc and∈ ′ ′ ′ ∈ ′ ′ ′ a′R + b′R + c′R = R. Given A, B M (R), we say that A and B are equivalent, which is ∈ m,p denoted by A B if there exist X GLp (R), Y GLm (R) such that B = ∼ ∈ ∈ 4 Y AX. We say that R is an elementary divisor ring if for any m, p Z 1, any A M (R) is equivalent to a diagonal matrix diag(b , b , . . . ,∈ b ) ≥ , ∈ m,p 1 2 n m,p where b1 b2 . . . bn. In this case, b1 divides every coefficient of the matrix A. Any elementary| | | divisor ring is a K-Hermite ring. The following property will be a crucial tool. Lemma 2.2 ([MM82, Proposition 8]). Let A be a right K-Hermite ring. Then the stable rank of A satisfies sr A 2. ≤ If R is a PID, then it is an elementary divisor ring, so a K-Hermite ring. You can refer to [Jac85, Section 3.7] for details. The reductions of matrices with coefficients in R into triangular or diagonal ones can be computed, provided that for any a, b R, we know how to compute λ, µ, d R such that ∈ ∈ aR + bR = dR, aλ + bµ = d. n Besides, recall that for any M Mn (R), there exist (bi)1 i n R such ∈ ≤ ≤ ∈ that b1 b2 . . . bn and | | | M diag(b , . . . , b ). ∼ 1 n In this case, the elements (bi)1 i n are unique up to multiplication by a unit of R and are called the invariant≤ ≤ factors of M. Such a reduction is called the 4 Kaplansky’s definition is not limited to a commutative context, but we will not need such a generality. ON SOME EUCLIDEAN PROPERTIES OF MATRIX ALGEBRAS 5
Smith normal form of M. The largest integer r such that br = 0 is the rank rk(M) of M. It corresponds to the classical notion of rank in6 vector spaces, because R can be embedded into its field of fractions. In particular, a matrix M Mn (R) has rank n if and only if det M = 0, and more generally, the rank∈ is equal to the maximal order of a nonzero6 minor. Smith normal form has some further properties.
Lemma 2.3. Let R be a PID and n Z 1. ∈ ≥ (a) Let M M (R)• and let b , b , . . . , b R be the invariant factors ∈ n 1 2 r ∈ of M. For any 1 l r, l b is a gcd of the l l minors of M. ≤ ≤ i=1 i × In particular, b1 is a gcd of the coefficients of M. a b Q (b) Let M = c d M2 (R). If the greatest common divisor of a, b, c, d is 1, then ∈ M diag(1, ad bc). ∼ − Proof. Item (a) is a reformulation of [Jac85, Theorem 3.9]. For (b), write M diag(b , b ), for elements b b in R. Thanks to (a), we can take ∼ 1 2 1| 2 b1 = 1. Besides, det M and b1b2 coincide up to multiplication by a unit, which completes the proof. 2.2. Basic remarks. In Definition 1.1, we have distinguished right and left Euclidean rings, but such a care will be useless in our context.
Proposition 2.4. Let R be a commutative ring and n Z>1. Then Mn (R) is right Euclidean if and only if it is left Euclidean. ∈
Proof. Let f be any function M (R)• . We define n −→ O T Mn (R)• f : −→ O T . M f(M ) 7−→ Then for any A,B,Q M (R), ∈ n T T T T T T T A = BQ A = Q B and f (A BQ)= f(A Q B ). ⇐⇒ − − Hence, f is a right Euclidean stathm if and only f is a left Euclidean stathm.
Therefore, to prove Theorem 4.1, we will only need to deal with right Euclidean stathms.
Proposition 2.5. Let R be a commutative ring and n Z 1. If every right ∈ ≥ ideal of Mn (R) is principal, then R is a principal ideal ring. Proof. This is certainly very classical, but we include the proof to emphasize its simplicity. Let I be an ideal of R. We define
I = (ai,j)1 i,j n Mn (R) , for any 1 j n, a1,j I . { ≤ ≤ ∈ ≤ ≤ ∈ } It is clear that I is a right ideal of Mn (R). Let us consider det I = det M, M I . Then, the fact that I is an ideal and Leibniz formula {for determinants∈ } imply that det I I. Besides, for any a I, define ⊆ ∈ A = diag(a, 1,..., 1) M (R) . ∈ n Then A I and a = det A det I. Therefore, we also have I det I, which proves that∈ det I = I. ∈ ⊆ 6 PIERRE LEZOWSKI
But there exists α Mn (R) such that I = αMn (R). Then I = det I = (det α) R, which completes∈ the proof. Remark 2.6. If we do not assume R to be commutative, this property is false in general: consider the Weyl algebra A1 in characteristic 0. Then A1 admits non-principal right ideals, but for any n 2, every right ideal of Mn (A1) is principal, see [MR01, 7.11.7 and 7.11.8, p.≥ 293].
2.3. Length in a PID. Let R be a PID and x R•. Since R is a unique factorization domain, x may be decomposed into∈ a finite product of prime elements n ei x = u pi , Yi=1 where n Z 0, u R×, for any 1 i n, pi R is prime and ei Z>0. The decomposition∈ ≥ ∈ is unique up to≤ multiplication≤ ∈ by units and order.∈ We n set ℓ(x)= i=1 ei, which defines a function
R• Z 0 P ℓ : ≥ . x −→ ℓ(x) 7−→ Remark that ℓ is invariant under multiplication by a unit. Besides, if a, b are elements of R such that a divides b and b = 0, then ℓ(a) ℓ(b) and the equality holds if and only if a and b are associate,6 that is to say≤ there exists u R such that b = au. ∈ × 2.4. Explicit stable rank 2 in a PID. The following classical lemma will be very useful. It can be seen as an easy consequence of Lemma 2.2, but we give its proof to see how explicit it is; this will prove useful for Section 5. Lemma 2.7. Let R be a PID. Then sr R 2. More precisely, for any a, b, c R which are not all equal to 0, there≤ exist z,t R, such that gcd(a +∈cz, b + ct)= gcd(a, b, c), which is nonzero and divides∈ c. Proof. We will proceed in two steps. (a) First, consider a, b, c R such that b = 0 and gcd(a, b, c) = 1, we will prove that there exists z ∈ R such that gcd(6 a + cz, b) = 1. If a and b are coprime, we can take z = 0∈. If not, write a decomposition of b = 0 as follows: 6 l αi (4) b = di z, Yi=1 where l Z>0, (di)1 i l is a family of distinct and non-associated primes in ∈ ≤ ≤ R, for any 1 i l, αi Z>0, di divides a, but di does not divide z. Take a prime p R≤ such≤ that∈p divides b. If p divides a, then p is associated to ∈ some di, for 1 i l, and p does not divide z. Therefore, if p divides a+cz, then it necessarily≤ ≤ divides a + cz a = cz, so it divides c. Then p divides a, b, and c, which are coprime. This− is impossible, so p does not divide a + cz. If p does not divide a, then p divides z. Thus, p does not divide a + cz in this case either. Consequently, gcd(a + cz, b) = 1. (b) Now, we consider a, b, c R which are not all equal to 0. If c = 0, take z = t = 0. From now on,∈ assume that c = 0. Take t 0, 1 such 6 ∈ { } that b + ct = 0. Set d = gcd(a, b, c), consider a = a , b = b+tc , c = c 6 ′ d ′ d ′ d ON SOME EUCLIDEAN PROPERTIES OF MATRIX ALGEBRAS 7 and apply (a): there exists z R such that gcd(a′ + c′z, b′) = 1. Then gcd(a + cz, b + ct)= d. ∈ Remark 2.8. In the proof above, we can compute z in (4) without actually computing a decomposition of b into a product of primes, it is enough to compute some gcds. Proof. For a, b R, b = 0, we want to find a pair (d, z) R2 such that b = dz, gcd(d, z)∈ = 1, and6 for any prime p dividing b, p divides∈ a if and only if p divides d. We build inductively a pair (dm, zm) of elements of R. Write d1 = gcd(b, a) and b = d z for some z R. If d and z are coprime, then m = 1 and 1 1 1 ∈ 1 1 we are done. If not, assume that we have (di, zi) such that b = dizi. If gcd(di, zi), we are done, set m = i. If not, set di+1 = di gcd(di, zi) and write z = di+1zi+1. As at each step di is a divisor of b and a strict divisor of di+1, we are done in a finite number of steps: we obtain
z = dmzm, where gcd(dm, zm) = 1. Besides, it is straightforward that gcd(b, a) = d1 divides d . Notice that for any i 1, d = d gcd(d , z ), so any prime m ≥ i+1 i i i divisor of di+1 is a prime divisor of di. Consequently, any prime divisor of dm is a prime divisor of d1 = gcd(b, a). Take a prime p dividing b. If p divides a, then it divides gcd(b, a), so p divides dm. Conversely, if p divides dm, then it divides d1 = gcd(b, a), so p divides a. Hence the pair (d, z) = (dm, zm) is convenient.
2.5. Conventions and notations for ordinals. We follow the notation used by Clark [Cla15], that is to say we denote by ω the least infinite ordinal, and for ordinal arithmetic, we fix the ordinal addition so that ω + 1 >ω = 1+ω, and for the multiplication 2ω = ω+ω>ω2= ω. For short, for r Z>0, r ∈ a , b , 1 i r, we write a ωbi = a ωb1 + a ωb2 + + a ωbr . i i ∈ O ≤ ≤ i=1 i 1 2 · · · r We denote by the Hessenberg sum of ordinals, that is to say for k Z>0, ⊕ P ∈ (ai)0 i k , (bi)0 i k finite sequences of nonnegative integers, ≤ ≤ ≤ ≤ k k k k i k i k i aiω − biω − = (ai + bi)ω − . ! ⊕ ! ! Xi=0 Xi=0 Xi=0 For n Z , α , we write n α = α α α. ∈ >0 ∈ O ⊗ ⊕ ⊕···⊕ n times Consider a right Euclidean ring A. In Definition 1.1, the right Euclidean stathm ϕ is not defined at 0. Following| Clark,{z we define} ϕ(0) to be the smallest α such that for all a A , ∈ O ∈ • (5) ϕ(a) < α. Now, we associate to A the following ordinal number, called (right) Euclidean order type: e(A) = inf ϕ(0), ϕ : A , ϕ right Euclidean stathm . { −→ O } 8 PIERRE LEZOWSKI
In other words, e(A)= θ(0), where θ is the smallest right Euclidean stathm for A, as defined by Samuel [Sam71] (or “bottom Euclidean function”, with Clark’s terminology): it is the function defined by A θ : • x −→inf O φ(x), φ : A , φ right Euclidean stathm ; 7−→ { • −→ O } it is a right Euclidean stathm.
Remark 2.9. Let R be a commutative ring, n Z 1 so that Mn (R) is right ∈ ≥ Euclidean. Let θ be the smallest right Euclidean stathm for Mn (R). Then for any m,m′ Mn (R) such that m m′, we have θ(m) = θ(m′). In particular, if R∈is an elementary divisor∼ ring, θT = θ and θ is the smallest left Euclidean stathm for Mn (R). For this reason, we will write Euclidean order type instead of right Euclidean order type in what follows. Proof. It follows immediately from the fact that the function
M (R)• θ˜ : n m −→ Oinf θ(m ), m m 7−→ { ′ ′ ∼ } is a right Euclidean stathm verifying θ˜ θ. Therefore, θ˜ = θ. ≤ Lemma 2.10. Let A be a right Euclidean ring and θ : A be the • −→ O smallest right Euclidean stathm. Take x A• and A xA such that 0 is a system of representatives of A∈/xA. ThenS ⊆ \ S ∪ { } θ(x) sup inf θ(y + xa) + 1. ≤ y a A ∈S ∈ Proof. This is a consequence of Motzkin’s construction. For α , define 0 ∈ O Aα = z A, θ(z) α and Aα = β<αAβ. Then we have (see [Sam71] or {5 ∈ ≤ } ∪ [Cla15] ) . A = b A, the composite map A0 0 ֒ A ։ A/xA is onto α ∈ α ∪ { } −→ − Fix α = sup y infa A θ(y +xa)+1, we will prove that x Aα. Let yˆ+x A ∈S ∈ ∈ ∈ A/xA xA , there exists y such that yˆ+ xA = y + xA. By definition of \ { } ∈ S 0 α, there exists a A such that θ(y + xa) < α, therefore y + xa Aα, which concludes the proof∈ as y + xa + xA =y ˆ + xA. ∈
3. A left and right Euclidean stathm for matrix algebras over a PID
3.1. Statement and first remarks. The purpose of this section will be to establish the following result.
Theorem 3.1. Let R be a PID and n Z>1. Then Mn (R) is a left and right Euclidean ring. ∈ Let R be a PID and n Z . We define ∈ >1 Mn (R)• −→ O rk M rk M i ρn : M ℓ(b )ω if 7−→ i=1 i − b1, b2, . . . , brk M are the invariant factors of M. P 5 They deal with the commutative case but never use the commutativity hypothesis in this context. ON SOME EUCLIDEAN PROPERTIES OF MATRIX ALGEBRAS 9
The function ρn is well-defined because the function ℓ is invariant under multiplication by a unit. Now define
M (R)• ϕ : n n M −→( On rk M)ωn + ρ (M). 7−→ − n Notice that if A, B M (R)• satisfy A B, then ρ (A) = ρ (B) and ∈ n ∼ n n ϕn(A)= ϕn(B).
Proposition 3.2. The function ϕn is a right and left Euclidean stathm. The remainder of this section will be devoted to the proof of Proposi- tion 3.2. This will imply Theorem 3.1. T Remark 3.3. For any n> 1, ϕn = ϕn, so the proof of Proposition 2.4 implies that ϕn is a left Euclidean stathm if and only if it is a right Euclidean stathm. Remark. The Euclidean stathm is in no way unique. For instance, the fol- lowing function is a left and right Euclidean stathm:
M2 (Z)• −→ O det M if det m = 0, ψ : 2 M | | α6 0 7−→ ω + α if M . | | ∼ 0 0 See Proposition 6.5 for a more general construction. 3.2. Case of size 2 matrices. To prove Proposition 3.2, we will first deal with 2 by 2 matrices. Lemma 3.4. Let A = a 0 ,B = diag(b , b ) M (R), where b divides b c 1 2 ∈ 2 1 b2 = 0. If b1 divides a, b and c, but b2 does not divide b or b2 does not divide c, then6 there exists Q M (R) such that ∈ 2 A BQ diag(b , e), − ∼ 1 where e R is such that b e and e is a strict divisor of b . ∈ • 1| 2 Proof of Lemma 3.4. Set e = gcd(b, c, b2) = 0, which is a multiple of b1 and a strict divisor of b . Lemma 2.7 implies that6 there exists z,t R such that 2 ∈ gcd (c + b2z, b + b2t)= e. Therefore, there exist λ, µ R which are coprime and satisfy ∈ (6) λ(b + b2t)+ µ(c + b2z)= e.
a/b1 µ λ Set Q = t− z . Then − − µb λb A BQ = 1 − 1 . − b + b2t c + b2z Since λ and µ are coprime, the gcd of the coefficients of A BQ is b1. Besides, (6) implies that det(A BQ)= b e. As a result, thanks− to Lemma 2.3(b), − 1 A BQ diag(b , e). − ∼ 1 10 PIERRE LEZOWSKI
Lemma 3.5. Let A M2 (R) and B = diag(b1, b2) M2 (R), where b1 b2 and b = 0. Then there∈ exists Q M (R) such that ∈ | 2 6 ∈ 2 A = BQ or (rk(A BQ) = 2 and ρ (A BQ) < ρ (B)) . − 2 − 2 a 0 Proof of Lemma 3.5. Take T GL2 (R) such that AT = b c . We distin- guish two cases. ∈ a. Assume that b1 does not divide a, b or c. Fix λ, µ 0, 1 such that λ 1 ∈ { } Qλ,µ = 0 −µ satisfies a b1λ b1 det(AT BQλ,µ)= − µb2(a b1λ) = 0. − b c − − 6
Therefore, AT BQλ,µ diag( α, β) for α β R, β = 0. But − ∼ | ∈ 6 a b1λ b1 (7) AT BQλ,µ = − , − b c µb2 − so, thanks to Lemma 2.3(a), α = gcd(a b λ, b , b, c µb ) = gcd(a, b , b, c), − 1 1 − 2 1 since b1 divides b2. Then α is a strict divisor of b1. In particular, ℓ(α) <ℓ(b1). Consequently, by setting Q = Q T 1, we have rk(A BQ) = 2 and λ,µ − − ρ (A BQ)= ℓ(α)ω + ℓ(β) <ℓ(b )ω + ℓ(b )= ρ (B). 2 − 1 2 2 b. If b1 divides a, b, and c, we have two sub-cases. Either b2 divides b and 1 1 c and then Q = (B− AT )T − M2 (R) satisfies A = BQ, or b2 does not divide b or c, and then we apply∈ Lemma 3.4 to find Q M (R) such that ∈ 2 A BQ diag(b , e), − ∼ 1 where e R is such that b e and e is a strict divisor of b . Then ∈ • 1| 2 ρ (A BQ)= ℓ(b )ω + ℓ(e) <ℓ(b )ω + ℓ(b )= ρ (B). 2 − 1 1 2 2
3.3. Case of size n full-rank matrices. Now, we extend Lemma 3.5 to n by n matrices, where n Z . ∈ >1 Lemma 3.6. Let n Z , A M (R), and B = diag(b , . . . , b ) ∈ >1 ∈ n 1 n ∈ Mn (R), where b b . . . b = 0. 1| 2| | n 6 Then there exists Q M (R) such that ∈ n A = BQ or (rk(A BQ)= n and ρ (A BQ) < ρ (B)) . − n − n Proof of Lemma 3.6. We prove it by induction on n 2. The case n = 2 is Lemma 3.5. Set n 3 and assume that Lemma 3.6≥ holds for all strictly smaller dimensions. ≥ Take A Mn (R). Consider T GLn (R) such that AT = (ai,j)1 i,j n is lower triangular.∈ ∈ ≤ ≤ 1st step. Assume that there exist 1 i , j n such that b does not divide ≤ 0 0 ≤ 1 ai0,j0 . For any 1 i n, take µi 0, 1 such that ai,i µibi = 0. Then take λ 0, 1 such≤ that≤ ∈ { } − 6 ∈ { } (a µ b )(a µ b ) b (a λb ) = 0. 1,1 − 1 1 2,2 − 2 2 − 1 2,1 − 2 6 ON SOME EUCLIDEAN PROPERTIES OF MATRIX ALGEBRAS 11
Now consider
µ1 1 − 02,n 2 Q′ = λ µ2 − M (R) , ∈ n 0n 2,2 diag(µ3,...,µn) − then AT BQ′ has rank n, its invariant factors are b1′ . . . bn′ = 0, and b1′ divides all− coefficients of AT BQ (cf. Lemma 2.3(a)).| In| particular,6 b − ′ 1′ divides b1 and ai0,j0 , so it is a strict divisor of b1. Therefore, ℓ(b1′ ) < ℓ(b1). Taking Q = Q T 1, we find rk(A BQ)= n and ′ − − n n n i n i ρ (A BQ)= ρ (AT BQ′)= ω − ℓ(b′ ) < ω − ℓ(b )= ρ (B). n − n − i i n Xi=1 Xi=1 2nd step. From now on, we assume that b1 divides all coefficients of AT . Take A′ Mn 1 (R) such that ∈ − b1 01,n 1 − a1,1 1 01,n 1 a2,1 AT B b1 − − = . . − 0n 1,1 0n 1 . A′ − − an,1 Set B = diag(b , . . . , b ). By the induction hypothesis, there exists Q ′ 2 n ′ ∈ Mn 1 (R) such that R′ = A′ B′Q′ satisfies − − R′ = 0n 1 or rk R′ = n 1 and ρn 1(R′) < ρn 1(B′) . − − − −
1. Assume that R′ = 0n 1 . Its invariant factors (b2′ , . . . , bn′ ) are all divisible 6 − by b1, as all coefficients of A′ and B′ are divisible by b1 (see Lemma 2.3(a)). There exist X,Y GLn 1 (R) such that YR′X = diag(b2′ , . . . , bn′ ). Then ∈ − b1 01,n 1 − a2,1 0 01,n 1 . B − . A′ − 0n 1,1 Q′ − an,1 b1 01,n 1 1 − 1 1 01,n 1 − a2,1 1 01,n 1 − = − . − 0n 1,1 Y . YR′X 0n 1,1 X − − an,1 diag(b , b′ , . . . , b′ ). ∼ 1 2 n
Thus, there exists Q′′ Mn (R) such that AT BQ′′ diag(b1, b2′ , . . . , bn′ ). Taking Q = Q T 1 ∈M (R), we obtain − ∼ ′′ − ∈ n A BQ diag(b , b′ , . . . , b′ ). − ∼ 1 2 n
Since b1 b2′ . . . bn′ = 0, rk(A BQ) = n and they are the invariant factors of A BQ| |. Hence| 6 − − n 1 n 1 ρn(A BQ)= ℓ(b1)ω − + ρn 1(R′) <ℓ(b1)ω − + ρn 1(B′)= ρn(B). − − − 12 PIERRE LEZOWSKI
2. Assume now that R′ = 0n 1. We distinguish two subcases. − 2.a. First, assume that for all l> 1, bl divides al,1, then
b1 01,n 1 1 01,n 1 − − a2,1 a2,1/b2 . B . = 0n. . A′ − . Q′ an,1 an,1/bn Therefore, there exists Q′′ Mn (R) such that AT BQ′′ = 0n. Take 1 ∈ − Q = Q′′T − , then A BQ = 0 . − n 2.b. Now assume that there exists l > 1 such that bl does not divide al,1. Define b1 01,n 1 − a2,1 1 01,n 1 A′′ = . B − + B . A′ − 0n 1,1 Q′ − an,1 b1 01,n 1 − a2,1 = . . . diag(b2, . . . , bn) an,1 Take l> 1 to be the smallest integer such that there exists m l such that ≥ bl does not divide am,1. 2.b.i. If l
Now, we have all the tools required to prove Proposition 3.2. 14 PIERRE LEZOWSKI
3.5. Proof of Proposition 3.2. Recall that R is a PID and n Z . ∈ >1 Thanks to Remark 3.3, if suffices to prove that ϕn is a right Euclidean stathm. Let A, B Mn (R), B = 0. We want to find Q Mn (R) such that A = BQ or ϕ (A ∈ BQ) < ϕ (B6 ). ∈ n − n Set r = rk B. We take X,Y,T GLn (R), and b1 b2 . . . br R• such that ∈ | | | ∈ Y BX = diag(b , . . . , b , 0,..., 0) M (R) , 1 r ∈ n and Y AT = (ai,j)1 i,j n is lower triangular. ≤ ≤ 1. If r = n, then, thanks to Lemma 3.6, there exists Q′ Mn (R) such that Y AT = YBXQ , or rk(Y AT YBXQ ) = n and ρ (Y∈ AT YBXQ ) < ′ − ′ n − ′ ρn(Y BX). Setting Q = XQ′, we have as required A = BQ or ϕ (A BQ) < ϕ (B). n − n 2. From now on, we assume that r A1 0r,n r Y AT YBXD = − − A2 A3 where A1 Mr (R) is lower triangular, A2 Mn r,r (R), A3 Mn r (R). ∈ ∈ − ∈ − By construction, rk A1 = r. M (1) M (2) Notation. Let M = , where 1 r rk(Y AT YBXQ′) rk Extr (Y AT YBXQ′; i , 1) > r. − ≥ r − 0 It follows that rk(A BQ) > r = rk B, for Q = XQ T 1, which implies − ′ − ϕ (A BQ) < ϕ (B). n − n 3. Now we can assume that A2 = 0n r,r and A3 = 0r. If r = 1, then we − can apply Lemma 3.7 to find Q′ Mn (R) such that Y AT = YBXQ′ or ϕ (Y AT YBXQ ) < ϕ (Y BX).∈ Set Q = XQ T 1, then n − ′ n ′ − A = BQ or ϕ (A BQ) < ϕ (B). n − n It remains to consider r> 1. We set B = diag(b , . . . , b ) M (R). Thanks 1 1 r ∈ r to Lemma 3.6, there exists Q1 Mr (R) such that for R1 = A1 B1Q1, we have ∈ − R1 = 0r or (rk R1 = r and ρr(R1) < ρr(B1)) . In any case, set Q1 0r,n r 1 Q = X − + D T − , 0n r,r 0n r − − then 1 R1 0r,n r 1 A BQ = Y − − T − . − 0n r,r 0n r − − If R = 0 , then A = BQ. If R = 0 , then rk R = r, so 1 r 1 6 r 1 ϕ (A BQ) = (n r)ωn + ρ (R ) < (n r)ωn + ρ (B )= ϕ (B). n − − r 1 − r 1 n 4. The Euclidean property for matrix algebras over a PIR The aim of this section will be to prove the following property. Theorem 4.1. Let R be a commutative ring and n Z>1. Then Mn (R) is right and left Euclidean if and only if R is a principal∈ ideal ring. 4.1. Some general properties of Euclidean and principal ideal rings. A PIR R is said to be special if R has a unique prime ideal and this ideal is nilpotent. To infer Theorem 4.1 from Theorem 3.1, we will use the following property due to Samuel and Zariski. Proposition 4.2 ([ZS75, Theorem 33, p. 245]). Let R be a PIR, then it can l be written as a direct product R = i=1 Ri, where each Ri is a PID or a special PIR. Q Special PIRs can be easily dealt with. Lemma 4.3. Let S be a special PIR. Then for any n Z 1, Mn (S) is ∈ ≥ Euclidean, and e(Mn (S)) <ω. 16 PIERRE LEZOWSKI Proof. Let p be the prime ideal of S. Take m Z such that pm = 0. We ∈ >0 Mn (S)• 0, 1,...,m 1 prove that the function ψn : −→ { − } is a right M vp det(M) Euclidean stathm. 7−→ ◦ We know that S is the homomorphic image of a PID R (see [Hun68]). Consider a surjective homomorphism π : R S, which we extend to a surjective homomorphism π : M (R) M −→(S). n −→ n Take A, B M (S), B = 0 . There exist Aˆ, Bˆ M (R) such that ∈ n 6 n ∈ n π(Aˆ)= A and π(Bˆ)= B. Besides, there exist Dˆ , Aˆ , Bˆ M (R) such that ′ ′ ∈ n Aˆ = Dˆ Aˆ′, Bˆ = DˆBˆ′, and Aˆ′Mn (R)+ Bˆ′Mn (R) = Mn (R). Set D = π(Dˆ), A′ = π(Aˆ′), and B′ = π(Bˆ′). Then we have A = DA′, B = DB′, and A′Mn (S)+ B′Mn (S) = Mn (S) . But sr S = 1, so sr Mn (S) = 1, and there exists Q′ Mn (S) such that U = A B Q GL (S). ∈ ′ − ′ ′ ∈ n If vp det(B ) > 0, set Q = Q , then ◦ ′ ′ vp det(A BQ)= vp det(DU)= vp det D < vp det(DB′)= vp det B. ◦ − ◦ ◦ ◦ ◦ 1 If vp det(B ) = 0, then B GL (S). Set Q = B A , we have A BQ = ◦ ′ ′ ∈ n ′− ′ − 0n. Lemma 4.4. Let l Z 1 and Ai, 1 i l be right Euclidean rings. Then l ∈ ≥ ≤ ≤ the product ring i=1 Ai is a right Euclidean ring. Proof. You can referQ to [Sam71, Proposition 6] or [Cla15, Theorem 3.13], where the commutativity hypothesis is not used. We give some details to get some insight into Remark 4.5. Note that an immediate induction shows that it is enough to consider the product of two right Euclidean rings A1 and A2. In that case, we consider two right Euclidean stathms ϕi : Ai• , extended at 0 as in (5), for i = 1, 2, and we can prove that −→ O (A A ) ϕ : 1 2 • (r×, r ) −→ϕ O (r ) ϕ (r ) 1 2 7−→ 1 1 ⊕ 2 2 is a right Euclidean stathm. In fact, we also have the following property. Remark 4.5 ([Cla15, Theorem 3.406]). With the above hypotheses, l l l e(Ai) e Ai e(Ai). ≤ ! ≤ Xi=1 Yi=1 Mi=1 The upper bound, which is easily implied by the proof of Lemma 4.4 above, has the following consequence: if for any i, Ai is right Euclidean and e(Ai) < l l ω), then i=1 Ai is also right Euclidean and e i=1 Ai <ω. Q Q 6 As in Footnote 5, Clark sets himself in a commutative context, but this property does not rely on the commutative hypothesis. ON SOME EUCLIDEAN PROPERTIES OF MATRIX ALGEBRAS 17 4.2. Proof of Theorem 4.1. Let R be a commutative ring and n Z . ∈ >1 If Mn (R) is right Euclidean, then every right ideal of Mn (R) is principal. Therefore, thanks to Proposition 2.5, R is a PIR. Conversely, assume that R is a PIR. Thanks to Proposition 4.2, R can be written as R = l R , such that for any 1 i l, R is either a PID or i=1 i ≤ ≤ i a special PIR. But now, M (R) is isomorphic to l M (R ). Thanks to Q n i=1 n i Theorem 3.1 and Lemma 4.3, for any 1 i l, Mn (Ri) is Euclidean. Using l ≤ ≤ Q Lemma 4.4, we can conclude that i=1 Mn (Ri) is Euclidean. 5. EffectivityQ 5.1. General result. The fact that a ring A is right Euclidean for some right Euclidean stathm ϕ : A does not mean that we know how to • −→ O compute a quotient (or equivalently a remainder) for each pair (a, b) A A•, that is to say an element q A such that ∈ × ∈ a = bq or ϕ(a bq) < ϕ(b). − Nevertheless, in the case when A = Mn (R), with a further condition on R, we can effectively compute it. More precisely, we have the following property. Proposition 5.1. Let n> 1 and R be a PID. The following statements are equivalent. (a) For any a, b R, we can compute d = gcd(a, b), and elements u, v R such that∈ au + bv = d. ∈ (b) For any A, B Mn (R), we can compute some Q Mn (R) such that ϕ (A BQ∈) < ϕ (B). ∈ n − n Proof. Assuming (a), a careful reading of the proof in Section 3 shows that we may compute (b). Indeed, all constructions rely on gcds (see Remark 2.8), and reduction of matrices into echelon form or Smith normal form. These reductions can be explicitly computed assuming (a). Conversely, take a, b R. Set A = diag(1,..., 1, a) Mn (R) and B = diag(1,..., 1, b) M (R∈). As we can compute quotients∈ (and remainders) ∈ n of Euclidean divisions in Mn (R), we can apply the Euclidean algorithm to A and B: Algorithm 5.2. Input: A, B M (R). ∈ n Set D = A, U = 1n, V = 0n, D1 = B, U1 = 0n, V1 = 1n. (i) If D1 = 0n, return [U, V, D]. (ii) Compute Q, R M (R) such that D = D Q + R. Set ∈ n 1 (D, U, V, D ,U ,V ) = (D ,U ,V ,R,U U Q, V V Q) 1 1 1 1 1 1 − 1 − 1 and go to Step (i). We obtain U, V M (R) such that ∈ n (8) AU + BV = D, where D is a gcld of A and B in Mn (R). But a gcld of A and B is diag(1,..., 1, d) where d is a gcd of a and b. Therefore, there exists T ∈ GLn (R) such that D = diag(1,..., 1, d)T . In particular, det D = dε for 18 PIERRE LEZOWSKI some ǫ R×. We replace d by dε such that det D = d. Then, denoting by C the cofactor∈ matrix of D, (8) implies that T T AUC + BVC = dIn. Identifying the coefficients in position (n,n), we obtain au + bv = d. 5.2. Computation of gcd of matrices, an example. It is straightfor- ward to remark that Algorithm 5.2 above will return a gcld of A and B in a finite number of steps. We can similarly compute gcrds. The following example illustrates such computations, and the fact that gcrd and gcld may be unconnected. Example 5.3. In M3 (Z), consider the matrices 1 1 0 1 1 1 A = −2 2 0 and B = −2− 2− 2 . 1 −1 2 2 1 0 − − 2 0 2 Then A = B −3 1 2 , therefore a gcld for A and B is B, which has 0− 0− 0 rank 2. Besides, 0 0 0 1 1 0 A = 0 0 1 B + −4 1 0 , 0 0− 0 1 −1 2 − − 2 1 1 1 1 0 2 1 3 − − − − B = 7 3 4 4 1 0 + 3 2 10 , −0− 0− 0 1 −1 2 2 1 0 − − 1 1 0 10 3 5 2 1 3 −4 1 0 = 20 6 9 −3− 2− 10 . 1 −1 2 − 4 −1 −3 2 1 0 − − − − − 2 1 3 Consequently, a gcrd of A and B is −3− 2− 10 , which has rank 3. 2 1 0 Remark 5.4. In a PID R, under the effectivity conditions of Proposition 5.1, there is a more direct way to compute gclds and gcrds of matrices in Mn (R). Indeed, given A, B Mn (R), we can proceed as follows. Consider the matrix A B M (R∈) and compute R M (R), U GL (R) such that ∈ n,2n ∈ n ∈ 2n