Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
arXiv:1705.06122v1 [math.NT] 17 May 2017 Multidimensional p-adic continued fraction algorithms Asaki Saito1 , Jun-ichi Tamura2 , and Shin-ichi Yasutomi3, 1 ∗ Future University Hakodate, Hakodate, Hokkaido 041-8655, Japan 2 Tsuda College, Kodaira, Tokyo 187-8577, Japan 3 Toho University, Funabashi, Chiba 274-8510, Japan May 18, 2017 Abstract We give a new class of multidimensional p-adic continued fraction algorithms. We propose an algorithm in the class for which we can expect that multidimensional p-adic version of Lagrange’s Theorem holds. 1 Introduction Throughout the paper, p denotes a fixed prime number, Qp (resp., Zp ) the closure of Q(resp., Z) with respect to the p-adic topology. For x ∈ Qp , ordp (x) denotes the highest power of p by which x is divided to be a p-adic integer. Schneider [11] has introduced the following p-adic continued fraction algorithm. Let ξ ∈ Zp . We define ξ1 ∈ pZp by ξ1 = ξ − a0 , where a0 ∈ {0, 1, . . . , p − 1}. We define ξn (n ≥ 2) recursively by ξn = pordp (ξn−1 ) − an−1 , ξn−1 where an−1 ∈ {1, . . . , p − 1} is chosen such that ξn ∈ pZp . Then, we have pordp (ξ1 ) ξ = a0 + pordp(ξ2 ) a1 + a2 + ∗ Author to whom correspondence mail:[email protected] pordp (ξ3 ) a3 + . . . should 1 . be addressed. Electronic Weger [13] has shown that some quadratic elements have not eventually periodic expansion by Schneider’s algorithm. Ruban [7] has proposed another p-adic continued fraction algorithm different from Schneider’s. Ooto [6] has shown a similar result concerning the algorithm given by Ruban. Although Browkin [2] proposed some p-adic continued fraction algorithms, it has not been proved that the continued fraction expansion of every quadratic element obtained by his algorithm is eventually periodic. We [8] have introduced some new p-adic continued fraction algorithms, and have shown p-adic version of Lagrange’s theorem, i.e., if α ∈ Qp is a quadratic element over Q, the continued fractions for α obtained by our algorithms become periodic. Bekki [1] has shown p-adic version of Lagrange’s theorem for imaginary irrationals on his continued fraction algorithm. It seems that there are only a few results on multidimensional p-adic continued fractions. By disclosing a link between the hermitian canonical forms of certain integral matrices and p-adic numbers, Tamura [12] has showed that a multidimensional p-adic continued fraction converges to (x, x2 , . . . , xn−1 ) in the p-adic sense without considering algorithms of continued fraction expansion, where x is the root of some polynomial with degree n. We [9] considered some new class of multidimensional p-adic continued fractions and constructed some explicit formulae of multidimensional continued fractions related to algebraic elements of Qp over Q. In this paper, we propose a class A of multidimensional p-adic continued fraction algorithms such that we can expect that Conjecture. For any Q-basis {1, α1 , . . . , αs } of any given field K ⊂ Qp with degree [K : Q] = s + 1, the continued fraction expansion of α = (α1 , . . . , αs ) by the s-dimensional continued fraction algorithm in the class A always becomes eventually periodic. In particular we can give a proof of the conjecture for s = 1. For s > 1, we show that there exist infinitely many α ∈ K s having eventually periodic continued fraction expansion of by the algorithm. We give numerical experiments for s = 2, 3, 4, 5, which support the conjecture, cf. Tables 9-12 in marked contrast with Tables 1-8 concerning some other algorithm. 2 Notations and some lemmas We denote by Qp the algebraic closure of Qp , and by Ap the set of algebraic elements over Q in Qp . We put Vp := Zp \pZp , C := {0, 1, . . . , p − 1}. For α= X cn pn ∈ Qp \{0} (cn ∈ C), n∈Z we define ordp (α) := min{n|cn 6= 0} (ordp (0) := ∞), |α|p := p−ordp (α) (|0|p = 0) , ωp (α) := c0 and ⌊α⌋ := Σn∈Z≤0 cn pn . 2 For an integer s > 0 and α = (α1 , . . . , αs )T ∈ Qsp , we define ordp (α) := min ordp (αi ), 1≤i≤s |α|p := p−ordp (α) , ⌊α⌋ := (⌊α1 ⌋, . . . , ⌊αs ⌋)T . We define a transformation Tb on Qp as follows: for α ∈ Qp with α 6= 0, Tb (α) := pordp(α) − ωp α Tb (0) := 0. pordp(α) α , The following Condition (H) will play a key role along the paper. Condition (H) For β ∈ (Ap \Q) ∩ pZp , if β has a minimal polynomial xn + a1 xn−1 + . . . + an over Q such that ai ∈ Zp for 1 ≤ i ≤ n, ordp (an−1 ) = 0 and ordp (an ) > 0, then we say that β satisfies the condition (H). pm x+a for x ∈ pZp , where m is a positive integer and a ∈ Vp . Then, for α, β ∈ pZp , |T (α) − T (β)|p = p−m |α − β|p . Lemma 2.1. Let T be the transformation on pZp defined by T (x) := Proof. For α, β ∈ pZp , pm (α − β) = p−m |α − β|p . |T (α) − T (β)|p = (α + a)(β + a) p Lemma 2.2. Let β ∈ (Ap \Q) ∩ pZp . Then, there exists a positive integer m such that for every integer n ≥ m, Tbn (β) satisfies the condition (H). Proof. First, we assume that every algebraic conjugate of β is in Qp . We denote by β1 (= β), . . . , βn the algebraic conjugates of β and by σ1 (= identity), . . . , σn the embeddings of Q(β) into Qp . Then, we have for 1 ≤ i ≤ n pordp (β) σi (Tb (β)) := − ωp σi (β) pordp(β) β , which implies that for 2 ≤ i ≤ n, ordp (σi (Tb (β))) = 0 if ordp (σi (β)) < ordp (β), ordp (σi (Tb (β))) < 0 if ordp (σi (β)) > ordp (β), and ordp (σi (Tb (β))) ≥ 0 if ordp (σi (β)) = ordp (β). Therefore, for 2 ≤ i ≤ n if for some integer m > 0 ordp (σi (Tbm (β))) 6= ordp (Tbm (β)) holds, then ordp (σi (Tbk (β))) = 0 for k ≥ m+2. We assume that there exists some j with 2 ≤ j ≤ n such that ordp (σj (Tbm (β))) = 3 ordp (Tbm (β)) holds for every integer m ≥ 0. Then, it is not difficult to see that for every integer m ≥ 0 ordp(T m ((σj (β))) ordp (T m (β)) b b p p ωp = ω . p m m Tb (σj (β)) Tb (β) For m ∈ Z≥0 we set a transformation Tm on pZp by for x ∈ pZp m pordp(Tb (β)) Tm (x) := ordp (T m (β)) . b p x + ωp m Tb (β) By virtue of Lemma 2.1 we see that |β − σj (β)|p =|T0 ◦ · · · ◦ Tm−1 (Tbm (β)) − T0 ◦ · · · ◦ Tm−1 (Tbm (σj (β)))|p m−1 = p−Σk=0 ordp (Tbk (β)) m−1 |Tbm (β) − Tbm (σj (β))|p < p−Σk=0 ordp (Tbk (β)) . Therefore, as m → ∞ we have |β − σj (β)|p = 0, which contradicts β 6= σj (β). Thus, we see that there exists some integer m′ > 0 such that for every integer ′′ ′′ m′′ > m′ ordp (σi (Tbm (β))) = 0 for i = 2, . . . , n. Therefore, Tbm (β) satisfies the condition (H) for m′′ > m′ . If some algebraic conjugates of β are not included in Qp , then considering a finite extension of Qp adjoining all conjugates of β we have a similar proof. We remark that for a polynomial p(x) = xn + a1 xn−1 + . . . + an ∈ Q[x] with ai ∈ Zp for 1 ≤ i ≤ n, ordp (an−1 ) = 0 and ordp (an ) > 0, by Hensel’s lemma there exists a α ∈ Qp such that p(α) = 0 and ordp (α) = ordp (an ). By Lemma 2.2 we have the following proposition. Proposition 2.3. Let K ⊂ Qp be a finite extension of Q. Then, there exists α ∈ K which satisfies the condition (H) and K = Q(α). 3 c-map In this section, we introduce a class of multidimensional p-adic continued fraction algorithms. Let K ⊂ Qp be a finite extension of Q and its degree over Q be deg. We define s ∈ Z by s := deg − 1 for K 6= Q and s := 1 for K = Q. Ind denotes the set {1, 2, . . . , s}. D denotes K s and E denotes (pZp )s ∩ K s . We denote L(D) the set of linear fractional transformations on D. We now introduce a map Φ : D → Ind × L(D) × GL(s, Zp ∩ Q). For α = (α1 , . . . , αs ) ∈ D, Φ(α) is defined by Φ(α) := (φ(α), Fα , Aα ), where Fα = (f1 , . . . , fs ) is given as follows: In the case of αφ(α) 6= 0: We define fi for (x1 , . . . , xs ) ∈ D with xφ(α) 6= 0 by ui pordp(αi ) − vi if i = φ(α), xi fi (x1 , . . . , xs ) := u′ pk xi − vi′ if i 6= φ(α), i xφ(α) 4 where ui , vi , u′i ∈ Vp ∩ Q, vi′ ∈ Zp ∩ Q, and k = max{ordp (αφ(α) ) − ordp (αi ),0}. We also assume that fi (α) ∈ pZp for 1 ≤ i ≤ s. We do not define fi for (x1 , . . . , xs ) ∈ D with xφ(α) = 0. In the case of αφ(α) = 0: We define fi for (x1 , . . . , xs ) ∈ D by fi (x1 , . . . , xs ) := xi . In what follows, we call Φ a c-map. Lemma 3.1. Let Φ: Φ(α) = (φ(α), Fα , Aα ) for α = (α1 , . . . , αs )T ∈ D be a c-map and α′ ∈ D with αφ(α′ ) 6= 0. Then, Aα′ (f1 (α′ ), . . . , fs (α′ ))T ∈ E holds, where F = (f1 , . . . , fs )T . Proof. The proof is easy. Lemma 3.2. Let Φ: Φ(α) = (φ(α), Fα , Aα ) for α ∈ D be a c-map, α′ = (α1 , . . . , αs )T ∈ D and (β1 , . . . , βs )T = Aα′ (f1 (α′ ), . . . , fs (α′ ))T , where F = (f1 , . . . , fs )T . Let 1, α1 , . . . , αs be linearly independent over Q. Then, 1, β1 , . . . , βs are linearly independent over Q, Proof. Let 1, α1 , . . . , αs be linearly independent over Q. It is not difficult to see that 1, f1 (α′ ), . . . , fs (α′ ) are linearly independent over Q. Since Aα′ ∈ GL(s, Zp ∩ Q), 1, β1 , . . . , βs are linearly independent over Q. Lemma 3.3. Let Φ : Φ(α) = (φ(α), Fα , Aα ) for α ∈ D be a c-map. If αφ(α) 6= 0 for α = (α1 , . . . , αs )T ∈ D, then Fα−1 : E → D exists and for (x1 , . . . , xs ) ∈ E Fα−1 (x1 , . . . , xs )T is given by as follows: Fα−1 (x1 , . . . , xs )T = (g1 (x1 , . . . , xs ), . . . , gs (x1 , . . . , xs ))T , where if i = φ(α), for some u′i , vi′ ∈ Vp ∩ Q gi (x1 , . . . , xs ) = u′i pordp(αi ) , xi + vi′ if i 6= φ(α), for some u′i , vi′ ∈ Vp ∩ Q, wi′ ∈ Zp ∩ Q gi (x1 , . . . , xs ) = u′i pk (xi + wi′ ) , ′ xφ(α) + vφ(α) where k = min{ordp (αφ(α) ), ordp (αi )}. Proof. The proof is easy. Lemma 3.4. Let Φ : Φ(α) = (φ(α), Fα , Aα ) for α ∈ D be a c-map and α′ = (α1 , . . . , αs )T ∈ E. Then, Fα−1 (A−1 α (E)) ⊂ E holds. 5 −1 Proof. It is not difficult to see A−1 α (E) = E. We will show Fα (E) ⊂ E. If −1 αφ(α) = 0, then F is the identity map and we have Fα (E) ⊂ E. We suppose αφ(α) 6= 0. Let (x1 , . . . , xs ) ∈ E and Fα−1 (x1 , . . . , xs )T = (g1 (x1 , . . . , xs ), . . . , gs (x1 , . . . , xs ))T . By Lemma 3.4 we have for i = φ(α), ′ ord (α ) u p p i = p−ordp(αi ) , |gi (x1 , . . . , xs )|p = i xi + vi′ p where u′i , vi′ ∈ Vp ∩ Q, and for i 6= φ(α) u′ pk (x + w′ ) i i i |gi (x1 , . . . , xs )|p = ≤ p−k , ′ xφ(α) + vφ(α) p where k = min{ordp (αφ(α) ), ordp (αi )}, u′i , vi′ ∈ Vp ∩ Q and wi′ ∈ Zp ∩ Q. Lemma 3.5. Let Φ : Φ(α) = (φ(α), Fα , Aα ) for α ∈ D be a c-map. If αφ(α) 6= 0 for α = (α1 , . . . , αs )T ∈ E, then for x = (x1 , . . . , xs )T , y = (y1 , . . . , ys )T ∈ E −1 −j (x) − Fα−1 A−1 |Fα−1 Aα α (y)|p ≤ p |x − y|p , where j = min{ordp (αi )|1 ≤ i ≤ s}. If αφ(α) = 0, for x = (x1 , . . . , xs )T , y = (y1 , . . . , ys )T ∈ E |Fα−1 A−1 α (x) − Fα−1 A−1 (y)| = |x − y| . p p α Proof. Let α = (α1 , . . . , αs )T ∈ E. We assume that x = (x1 , . . . , xs )T , y = (y1 , . . . , ys )T ∈ E. First, we suppose αφ(α) 6= 0. Let Fα−1 = (g1 , . . . , gs )T . By Lemma 3.4 we see that for i = φ(α), ′ ord (α ) ui p p i u′i pordp (αi ) − |gi (x) − gi (y)|p = xi + vi′ yi + vi′ p = |u′i pordp (αi ) |p |xi − yi |p = p−ordp(αi ) |xi − yi |p , |xi + vi′ |p |yi + vi′ |p where u′i , vi′ ∈ Vp ∩ Q for i 6= φ(α) and 1 ≤ i ≤ s, u′ pk (x + w′ ) u′ pk (y + w′ ) i i i i i i |gi (x) − gi (y)|p = − ′ ′ xφ(α) + vφ(α) yφ(α) + vφ(α) p = |u′i pk |p |(xi − ′ yi )(vφ(α) + yφ(α) ) − (xφ(α) − yφ(α) )(wi′ + yi )|p ′ ′ |p |yφ(α) + vφ(α) |p |xφ(α) + vφ(α) ≤ p−k |x − y|p , where k = min{ordp (αφ(α) ), ordp (αi )}, u′i , vi′ ∈ Vp ∩ Q and wi′ ∈ Zp ∩ Q. Since −1 −1 −1 −1 −1 A−1 α (x), Aα (y) ∈ E and |Aα (x) − Aα (y)|p = |x − y|p , we have |Fα Aα (x) − −1 −1 −j Fα Aα (y)|p ≤ p |x − y|p , where j = min{ordp (αi )|1 ≤ i ≤ s}. Secondly, we suppose αφ(α) = 0. Then, F is the identity map. Therefore, we have −1 −1 |Fα−1 Aα (x) − Fα−1 ′ Aα (y)|p = |x − y|p . 6 We introduce a p-adic continued fraction via c-map. Let Φ : Φ(α) = (φ(α), Fα , Aα ) for α = (α1 , . . . , αs )T ∈ D be a c-map. Let α = (α1 , . . . , αs )T ∈ D. We set α(0) := α. α(1) , α(2) , . . . are defined inductively as follows: We suppose that α(n) for n ∈ Z≥0 is defined. We set α(n+1) := Aα(n) Fα(n) (α(n) ). We say that α has a Φ continued fraction expansion {Φ(α(0) ), Φ(α(1) ), . . .}. (n) (n) α(n) = (α1 , . . . , αs ) is referred to the n-th remainder of α. For an integer n ≥ 0 we define the n-th convergent π(α; n) by −1 −1 −1 π(α; n) := Fα−1 (0) Aα(0) · · · Fα(n−1) Aα(n−1) (0). We remark that π(α; n) is a rational number for n ∈ Z≥0 . We say that α has a periodic Φ continued fraction expansion if α(m1 ) = α(m2 ) holds for some m1 , m2 ∈ Z≥0 with m1 6= m2 . We define that α has a finite Φ continued fraction expansion if α(m) = 0 holds for some m ∈ Z≥0 . We define that α has an infinite Φ continued fraction expansion if α does not have a finite Φ continued fraction expansion. Theorem 3.6. Let Φ : Φ(α) = (φ(α), Fα , Aα ) for α ∈ D be a c-map and α′ ∈ D. Let φ(α′ (n) ) be not equal to 0 for infinitely many n. Then, lim π(α′ ; n) = α′ . n→∞ Proof. By Lemma 3.1 there exists an integer m ≥ 0 such that α′ By Lemma 3.5 we have |α′ (m) (m) ∈ E holds. (m) − π(α′ ; n)|p −1 −1 −1 ′ (m+n) )− = F −1 (m) A ′ (m) · · · F ′ (m+n−1) A ′ (m+n−1) (α α′ α α α A−1 · · · F −1 A−1 (0) F −1 ′ (m) ′ (m) ′ (m+n−1) ′ (m+n−1) α ≤p α α −(j0 +...+jn−1 ) α ′ (m+n) |α |p < p p −(j0 +...+jn−1 ) , (m+i) (m+i) where for i ∈ Z≥0 ji := min{ordp (αk )|1 ≤ k ≤ s} if α ′ (m+i) 6= 0 and φ(α ) P∞ (m) (m) (m+i) ′ ; n) = α′ . ji := 0 if α ′ (m+i) = 0. Since k=0 jk = ∞, we have lim π(α φ(α ) n→∞ Therefore, we have lim π(α′ (0) n→∞ ; m + n) A−1 · · · F −1 A−1 (π(α′ = lim F −1 ′ (0) ′ (0) ′ (m−1) ′ (m−1) n→∞ α α α = F −1 A−1 · · · F −1 A−1 (α′ ′ (0) ′ (0) ′ (m−1) ′ (m−1) α α α (m) α (m) α ; n)) ) = α′ . Proposition 3.7. Let Φ : Φ(α) = (φ(α), Fα , Aα ) for α ∈ D be a c-map. If for α ∈ D, α has a finite Φ continued fraction, then lim π(α; n) = α. n→∞ 7 Proof. Let α ∈ D and α have a finite Φ continued fraction. Then, there exists an integer m ≥ 0 such that α(m) = 0. It is clear that α(n) = 0 for every n ≥ m. −1 −1 −1 Therefore, for every n ≥ m we see that α(0) = Fα−1 (0) Aα(0) · · · Fα(n−1) Aα(n−1) (0) = π(α; n). Thus, we obtain the proposition. 4 Algorithms In this section, we introduce some p-adic continued fraction algorithms. Let K ⊂ Qp be a finite extension of Q and its degree over Q be deg. We define s ∈ Z as s := deg − 1 for K 6= Q and s := 1 for K = Q. By Proposition 2.3 there exists an element z in K which satisfies the condition (H) and K = Q(z). Hereafter, we suppose that z ∈ K satisfy the condition (H), K = Q(z) and ǫ ∈ {−1, 1}. We need some definitions to define a c-map. For α = (α1 , . . . , αs )T ∈ D, and j ∈ [α,ǫ];(j) [α,ǫ];(j) [α,ǫ] = (g1 , . . . , gs ) Ind, we define linear fractional transformations Gj on D as follows: if αj 6= 0, for x := (x1 , . . . , xs )T ∈ D and i ∈ Ind, ord (α ) ordp (αj ) ǫp p j ǫp if i = j, − ωp [α,ǫ];(j) xj k αj gi (x) := k ǫp αi ǫp xi if i 6= j, − ωp xj αj where k = max{ordp (αj ) − ordp (αi ), 0}, and if αj = 0, for x = (x1 , . . . , xs )T ∈ D [α,ǫ] Gj (x) := (x). We define 0 φ : D → Ind as 0 φ(α) := 1 for α ∈ D. S = (sij ) ∈ GL(s, Zp ∩ Q) denotes the matrix ( δ(i+1)j for 1 ≤ i ≤ s − 1, 1 ≤ j ≤ s, sij := δ1j for i = s, 1 ≤ j ≤ s, where δii := 1 and δij := 0 for i 6= j for i, j ∈ Ind. Define the mapping D → Ind [ǫ] [α,ǫ] Φ0 (α) := (1, G1 , S), [1] for α ∈ D. We remark that for s = 1 the Φ0 continued fraction algorithm coincides with Schneider’s continued fraction algorithm. [α,ǫ,z] We assume that K 6= Q. We define linear fractional transformations Hj = [α,ǫ,z];(j) [α,ǫ,z];(j) , . . . , hs ) on D for j ∈ Ind as follows: (h1 [α,ǫ];(j) [α,ǫ];(j) (α) is written uniquely gi (α) = a0 + For α = (α1 , . . . , αs )T ∈ D gi s ′ a1 z + . . . + as z where ai ∈ Q for i ∈ Ind. Let a > 0 be the greatest common factor of the numerators of ai (1 ≤ i ≤ s) which is not divided by p. For 8 [α,ǫ,z];(j) x = (x1 , . . . , xs ) ∈ D we define hi [α,ǫ,z];(j) hi (x) for j ∈ Ind as [α,ǫ];(j) gi (x) := (x) a′ − a0 j a0 k + ′ . a′ a p Define [ǫ,z] Φ1 [α,ǫ,z] (α) := (1, H 1 , S), [ǫ,z],(n) for α ∈ D. We need some definitions to introduce the c-map Φ2 for n ∈ Z≥1 . For α ∈ K α is written uniquely α = a0 + a1 z + . . . + as z s where ai ∈ Q for 0 ≤ i ≤ s, and denomz (α) is defined by denomz (α) := min{|d| |d ∈ Z, d(a0 + a1 x + . . . + as xs ) ∈ Z[x]}, and for α = (α1 , . . . , αs )T ∈ K s we define denomz (α) := max{denomz (αi )|1 ≤ i ≤ s}. (1) We define v[ǫ,z] : D → Z for α = (α1 , . . . , αs )T ∈ D (1) [α,ǫ,z] v[ǫ,z] (α) := min{denomz (Hi (α))|1 ≤ i ≤ s}. (n) We define v[ǫ,z] : D → Z for n = 2, 3, . . . recursively as follows (n) [α,ǫ,z] v[ǫ,z] (α) := min{denomz (Hi [α,ǫ,z] (n−1) (α))v[ǫ,z] (Hi (α))|1 ≤ i ≤ s}. (n) Define φ[ǫ,z] : D → Ind for n ∈ Z≥1 (n+1) (n) [α,ǫ,z] φ[ǫ,z] (α) := min{i ∈ Ind|v[ǫ,z] (α) = denomz (Hi [ǫ,z],(n) for α = (α1 , . . . , αs )T ∈ D. Define Φ2 [ǫ,z],(n) Φ2 (n) (n) [α,ǫ,z] (α))v[ǫ,z] (Hi (α))}, for n ∈ Z≥1 (α) := (φ[ǫ,z] (α), H [α,ǫ,z] (n) φ[ǫ,z] (α) , id), where α ∈ D and id is the identity map. [ǫ,z] Let us define a c-map Φ3 . We need some definitions to define the c-map. We define ⌊α : m⌋p and hα : mip by ⌊α : m⌋p := Σn≤m,n∈Z cn pn , hα : mip := Σn>m,n∈Z cn pn . Mn (Q) denotes all n × n Q matrices for n ∈ Z≥1 . We define that M = (mij ) ∈ Mn (Q) is p-reduced, if M satisfies that for every integer i with 1 ≤ i ≤ n there exist uniquely an integer j(:= u(i)) with 0 ≤ j ≤ n such that (1) for every integer k with 1 ≤ k ≤ u(i) mik = 0, 9 (2) if u(i) 6= n, then miu(i)+1 ∈ {pl |l ∈ Z}, and for every integer k with i < k ≤ n mku(i)+1 = 0 (3) if i > 1, u(i) ≥ u(i − 1), and (4) if u(i) 6= n, then for every integer j with 1 ≤ j < i hmj u(i)+1 : ordp (miu(i)+1 ) − 1ip = 0. A matrix of Mn (Q) is converted to a p-reduced matrix by using the following row operations: (a) Switch two rows, (b) Multiply a row by an element of Vp ∩ Q, (c) Add a row to another one by multiplied by an element of Zp ∩ Q. For example, We introduce a following algorithm for M = (cij ) ∈ Mn (Q), which is given in Fig.4.1. The algorithm use elementally row operations (a), (b), (c). We call the algorithm p-reduced algorithm. The algorithm computes p-reduced matrices. Fig.4.1 Input M = (cij ) ∈ Mn (Q); Output M ′ = (cij ) ∈ Mn (Q); Let bi := (ci1 , . . . , cin ) for 1 ≤ i ≤ n; k1 = 1; k2 = 1; Do If there exists i with cik2 6= 0 for k1 ≤ i ≤ n, then Let m(k1 ≤ m ≤ n) be the least integer which satisfies |cmk2 |p = max{|cik2 |p |k1 ≤ i ≤ n, cik2 6= 0}; If m 6= k1 , swap bk1 and bm ; For i = k1 + 1, . . . , n cik2 bi = bi − bk ; ck1 k2 1 1 bk1 = bk ; |ck1 k2 |p ck1 k2 1 For i = 1, . . . , k1 − 1 hcik2 : ordp (ck1 k2 ) − 1ip bi = bi − bk1 ; ck1 k2 k1 = k1 + 1; k2 = k2 + 1; else k2 = k2 + 1; while(k2 ≤ n); Output M ′ = (cij ); When M ∈ Mn (Q) is converted by the p-reduced algorithm to M ′ ∈ Mn (Q), there exists N ∈ GL(n, Zp ∩ Q) associated with the algorithm such that M ′ = 10 N M and we denote by pr(M ) N . We can prove the following lemma in the usual way. Lemma 4.1. For M ∈ Mn (Q) and N1 , N2 ∈ GL(n, Zp ∩ Q) N1 M and N2 M are p-reduced, then N1 M = N2 M holds. We give an example. 10 3/2 Example. Let p = 2. M = is converted by the p-reduced algorithm −5 7 1 0 14/155 −3/155 to M ′ = and pr(M ) := . 0 1/2 1/31 2/31 We recall that z ∈ K satisfy the condition (H) and K = Q(z). For α = (α1 , . . . , αs )T ∈ D the Mα = (mij ) ∈ Ms×(s+1) (Q) is defined by α = Mα (z, z 2 , . . . , z s , 1)T and Mα′ ∈ Ms (Q) defines Mα′ := (mij )1≤i≤s,1≤j≤s . We define the map τz : D → Ms (Zp ∩ Q) τz (α) := pr(Mα′ ), for α ∈ D. We define the c-map [ǫ,z] Φ3 [α,ǫ,z] (α) := (1, H1 [α,ǫ,z] , τz (H1 )), for α ∈ D. 5 Extension of c-map Let Φ(α) = (φ(α), Fα , Aα ) for α = (α1 , . . . , αs )T ∈ D be a c-map and γ : D → pZp ∩Q be a map. We propose a p-adic continued fraction via Φ and γ as follows: Let α = (α1 , . . . , αs )T ∈ D. We set α(0) := α. α(1) , α(2) , . . . are defined inductively as follows: We suppose that α(n) for n ∈ Z≥0 is defined. We set α(n+1) := Aα(n) Fα(n) (α(n) ) + γ(α(n) ). We say that α has the affine c-map continued fraction expansion {Φ′ (α(0) ), Φ′ (α(1) ), . . .}, where Φ′ (β) = (φ(β), Fβ , Aβ , γ(β)) for β ∈ D. We can easily see that Theorem 3.6 holds for an affine c-map. [·,ǫ,z] [·,ǫ,z] [·,ǫ,z] For the c-map Φ3 := (1, H1 , τz (H1 )) we will define an affine c-map. [α,ǫ,z] [α,ǫ,z] Let α = (α1 , . . . , αs )T ∈ D and let β = (β1 , . . . , βs )T := τz (H1 (α))(H1 (α)). Then, there exists a unique M = (mij ) ∈ Ms×(s+1) (Q) such that β = M (z, z 2 , . . . , z s , 1)T . We define γ ′ (α) ∈ E as γ ′ (α) := (−hm1s+1 ip , . . . , −hmss+1 ip )T . [ǫ,z] [·,ǫ,z] We define the affine c-map Φ3A := (1, H1 11 [·,ǫ,z] , τz (H1 ), γ ′ (·)). 6 Rational and Quadratic Cases In this section, we consider a periodicity of our algorithms in specific cases. It is well-known that every rational number has a finite continued fraction expansion. [−1] We will see that a similar result holds for the Φ0 continued fraction algorithm. Proposition 6.1. Let K be Q. Then, for every rational number α α has a [−1] finite Φ0 continued fraction expansion. Proof. We define the height of rational number by the summation of the absolute value of its numerator and the absolute one of its denominator, which is denoted by height(). Let α be a rational number. By Lemma 3.1 we assume α ∈ pZp . [−1] If α = 0, then α has apparently a finite Φ0 continued fraction expansion. k Let α = mm1 p2 6= 0 where m1 , m2 are relatively prime integers with ordp (m1 ) = ordp (m2 ) = 0 and k ∈ Z≥1 . Then, we see that the next remainder α1 is −m2 −jm1 2 where j = ωp (− m m1 m1 ). We have height(α) = |m1 pk | + |m2 | ≥ p|m1 | + |m2 | ≥ | − m2 − jm1 | + |m1 | = height(α1 ), in which the equality holds for α > 0, j = p − 1 and k = 1. If the sign of α is positive, the sign of the next remainder is minus. Therefore, if α has an [−1] infinite Φ0 continued fraction expansion, we see height(α) > height(α2 ) > height(α4 ) > . . ., which is a contradiction. [1] For the Φ0 continued fraction algorithm on Q, which coincides with Schneider’s continued fraction algorithm Bundschuh [3] gave results, in which every rational number has an infinite periodic expansion or a finite expansion. Next, we consider quadratic cases. Weger [13] showed that some quadratic elements have a non periodic Schneider’s continued fraction expansion. For the [−1] Φ0 continued fraction we can not expect that every quadratic element have a periodic continued fraction expansion from the numerical experiments (see Section 8), but we can not give the proof. [ǫ,z] For the Φ1 continued fraction we give following Theorem. Lemma 6.2. Let K be a quadratic field over Q. Let z ∈ K satisfy the condition [ǫ,z] (H) and K = Q(z). For q ∈ Q with ordp (qz) > 0, qz has a periodic Φ1 continued fraction expansion. Proof. Let x2 +ux+vpk be a minimal polynomial of z, where u, v ∈ Q, ordp (u) = v1 m1 , where v1 , v2 are relatively ordp (v) = 0 and k > 0. Let v = and q = m 2 v2 prime integers with v2 > 0 and m1 , m2 are relatively prime integers with m2 > 0. First, we suppose that ordp (q) = 0. Since ordp (z) = k, we have G[qz,ǫ] (qz) = ǫ m2 p k m2 uv2 m2 v2 z −ǫ − j, − j = −ǫ m1 z m1 v1 m1 v1 12 k where j = ωp (ǫ mm21pz ). Let a′ ∈ Z>0 be the numerator of −ǫ have H [qz,ǫ,z] (qz) = −ǫ m2 v2 . Then, we m1 v1 m2 uv2 j m2 v2 z + −ǫ − . a′ m1 v1 a′ m1 v1 a′ p m2 zv2 2 uv2 ∈ pZp , we have −ǫ am′ m − aj′ ∈ pZp ∩ Q, which 1 v1 a′ m1 v1 m2 uv2 j implies ⌊−ǫ ′ − ⌋p = 0. Let b be the greatest common factor of m2 and a m1 v1 a′ v1 and b′ be the greatest common factor of m1 and v2 . Then, we have z . H [qz,ǫ,z] (qz) = −ǫ ′ (m1 /b )(v1 /b) Since H [α,ǫ,z] (qz), −ǫ Therefore, the nth remainder of qz for n ≥ 2 is δvnnz , where δn ∈ {1, −1} and vn is a factor of v1 . Since the set { δvnnz |n ∈ Z≥2 } is a finite set, qz has a periodic [ǫ,z] Φ1 continued fraction expansion. Next, we suppose that ordp (q) = g > 0. Let m1 = m′1 pg , where m′1 ∈ Z. Then, we have G[qz,ǫ] (qz) = ǫ k+g p where j = ω(ǫ m2nz m2 pk+g m2 uv2 m2 v2 z −ǫ ′ − j, − j = −ǫ ′ m1 z m1 v1 m1 v1 ). Similarly, we have H [qz,ǫ,z] (qz) = −ǫ z , (m′1 /b′ )(v1 /b) where b is the greatest common factor of m2 and v1 and b′ is the greatest common factor of m′1 and v2 . Since ordp ( (m′ /b′1)(v1 /b) ) = 0 holds, by the previous 1 argument it follows Lemma. Finally, we suppose that ordp (q) = −g < 0. Let m2 = m′2 pg , where m′2 ∈ Z. Then, we have G[qz,ǫ] (qz) = ǫ m2 pk−g m′ uv2 m′ v2 z −ǫ 2 − j, − j = −ǫ 2 m1 z m1 v1 m1 v1 k−g 2p where j = ωp (ǫ mm ). Similarly, we have 1z H [qz,ǫ,z] (qz) = −ǫ z , (m1 /b′ )(v1 /b) where b is the greatest common factor of m′2 and v1 and b′ is the greatest common factor of m1 and v2 . Since ordp ( (m1 /b′1)(v1 /b) ) = 0 holds, by the previous argument it follows Lemma. Theorem 6.3. Let K be a quadratic field over Q. Let z ∈ K satisfy the condition (H) and K = Q(z). Then, for every rational number α α has a finite [ǫ,z] Φ1 continued fraction expansion. For every α ∈ K with α ∈ / Q, α has a [ǫ,z] periodic Φ1 continued fraction expansion. 13 Proof. Let x2 +ux+vpk be a minimal polynomial of z, where u, v ∈ Q, ordp (u) = v1 ordp (v) = 0 and k > 0. Let v = , where v1 , v2 are relatively prime integer. v2 Let α 6= 0 be a rational number. By Lemma 3.1 we assume α ∈ pZp . Let k1 1p 6= 0 such that m1 , m2 are relatively prime integers with ordp (m1 ) = α = mm 2 2 −jm1 ordp (m2 ) = 0 and k1 ∈ Z≥1 . Then, we see that G[α,ǫ] (α) = −ǫmm where 1 −ǫm2 −jm1 ǫm2 [α,ǫ,z] j = ωp (− m1 ). Since ∈ pZp ∩ Q, H (α) = 0. Therefore, α has a m1 [ǫ,z] finite Φ1 continued fraction expansion. Let α ∈ K be not rational. Then, we have H [α,ǫ,z] (α) = qz or H [α,ǫ,z] (α) = z + q′ , m where m 6= 0 ∈ Z and q, q ′ ∈ Q with |q ′ |p ≥ 1. If H [α,ǫ,z] (α) = qz holds, by [ǫ,z] Lemma 6.2 α has a periodic Φ1 continued fraction expansion. We suppose z z + q ′ . Let k2 := ordp ( + q ′ ) > 0. First, we suppose that H [α,ǫ,z] (α) = m m z ordp (q ′ ) = 0. Then, since we see ordp ( ) = 0, we have ordp (m) = k. Let m m = pk m′ , where m′ ∈ Z and ordp (m′ ) = 0. Then, we get ǫm′ pk ǫ = (−z + m′ pk q ′ − u). z ′ k ′ 2 ′ pk q ′ + pk v (m p q ) − um ′ +q m We consider the coefficient of the term z of the above formula. We have |(m′ pk q ′ )2 − um′ pk q ′ + v|p = |(m′ pk q ′ )2 − u(m′ pk q ′ + z) + uz + pk v|p = |(m′ pk q ′ )2 − u(m′ pk q ′ + z) − z 2 |p ′ k ′ 2 z 2 ′ k ′ = (m p q ) − z − m p u q + ′ k mp p ′ k 2 ′ z z z = (m p ) q + ′ k q ′ − ′ k − m′ pk u q ′ + ′ k mp mp mp p ′ k z = m p u q ′ + ′ k = p−(k+k2 ) . mp p Therefore, we have −ǫm′ pk+k2 z (m′ pk q ′ )2 − um′ pk q ′ + v < 1. p ǫpk2 ǫpk2 < 1 we have By (1) and z − ωp z ′ ′ + q + q m m p ′ k+k2 k2 ǫm p ǫp ′ k < 1, z (m′ pk q ′ )2 − um′ pk q ′ + v (m p j − u) − ωp + q′ m p 14 (1) which yields z ′ H [ m +q ,ǫ,z] z −ǫm′ pk+k2 + q′ = z. ′ k ′ m (m p q )2 − um′ pk q ′ + pk v [ǫ,z] By Lemma 6.2 α has a periodic Φ1 continued fraction expansion. Next, we z suppose ordp (q ′ ) = −l < 0. Then, since ordp ( m ) = −l, we have ordp (m) = k+l. k+l ′ ′ Let m = p m , where m ∈ Z and ordp (m′ ) = 0. Then, by the previous argument there exists a integer j ′ 6= 0 such that z z [ ′z k +q′ pl ,ǫ,z] ′ l m p H + q p = ′. m′ p k j z ′ l z ′ z z + q ′ ) = H [p ( m +q ),ǫ,z] (pl ( m + q ′ )) holds, by Lemma 6.2 α Since H [ m +q ,ǫ,z] ( m [ǫ,z] has a periodic Φ1 continued fraction expansion. By the proof of Theorem 6.3 we have Corollary 6.4. Let K be a quadratic field over Q. Let z ∈ K satisfy the condition (H) and K = Q(z). Let x2 + ux + vpk be a minimal polynomial of z, v1 where u, v ∈ Q , ordp (u) = ordp (v) = 0 and k > 0. Let v = , where v1 , v2 v2 are relatively prime integer with v2 > 0. For α ∈ K with α ∈ / Q, α has a purely [ǫ,z] periodic Φ1 continued fraction expansion if and only if α∈{ δ |m ∈ Z>0 is a factor of v1 and δ ∈ {−1, 1}}. m [ǫ,z] [ǫ,z],(n) Remark. We remark that if K is a quadratic field, then c-map Φ1 = Φ2 hold for n ∈ Z≥1 . Therefore, these continued fraction algorithms coincide with each other. In the similar manner, we have Theorem 6.5. Let K ⊂ Q be a quadratic field over Q. Then, for every rational [ǫ,z] number α α has a finite Φ3 continued fraction expansion. For every α ∈ K [ǫ,z] with α ∈ / Q, α has a periodic Φ3 continued fraction expansion. 7 Multidimensional Cases We can expect that a modification of Lagrange’ Theorem holds for some our algorithm from numerical experiments (see Section 8). But we can not give the proofs. Dubois, Paysant-Le [5] showed that for every real cubic number field there is a pair of numbers which has a periodic Jacobi-Perron expansion. In this [ǫ,z] section, we show that a modified result holds for the Φ1 continued fraction [ǫ,z] algorithm and there exist many elements which have periodic Φ3 continued fraction expansions for any finite extension of Q included in Qp . 15 Let K be a cubic field over Q and K ⊂ Qp . Let z ∈ K satisfy the condition (H) and K = Q(z). Since mz satisfies above conditions for an arbitrary integer m which is relatively prime to p, we can choose z which is integral over Q. Let z be integral over Q and x3 + a1 x2 + a2 x + a3 pk (2) be the minimal polynomial of z, where ai ∈ Z for 1 ≤ i ≤ 3 , ordp (a2 ) = ordp (a3 ) = 0 and k ∈ Z>0 . [ǫ,z] Theorem 7.1. Let α := (z, z 2 )T . Then, α has a periodic Φ1 fraction expansion. Proof. We have [α,ǫ];(1) g1 (α) = ǫpk − ωp z ǫpk z −ǫz 2 −ǫa1 z −ǫa2 = + + − ωp a3 a3 a3 ǫpk z . Since the greatest common factor of −ǫ and −ǫa1 is 1, we have k −ǫz 2 −ǫa1 z −ǫa2 ǫp [α,ǫ,z];(1) h1 . + + − ωp (α) = a3 a3 a3 z p −ǫz 2 −ǫa1 z [α,ǫ,z];(1) < 1 imply (α)|p < 1 and + |h1 a3 a3 p k ǫp −ǫa2 = 0. − ωp a3 z p Therefore, we have [α,ǫ,z];(1) h1 (α) = −ǫz 2 −ǫa1 z + . a3 a3 Thus, we have [α,ǫ,z] α1 = SH1 (α) = T −ǫz 2 −ǫa1 z ǫz, + . a3 a3 Similarly, we have T −ǫz −z 2 −a1 z α2 = = , + , a3 a3 a3 T [α ,ǫ,z] α3 = SH1 2 (α2 ) = z, z 2 + a1 z , T −ǫz 2 −ǫa1 z [α3 ,ǫ,z] α4 = SH1 + (α3 ) = ǫz, . a3 a3 [α ,ǫ,z] SH1 1 (α1 ) [ǫ,z] Since α1 = α4 , α has a periodic Φ1 continued fraction expansion. 16 continued Let K ⊂ Qp be a finite extension of Q and K 6= Q. Let z ∈ K satisfy the condition (H) and K = Q(z). Let xs+1 + a1 xs + . . . + as+1 pk (3) be the minimal polynomial of z, where ai ∈ Q ∩ Zp for 1 ≤ i ≤ s + 1 and ordp (as ) = ordp (as+1 ) = 0. We set a0 := 1. P j Theorem 7.2. Let u1 := z and ui := i≤j≤s aij z for 2 ≤ i ≤ s, where aij ∈ Q ∩ Zp for 1 ≤ i ≤ s, i ≤ j ≤ s and ordp (aii ) ∈ Vp for 1 ≤ i ≤ s. Then, [ǫ,z] α := (u1 , . . . , us )T has a periodic Φ3 continued fraction expansion. Proof. We have ǫpk z k X −ǫas−i z i ǫp = . − ωp as+1 z [α,ǫ];(1) (α) g1 ǫpk − ωp = z 0≤i≤s s−i The greatest common factor of the numerators of −ǫa as+1 for 1 ≤ i ≤ s is 1. In the same way as the proof of Theorem 7.1 we have [α,ǫ,z];(1) h1 (α) = X −ǫas−i z i . as+1 1≤i≤s For 2 ≤ i ≤ s we have ǫu ǫui i − ωp z z ǫu X i . = −ǫaij z j−1 − ωp z [α,ǫ];(i) g1 (α) = i≤j≤s Let a′i > 0 ∈ Z be the greatest common factor of the numerators of −ǫaij for i ≤ j ≤ s. Similarly, for 2 ≤ i ≤ s we have [α,ǫ,z];(i) h1 X −ǫaij z j−1 . a′i (α) = i≤j≤s Then, Mα′ ∈ Ms (Q) is given by −ǫas−1 as+1 −ǫa22 a′ 2 ′ Mα = 0 .. . 0 ... ... ... ... −ǫa33 a′3 ... 0 .. . ... 17 ... −ǫa2s a′2 −ǫa3s a′3 −ǫass a′s −ǫa0 as+1 0 0 . .. . 0 ii We see that every element of Mα′ is in Zp . ordp ( −ǫa a′ ) = 0 for 2 ≤ i ≤ s i ′ 0 and ordp ( −ǫa as+1 ) = 0 implies that Mα is converted to the unit matrix E by the ′ p-reduced algorithm. Since pr(Mα )Mα′ = E holds, we have [α,ǫ,z];(1) α2 = pr(Mα′ )(h1 [α,ǫ,z];(s) (α), . . . , h1 (α))T = (z, z 2 , . . . , z s )T . [ǫ,z] Similarly, we have α3 = (z, z 2 , . . . , z s )T . Therefore, α has a periodic Φ3 continued fraction expansion. 8 Numerical Experiments In this section, we give some numerical results about our algorithms. For the calculation of the tables, we used some computers equipped with GiNaC [4] on GNU C++. In Table 1 for the prime numbers p with 2 ≤ p ≤ 100 and {z ∈ Qp |x2 +ax+bp is the minimal polynomial of z, a, b ∈ Z, 0 < a ≦ 10, −10 ≤ b ≤ 10, ordp (a) = 0} and a set of 100 elements in D we observe periodicity by the [α,ǫ,z] Φ0 continued fraction algorithm. We generate the set of elements denoted by T estz in D by using the algorithm given by Saito, Yamaguchi [10] as follows: Let 0.d1 d2 · · · be the binary expansion of the real positive root of x2 + x − 1, where {d1 , d2 · · · } is generated by the algorithm [10] in which they showed that the sequence has good properties as pseudorandom numbers. We set ei = 8 X 2k−1 d8i+k for i ∈ Z≥0 , k=1 ti = (−1)e6i+2 e6i+4 e6i+1 + (−1)e6i+5 z for i ∈ Z≥0 . e6i + 1 e6i+3 + 1 Let m be the least integer such that ♯({ti |i ≤ m}\{ti |ti ∈ Q, i ≤ m}) = 100. We define T estz ⊂ D as T estz := {ti |i ≤ m}\{ti |ti ∈ Q, i ≤ m}. Let us define a notation. For α ∈ K, there exist ai ∈ Q for i = 0, 1, . . . , s such that α = a0 + a1 z + . . . + as z s . Heightz (α) is defined as the maximum of the heights of ai for i = 0, 1, . . . , s. For α = (α1 , . . . , αs ) ∈ K s , Heightz (α) is defined as the maximum of Heightz (αi ) for i = 1, . . . , s. In Table 1 we do not iterate the algorithm when the height of a remainder exceeds 10300 . 18 Table 1 prime number 2 3 5 7 11 13 17 19 23 29 31 37 41 1∗ 2∗ 3∗ 4∗ 0 0 0 0 0 0 0 0 0 0 0 0 0 7800 11700 14400 16600 19000 19200 9600 19800 20000 20000 20000 20000 20000 0 0 1 2 0 0 0 0 0 0 0 0 0 7800 11700 14399 16598 19000 19200 19600 19800 20000 20000 20000 20000 20000 prime number 43 47 53 59 61 67 71 73 79 83 89 97 1∗ 2∗ 3∗ 4∗ 0 0 0 0 0 0 0 0 0 0 0 0 20000 20000 20000 20000 20000 20000 20000 20000 20000 20000 20000 20000 0 0 0 0 0 0 0 0 0 0 0 0 20000 20000 20000 20000 20000 20000 20000 20000 20000 20000 20000 20000 [1,z] [−1,z] 1∗ (resp.,3∗ )is the number of elements which have a periodic Φ0 (resp.,Φ0 ) continued fraction expansion and 2∗ (resp.,4∗ ) is the number of elements such [1];z [−1,z] that their heights of some remainders with respect to Φ0 (resp.,Φ0 ) exceed 300 10 and their successive iterations terminate. In Table 2 for the prime numbers p with 2 ≤ p ≤ 100 and {z ∈ Qp |x3 + ax + bp is the minimal polynomial of z, a, b ∈ Z, 0 < a ≦ 10, −10 ≤ b ≤ (2) 10, ordp (a) = 0} and the set of 100 elements in D denoted by T estz (de[ǫ,z] scribed later) we observe periodicity by the Φ1 continued fraction algorithm. (1) (1) (2) (2) Let 0.d1 d2 · · · (resp., 0.d1 d2 · · · ) be the binary expansion of the real posi(j) (j) 2 tive root of x + 2x − 1(resp.,x2 + 2x − 2), where {d1 , d2 · · · } for j = 1, 2 are generated by the algorithm [10]. We set that for j = 1, 2 (j) ei = 8 X (j) 2k−1 d8i+k for i ∈ Z≥0 , k=1 (j) ti (j) e9i+1 (j) e9i+2 = (−1) (j) e9i + 1 (j) e9i+4 (j) e9i+5 + (−1) (j) (j) e9i+8 e9i+3 + 1 z + (−1) (j) e9i+7 (j) e9i+6 + 1 z 2 for i ∈ Z≥0 . Let m be the least integer such that (1) (1) (2) (2) (1) (2) ♯({(ti , ti )|i ≤ m}\{(ti , ti )|1, ti , ti are linearly independent over Q, i ≤ m}) = 100. (2) We define T estz ⊂ D as (1) (2) (1) (2) (1) (2) T est(2) z := {(ti , ti )|i ≤ m}\{(ti , ti )|1, ti , ti 19 are linearly independent over Q, i ≤ m}. Table 2 prime number 2 3 5 7 11 13 17 19 23 29 31 37 41 1∗ 2∗ 3∗ 4∗ 7514 11052 12534 14872 16424 17388 18063 17956 18099 17493 18627 19331 18395 886 1548 2266 2328 2776 1812 1337 1644 1701 2307 1173 469 1405 7533 11706 12576 14843 16441 17383 18083 17934 18115 17496 18636 19325 18395 867 1524 2224 2357 2759 1817 1317 1666 1685 2304 1164 475 1405 prime number 43 47 53 59 61 67 71 73 79 83 89 97 1∗ 2∗ 3∗ 4∗ 19183 18972 18593 19164 19776 19577 19598 18965 19203 19774 19573 19335 617 1028 1207 636 224 223 202 835 797 26 227 665 19179 18982 18599 19160 19766 19569 19601 18961 19203 19780 19576 19336 621 1018 1201 640 224 231 199 839 797 20 224 664 [1,z] [−1,z] 1∗ (resp.,3∗ )is the number of elements which have the periodic Φ1 (resp.,Φ1 ) continued fraction expansion and 2∗ (resp.,4∗ ) is the number of elements such [1,z] [−1,z] that their heights of some remainders with respect to Φ1 (resp.,Φ1 ) exceed 10300 and their successive iterations terminate. In Table 3 for the prime numbers p with 2 ≤ p ≤ 100 and {z ∈ Qp |x3 + ax + bp is the minimal polynomial of z, a, b ∈ Z, 0 < a ≦ 10, −10 ≤ b ≤ 10, ordp (a) = 0} and 100 elements in D given in the same way as Table 2 [ǫ,z],(1) we observe periodicity by the Φ2 continued fraction algorithm. Table 3 prime number 2 3 5 7 11 13 17 19 23 29 31 37 41 1∗ 2∗ 3∗ 4∗ 7986 12227 14214 16874 18909 18977 19226 19486 19726 19739 19791 19721 19766 414 373 586 326 291 223 174 114 74 61 9 79 34 7976 12228 14203 1616870 18909 18974 19228 19497 19728 19745 19790 19717 19767 424 372 597 330 291 226 172 103 72 55 10 83 33 prime number 43 47 53 59 61 67 71 73 79 83 89 97 1∗ 2∗ 3∗ 4∗ 19738 19992 19777 19773 19973 19775 19790 19787 19994 19799 19781 19989 62 8 23 27 27 25 10 13 6 1 19 11 19737 19994 19777 19774 19971 19773 19790 19786 19994 19799 19781 19990 63 6 23 26 29 27 10 14 6 1 19 10 [1,z],(1) 1∗ (resp.,3∗ )is the number of elements which have the periodic Φ2 20 [−1,z],(1) (resp.,Φ2 ) continued fraction expansion and 2∗ (resp.,4∗ ) is the number of elements such [1,z],(1) [−1,z],(1) that their heights of some remainders with respect to Φ2 (resp.,Φ2 ) exceed 10300 and their successive iterations terminate. In Table 4 for the prime numbers p with 2 ≤ p ≤ 100 and {z ∈ Qp |x3 + ax + bp is the minimal polynomial of z, a, b ∈ Z, 0 < a ≦ 10, −10 ≤ b ≤ 10, ordp (a) = 0} and 100 elements in D given in the same way as Table 2 [ǫ,z],(2) we observe periodicity by the Φ2 continued fraction algorithm. Table 4 prime number 2 3 5 7 11 13 17 19 23 29 31 37 41 1∗ 2∗ 3∗ 4∗ 8393 12584 14792 17196 19200 19200 19397 19600 19800 19800 19800 19800 19800 7 16 8 4 0 0 3 0 0 0 0 0 0 8390 12590 14790 17193 19200 19200 19396 19600 19800 19800 19800 19800 19800 10 10 10 7 0 0 4 0 0 0 0 0 0 prime number 43 47 53 59 61 67 71 73 79 83 89 97 1∗ 2∗ 3∗ 4∗ 19800 20000 19800 19800 20000 19800 19800 19800 20000 19800 19800 20000 0 0 0 0 0 0 0 0 0 0 0 0 19800 20000 19800 19800 20000 19800 19800 19800 20000 19800 19800 20000 0 0 0 0 0 0 0 0 0 0 0 0 [1,z],(2) [−1,z],(2) 1∗ (resp.,3∗ )is the number of elements which have the periodic Φ2 (resp.,Φ2 continued fraction expansion and 2∗ (resp.,4∗ ) is the number of elements such [1,z],(2) [−1,z],(2) that their heights of some remainders with respect to Φ2 (resp.,Φ2 ) 300 exceed 10 and their successive iterations terminate. In Table 5 for the prime numbers p with 2 ≤ p ≤ 100 and {z ∈ Qp |x3 + ax + bp is the minimal polynomial of z, a, b ∈ Z, 0 < a ≦ 10, −10 ≤ b ≤ 10, ordp (a) = 0} and 100 elements in D given in the same way as Table 2 [ǫ,z] we observe periodicity by the Φ3 continued fraction algorithm. 21 ) Table 5 prime number 2 3 5 7 11 13 17 19 23 29 31 37 41 1∗ 2∗ 3∗ 4∗ 8395 12600 14800 17200 19200 19200 19400 19600 19800 19800 19800 19800 19800 5 0 0 0 0 0 0 0 0 0 0 0 0 8394 12600 14800 17200 19200 19200 19400 19600 19800 19800 19800 19800 19800 6 0 0 0 0 0 0 0 0 0 0 0 0 prime number 43 47 53 59 61 67 71 73 79 83 89 97 1∗ 2∗ 3∗ 4∗ 19800 20000 19800 19800 20000 19800 19800 19800 20000 19800 19800 20000 0 0 0 0 0 0 0 0 0 0 0 0 19800 20000 19800 19800 20000 19800 19800 19800 20000 19800 19800 20000 0 0 0 0 0 0 0 0 0 0 0 0 [1,z] [−1,z] 1∗ (resp.,3∗ )is the number of elements which have the periodic Φ3 (resp.,Φ3 ) continued fraction expansion and 2∗ (resp.,4∗ ) is the number of elements such [1,z] [−1,z] that their heights of some remainders with respect to Φ3 (resp.,Φ3 ) exceed 10300 and their successive iterations terminate. In Table 6 for the prime numbers p with 2 ≤ p ≤ 100 and {z ∈ Qp |x4 + ax + bp is the minimal polynomial of z, a, b ∈ Z, 0 < a ≦ 10, −10 ≤ b ≤ 10, ordp (a) = 0} and 100 elements in D given in the same way as Table 2 [ǫ,z] we observe periodicity by the Φ3 continued fraction algorithm. Table 6 prime number 2 3 5 7 11 13 17 19 23 29 31 37 41 1∗ 2∗ 3∗ 4∗ 8100 12900 15200 17400 19600 19800 19700 19900 19900 19800 19800 19800 20000 0 0 0 0 0 0 0 0 0 0 0 0 0 8100 12899 15200 17400 19600 19800 19700 19900 19900 19800 19800 19800 20000 0 1 0 0 0 0 0 0 0 0 0 0 0 prime number 43 47 53 59 61 67 71 73 79 83 89 97 1∗ 2∗ 3∗ 4∗ 20000 20000 20000 19800 19800 19800 19900 19900 20000 20000 20000 20000 0 0 0 0 0 0 0 0 0 0 0 0 20000 20000 20000 19800 19800 19800 19900 19900 20000 20000 20000 20000 0 0 0 0 0 0 0 0 0 0 0 0 [1,z] 1∗ (resp.,3∗ )is the number of elements which have a periodic Φ3 22 [−1,z] (resp.,Φ3 ) continued fraction expansion and 2∗ (resp.,4∗ ) is the number of elements such [1,z] [−1,z] that their heights of some remainders with respect to Φ3 (resp.,Φ3 ) exceed 10300 and their successive iterations terminate. In Table 7 for the prime numbers p with 2 ≤ p ≤ 100 and {z ∈ Qp |x5 + ax + bp is the minimal polynomial of z, a, b ∈ Z, 0 < a ≦ 10, −10 ≤ b ≤ 10, ordp (a) = 0} and 100 elements in D given in the same way as Table 2 [ǫ,z] we observe periodicity by the Φ3 continued fraction algorithm. Table 7 prime number 2 3 5 7 11 13 17 19 23 29 31 37 41 1∗ 2∗ 3∗ 4∗ 8794 13200 15200 17600 19600 19800 19800 19800 19800 19800 20000 20000 19800 6 0 0 0 0 0 0 0 0 0 0 0 0 8796 13200 15200 17600 19600 19800 19800 19800 19800 19800 20000 20000 19800 4 0 0 0 0 0 0 0 0 0 0 0 0 prime number 43 47 53 59 61 67 71 73 79 83 89 97 1∗ 2∗ 3∗ 4∗ 19800 20000 20000 20000 20000 20000 20000 20000 20000 19800 19800 20000 0 0 0 0 0 0 0 0 0 0 0 0 19800 20000 20000 20000 20000 20000 20000 19999 20000 19800 19800 20000 0 0 0 0 0 0 0 1 0 0 0 0 [1,z] [−1,z] 1∗ (resp.,3∗ )is the number of elements which have the periodic Φ3 (resp.,Φ3 ) continued fraction expansion and 2∗ (resp.,4∗ ) is the number of elements such [1,z] [−1,z] that their heights of some remainders with respect to Φ3 (resp.,Φ3 ) exceed 300 10 and their successive iterations terminate. In Table 8 for the prime numbers p with 2 ≤ p ≤ 100 and {z ∈ Qp |x6 + ax + bp is the minimal polynomial of z, a, b ∈ Z, 0 < a ≦ 10, −10 ≤ b ≤ 10, ordp (a) = 0} and 100 elements in D given in the same way as Table 2 [ǫ,z] we observe periodicity by the Φ3 continued fraction algorithm. 23 Table 8 prime number 2 3 5 7 11 13 17 19 23 29 31 37 41 1∗ 2∗ 3∗ 4∗ 9000 13400 15499 17600 19600 19800 19900 19900 19900 19900 19900 19900 19900 0 0 1 0 0 0 0 0 0 0 0 0 0 9000 13400 15498 17600 19600 19800 19900 19900 19900 19900 19900 19900 19900 0 0 2 0 0 0 0 0 0 0 0 0 0 prime number 43 47 53 59 61 67 71 73 79 83 89 97 1∗ 2∗ 3∗ 4∗ 20000 20000 20000 20000 20000 20000 20000 20000 19900 19900 20000 20000 0 0 0 0 0 0 0 0 0 0 0 0 20000 20000 20000 20000 20000 20000 20000 20000 19900 19900 20000 20000 0 0 0 0 0 0 0 0 0 0 0 0 [1,z] [−1,z] 1∗ (resp.,3∗ )is the number of elements which have the periodic Φ3 (resp.,Φ3 ) continued fraction expansion and 2∗ (resp.,4∗ ) is the number of elements such [1,z] [−1,z] that their heights of some remainders with respect to Φ3 (resp.,Φ3 ) exceed 10300 and their successive iterations terminate. In Table 9 for the prime numbers p with 2 ≤ p ≤ 100 and {z ∈ Qp |x3 + ax + bp is the minimal polynomial of z, a, b ∈ Z, 0 < a ≦ 10, −10 ≤ b ≤ 10, ordp (a) = 0} and 100 elements in D given in the same way as Table 2 [ǫ,z] we observe periodicity by the Φ3A continued fraction algorithm. Table 9 prime number 2 3 5 7 11 13 17 19 23 29 31 37 41 1∗ 2∗ 3∗ 4∗ 8400 12600 14800 17200 19200 19200 19400 19600 19800 19800 19800 19800 19800 0 0 0 0 0 0 0 0 0 0 0 0 0 8400 12600 14800 17200 19200 19200 19400 19600 19800 19800 19800 19800 19800 0 0 0 0 0 0 0 0 0 0 0 0 0 prime number 43 47 53 59 61 67 71 73 79 83 89 97 1∗ 2∗ 3∗ 4∗ 19800 20000 19800 19800 20000 19800 19800 19800 20000 19800 19800 20000 0 0 0 0 0 0 0 0 0 0 0 0 19800 20000 19800 19800 20000 19800 19800 19800 20000 19800 19800 20000 0 0 0 0 0 0 0 0 0 0 0 0 [1,z] [−1,z] 1∗ (resp.,3∗ )is the number of elements which have the periodic Φ3A (resp.,Φ3A 24 ) continued fraction expansion and 2∗ (resp.,4∗ ) is the number of elements such [−1,z] [1,z] that their heights of some remainders with respect to Φ3A (resp.,Φ3A ) exceed 10300 and their successive iterations terminate. In Table 10 for the prime numbers p with 2 ≤ p ≤ 100 and {z ∈ Qp |x4 + ax + bp is the minimal polynomial of z, a, b ∈ Z, 0 < a ≦ 10, −10 ≤ b ≤ 10, ordp (a) = 0} and 100 elements in D given in the same way as Table 2 [ǫ,z] we observe periodicity by the Φ3A continued fraction algorithm. Table 10 prime number 2 3 5 7 11 13 17 19 23 29 31 37 41 1∗ 2∗ 3∗ 4∗ 8100 12900 15200 17400 19600 19800 19700 19900 19900 19800 19800 19800 20000 0 0 0 0 0 0 0 0 0 0 0 0 0 8100 12900 15200 17400 19600 19800 19700 19900 19900 19800 19800 19800 20000 0 0 0 0 0 0 0 0 0 0 0 0 0 prime number 43 47 53 59 61 67 71 73 79 83 89 97 1∗ 2∗ 3∗ 4∗ 20000 20000 20000 19800 19800 19800 19900 19900 20000 20000 20000 20000 0 0 0 0 0 0 0 0 0 0 0 0 20000 20000 20000 19800 19800 19800 19900 19900 20000 20000 20000 20000 0 0 0 0 0 0 0 0 0 0 0 0 [1,z] [−1,z] 1∗ (resp.,3∗ )is the number of elements which have the periodic Φ3A (resp.,Φ3A ) continued fraction expansion and 2∗ (resp.,4∗ ) is the number of elements such [1,z] [−1,z] that their heights of some remainders with respect to Φ3A (resp.,Φ3A ) exceed 300 10 and their successive iterations terminate. In Table 11 for the prime numbers p with 2 ≤ p ≤ 100 and {z ∈ Qp |x5 + ax + bp is the minimal polynomial of z, a, b ∈ Z, 0 < a ≦ 10, −10 ≤ b ≤ 10, ordp (a) = 0} and 100 elements in D given in the same way as Table 2 [ǫ,z] we observe periodicity by the Φ3 continued fraction algorithm. 25 Table 11 prime number 2 3 5 7 11 13 17 19 23 29 31 37 41 1∗ 2∗ 3∗ 4∗ 8800 13200 15200 17600 19600 19800 19800 19800 19800 19800 20000 20000 19800 0 0 0 0 0 0 0 0 0 0 0 0 0 8800 13200 15200 17600 19600 19800 19800 19800 19800 19800 20000 20000 19800 0 0 0 0 0 0 0 0 0 0 0 0 0 prime number 43 47 53 59 61 67 71 73 79 83 89 97 1∗ 2∗ 3∗ 4∗ 19800 20000 20000 20000 20000 20000 20000 20000 20000 19800 19800 20000 0 0 0 0 0 0 0 0 0 0 0 0 19800 20000 20000 20000 20000 20000 20000 20000 20000 19800 19800 20000 0 0 0 0 0 0 0 0 0 0 0 0 [1,z] [−1,z] 1∗ (resp.,3∗ )is the number of elements which have the periodic Φ3A (resp.,Φ3A ) continued fraction expansion and 2∗ (resp.,4∗ ) is the number of elements such [−1,z] [1,z] that their heights of some remainders with respect to Φ3A (resp.,Φ3A ) exceed 10300 and their successive iterations terminate. In Table 12 for the prime numbers p with 2 ≤ p ≤ 100 and {z ∈ Qp |x6 + ax + bp is the minimal polynomial of z, a, b ∈ Z, 0 < a ≦ 10, −10 ≤ b ≤ 10, ordp (a) = 0} and 100 elements in D given in the same way as Table 2 [ǫ,z] we observe periodicity by the Φ3A continued fraction algorithm. Table 12 prime number 2 3 5 7 11 13 17 19 23 29 31 37 41 1∗ 2∗ 3∗ 4∗ 9000 13400 15500 17600 19600 19800 19900 19900 19900 19900 19900 19900 19900 0 0 0 0 0 0 0 0 0 0 0 0 0 9000 13400 15500 17600 19600 19800 19900 19900 19900 19900 19900 19900 19900 0 0 0 0 0 0 0 0 0 0 0 0 0 prime number 43 47 53 59 61 67 71 73 79 83 89 97 1∗ 2∗ 3∗ 4∗ 20000 20000 20000 20000 20000 20000 20000 20000 19900 19900 20000 20000 0 0 0 0 0 0 0 0 0 0 0 0 20000 20000 20000 20000 20000 20000 20000 20000 19900 19900 20000 20000 0 0 0 0 0 0 0 0 0 0 0 0 [1,z] [−1,z] 1∗ (resp.,3∗ )is the number of elements which have the periodic Φ3A (resp.,Φ3A 26 ) continued fraction expansion and 2∗ (resp.,4∗ ) is the number of elements such [−1,z] [1,z] that their heights of some remainders with respect to Φ3A (resp.,Φ3A ) exceed 10300 and their successive iterations terminate. 9 Conjecture We give the following conjecture which is supported by the numerical experiments. Conjecture. Let p be any prime number, and K be any finite extension of Q with K ⊂ Qp . Let s + 1 be its degree over Q and z ∈ K be any element with the condition (H) and K = Q(z). For every α = (α1 , . . . , αs ) ∈ D which satis[−1,z] fies that 1, α1 , . . . , αs are linearly independent over Q α has a periodic Φ3A continued fraction expansion. References [1] H.Bekki;On periodicity of geodesic continued fractions,J. Number Theory 177 (2017), 181-210. [2] J. Browkin; Continued fractions in local fields. II. Math. Comp. 70 (2001), no. 235, 1281-1292. [3] P. Bundschuh; p-adische Kettenbrüche und Irrationalität p-adischer Zahlen, Elem. Math. 32 (1977), no. 2, 36-40. [4] C.Bauer, A.Frink, R.Kreckel; Introduction to the GiNaC framework for symbolic computation within the C++ programming language, J. Symbolic Comput. 33 (2002), no. 1, 112. [5] E.Dubois, R.Paysant-Le Roux; Algorithme de JacobiPerron dans les extensions cubiques. C. R. Acad. Sci. Paris Sr. A-B 280(1975), 183186. [6] T.Ooto; Transcendental p-adic doi:10.1007/s00209-017-1859-2. continued fractions, Math.Z(2017), [7] A. A. Ruban; Certain metric properties of p-adic numbers, (Russian), Sibirsk. Mat. Zh. 11 (1970), 222-227. [8] A. Saito, J.-I.Tamura, S. Yasutomi; p-adic continued fractions and Lagrange’s theorem, preprint [9] A. Saito, J.-I.Tamura, S. Yasutomi; Convergence of multi-dimensional p-adic continued fractions, preprint(2017). 27 [10] A. Saito, A.Yamaguchi; Pseudorandom number generation using chaotic true orbits of the Bernoulli map, Chaos 26, 063122 (2016) [11] T. Schneider; Über p-adische Kettenbrüche, Symp. Math. 4 (1968/69), 181189. [12] J.-I.Tamura; A p-adic phenomenon related to certain integer matrices, and p-adic values of a multidimensional continued fraction, in: Summer School on the Theory of Uniform Distribution, RIMS Kôkyûroku Bessatsu B29 (2012), 1-40. [13] B. M. M.de Weger; Periodicity of p-adic continued fractions. Elem. Math. 43 (1988), no. 4, 112-116. 28