Download arXiv:1705.06122v1 [math.NT] 17 May 2017

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