Next Article in Journal
Hybrid Fruit-Fly Optimization Algorithm with K-Means for Text Document Clustering
Next Article in Special Issue
On Some Properties of the Limit Points of (z(n)/n)n
Previous Article in Journal
Nadler’s Theorem on Incomplete Modular Space
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Parity of the Order of Appearance in the Fibonacci Sequence

Department of Mathematics, Faculty of Science, University of Hradec Králové, 500 03 Hradec Králové, Czech Republic
Mathematics 2021, 9(16), 1928; https://doi.org/10.3390/math9161928
Submission received: 15 July 2021 / Revised: 4 August 2021 / Accepted: 12 August 2021 / Published: 13 August 2021
(This article belongs to the Special Issue New Insights in Algebra, Discrete Mathematics and Number Theory II)

Abstract

:
Let ( F n ) n 0 be the Fibonacci sequence. The order of appearance function (in the Fibonacci sequence) z : Z 1 Z 1 is defined as z ( n ) : = min { k 1 : F k 0 ( mod n ) } . In this paper, among other things, we prove that z ( n ) is an even number for almost all positive integers n (i.e., the set of such n has natural density equal to 1).

1. Introduction

Let ( F n ) n be the sequence of Fibonacci numbers which is defined by the recurrence F n + 2 = F n + 1 + F n , with initial values F 0 = 0 and F 1 = 1 . For any integer n 1 , the order of apparition (or the rank of appearance) of n (in the Fibonacci sequence), denoted by z ( n ) , is the positive index of the smallest Fibonacci number which is a multiple of n (more precisely, z ( n ) : = min { k 1 : n F k } ). The arithmetic function z : Z 1 Z 1 is well defined (as can be seen in Lucas [1], p. 300) and, in fact, z ( n ) 2 n is the sharpest upper bound (as can be seen in [2]). A few values of z ( n ) (for n [ 1 , 50 ] ) can be found in Table 1 (see the OEIS [3] sequence A001177 and for more facts on properties of z ( n ) see, e.g., [4,5,6,7,8,9]).
Recall that the natural density of A Z > 0 is the following limit (if it exists):
δ ( A ) : = lim x # A ( x ) x ,
where A ( x ) : = A [ 1 , x ] for x > 0 . In 2020, Trojovský [10] proved that natural density of the set { n 1 : z ( n ) < ϵ n } is 1, for all ϵ > 0 (which generalized a result due to Marques [11] about lim inf n z ( n ) / n ). Let f : Z 1 Z 1 be an arithmetic function. We denote by E f the set all n for which f ( n ) is an even number (i.e., E f : = { n 1 : f ( n ) 0 ( mod 2 ) } ). For example, if ϕ ( n ) and d ( n ) are the Euler totient function and the number-of-divisors function, respectively, then E ϕ = Z 3 and E d = Z 1 { n 2 : n 1 } . In particular, both sets E ϕ and E d have natural density equal to 1. Another well-known arithmetic function related to Fibonacci numbers is the Pisano period π : Z 1 Z 1 for which π ( n ) is the smallest period of ( F k ( mod n ) ) k . The first few values of π ( n ) (for n [ 1 , 50 ] ) are in Table 2 (see sequence A001175 in OEIS):
Since F i = i , for i { 1 , 2 } , the Pisano period may be defined as π ( n ) : = min { k 1 : F k 0 ( mod n ) and F k + 1 1 ( mod n ) } . Thus, for n 3 , π ( n ) is the order of the Fibonacci Q-matrix:
Q = 1 1 1 0
in the general linear group G L 2 ( Z n ) (i.e., the invertible 2 × 2 matrices modulo n). In particular, 1 = det ( Q π ( n ) ) = ( det ( Q ) ) π ( n ) = ( 1 ) π ( n ) , yielding that π ( n ) is even, for all n 3 . In conclusion, E π = Z 3 .
Observe that π ( n ) and z ( n ) have similar definitions (in fact, these functions are strongly connected, as can be seen in Lemma 3). However, they have very distinct behavior related to their parity. For example, Z 1 E z is an infinite set (since z ( 5 k ) = 5 k is an odd number, for all k 0 ).
Therefore, the aim of this paper is to study some arithmetical and analytical properties of E z : = { n 1 : z ( n ) 0 ( mod 2 ) } . Our first result provides three infinite arithmetic progressions belonging to this set. More precisely,
Theorem 1.
We have that { 3 n , 4 n 1 , 4 n } E z , for all n 1 .
Remark 1.
We remark that the arithmetic progressions of Theorem 1 are the largest ones (in the sense of natural density) belonging to E z . In fact, for other arithmetic progressions of density of at least 1 / 4 , we have infinitely many terms which are outside of E z . Some of these terms can be made explicit. For example:
  • If n = 5 k , then z ( 2 n ) = 3 · 5 k , for all k 0 .
  • If n = ( 5 k 1 ) / 2 , then z ( 2 n + 1 ) = 5 k , for all k 0 .
  • If n = ( 5 2 k 1 ) / 3 , then z ( 3 n + 1 ) = 5 2 k , for all k 0 .
  • If n = ( 5 2 k + 1 2 ) / 3 , then z ( 3 n + 2 ) = 5 2 k + 1 , for all k 0 .
  • If n = ( 5 k 1 ) / 4 , then z ( 4 n + 1 ) = 5 k , for all k 0 .
  • If n = ( 5 k 1 ) / 2 , then z ( 4 n + 2 ) = 3 · 5 k , for all k 0 .
We point out that Theorem 1 implies, in particular, that the natural density of E z is at least 1 / 2 (since ( 4 n 1 ) n and ( 4 n ) n are disjoint arithmetic progressions with natural density equal to 1 / 4 ). In fact, it is possible to obtain a better lower bound for δ ( E z ) , by noting that E z also contains ( 3 n ) n which has density 1 / 3 . Since the last sequence is not disjoint from the previous ones, we may use the inclusion–exclusion principle to deduce that the natural density of E z is at least
1 4 + 1 4 + 1 3 1 12 1 12 = 2 3 0.66 .
Table 3 suggests that for almost all positive integers n, one has that z ( n ) is even (i.e., one expects that δ ( E z ) = 1 ). Note that a sieve method of searching by arithmetic progressions lying in E z does not seem to reach the value 1.
However, we were able to surmount this obstacle to confirm the expectation:
Theorem 2.
There exists an effective computable constant c > 0 such that
# E z ( x ) x c x ( log x ) 1 / 2 ,
for all x > 6.3 × 10 27 . In particular, one has that δ ( E z ) = 1 .
Remark 2.
Note that all previously provided infinite sequences belonging to Z 1 E z have exponential growth (e.g., ( 5 n ) n ) . In fact, the previous theorem ensures, in particular, that no infinite arithmetic progression can be within Z 1 E z (since its natural density is 0 ) .
The proof of the theorems combines Diophantine properties of z ( n ) with analytical tools concerning primes in arithmetic progressions.

2. Auxiliary Results

In this section, we shall present some results which will be very important tools in the proof. The first ingredient is related to the value of z ( p k ) , for a prime number p and k 1 :
Lemma 1
(Theorem 2.4 of [12]). We have that z ( 2 k ) = 3 · 2 k 1 , for all k 2 , and z ( 3 k ) = 4 · 3 k 1 , for all k 1 . In general, it holds that
z ( p k ) = p max { k e ( p ) , 0 } z ( p ) ,
where e ( p ) : = max { k 0 : p k F z ( p ) } .
The next result provides some divisibility properties of the Pisano period for prime numbers.
Lemma 2
(Theorem 2.2 of [12]). Let p be a prime number. We have that
(i)
If p ± 1 ( mod 5 ) , then π ( p ) divides p 1 .
(ii)
If p ± 2 ( mod 5 ) , then π ( p ) divides 2 ( p + 1 ) . Furthermore, π ( p ) = 2 ( p + 1 ) / t , for some odd number t.
Observe that F π ( n ) F 0 0 ( mod n ) and then z ( n ) divides π ( n ) . Our next tool provides a characterization of the quotient π ( n ) / z ( n ) . In fact,
Lemma 3
(Theorem 1 of [13]). We have that π ( n ) / z ( n ) { 1 , 2 , 4 } , for all n 1 . Moreover, π ( n ) = 4 z ( n ) if and only if z ( n ) is odd.
The next tool is a kind of “formula” for z ( n ) depending on z ( p a ) for all primes p dividing n. The proof of this fact can be found in [14].
Lemma 4
(Theorem 3.3 of [14]). Let n > 1 be an integer with prime factorization n = p 1 a 1 p k a k . Then
z ( n ) = l c m ( z ( p 1 a 1 ) , , z ( p k a k ) ) .
In general, it holds that
z ( l c m ( m 1 , , m k ) ) = l c m ( z ( m 1 ) , , z ( m k ) ) .
In order to prove Theorem 2, we need an analytic tool related to the profusion of integers having factorization allowing only some classes of primes. The following notation will be used throughout this work: let P be the set of prime numbers and for an integer q 2 , set P ( a , q ) as the set of all prime numbers of the form a + k q , for some integer k 0 (Dirichlet’s theorem on arithmetic progressions) that ensures that P ( a , q ) is an infinite set whenever gcd ( a , q ) = 1 ). Let N a , q = { n 1 : p n p P ( a , q ) } be the set of all positive integers whose prime factors lie exclusively in P ( a , q ) . Additionally, denote β : = 1 / ϕ ( q ) and
G a , b ( s ) : = ζ ( s ) β p P ( a , q ) ( 1 p s ) 1 ,
which has an analytic continuation to a neighborhood of s = 1 . Here, as usual, ζ ( s ) denotes the Riemann zeta function.
Our last auxiliary lemma is a particular case of a theorem due to Chang and Martin [15]. We state this more precisely in the next lemma.
Lemma 5
(Theorem 3.4 of [15]). For any integer q 3 , there exists a positive absolute constant C such that uniformly for q ( log x ) 1 / 3 , we have
# { n x : n N a , q } = x ( log x ) 1 β G a , q ( 1 ) Γ ( β ) + O C ( log x ) 1 / 4 ,
where Γ ( z ) = 0 t z 1 e t d t is the Gamma function.
Now, we are ready to deal with the proof of the theorems.

3. The Proofs

3.1. The Proof of Theorem 1

  • The Case z ( 3 n )
Proof. 
First, in order to use Lemma 4, one may write n = 3 k m , where k 0 and gcd ( m , 3 ) = 1 . So, Lemma 4 yields that
z ( 3 n ) = z ( 3 k + 1 m ) = lcm ( z ( 3 k + 1 ) , z ( m ) ) = lcm ( 4 · 3 k , z ( m ) ) ,
where we used Lemma 1. Hence we deduce that 4 divides z ( 3 n ) . □
  • The Case z ( 4 n )
Proof. 
Let us write n = 2 k m , where m is odd and k 0 . Thus, by Lemma 4, we have
z ( 4 n ) = z ( 2 k + 2 m ) = lcm ( z ( 2 k + 2 ) , z ( m ) ) = lcm ( 3 · 2 k + 1 , z ( m ) ) ,
where we used Lemma 1. Since k + 1 1 , then 2 divides lcm ( 3 · 2 k + 1 , z ( m ) ) = z ( 4 n ) , as desired. □
  • The Case z ( 4 n + 3 )
This case is more delicate, since 4 n + 3 is a prime number for infinitely many values of n. To overcome this difficulty, we deal first with these prime values.
Claim 1.
If p = 4 k + 3 is a prime number, then z ( p ) is even. The proof splits into two cases concerning the residues of p modulo 5. In fact, if p ± 1 ( mod 5 ) , then Lemma 2 (i) implies that π ( p ) divides p 1 . Striving for a contradiction, suppose that z ( p ) is odd. Then, by Lemma 3, we have that π ( p ) = 4 z ( p ) leading to the absurdity that 4 z ( p ) divides p 1 = 4 n + 2 . Thus, z ( p ) is an even number. Now, for the case in which p ± 2 ( mod 5 ) , one has, by Lemma 2 (ii), that π ( p ) = 2 ( p + 1 ) / t , for some odd integer t. Thus, by Lemma 3, we obtain
2 i z ( p ) = π ( p ) = 2 ( p + 1 ) t = 8 ( n + 1 ) t ,
for some i { 0 , 1 , 2 } . Thus, z ( p ) = 2 3 i ( n + 1 ) / t is an even number, since 3 i 1 and t 1 ( mod 2 ) . This finishes the proof of the claim.
Proof. 
Let 4 n + 3 = p 1 a 1 p k a k be the prime factorization of the composite number 4 n + 3 , with a i 1 . By Lemma 4, we obtain
z ( 4 n + 3 ) = lcm ( z ( p 1 a 1 ) , , z ( p k a k ) ) .
Since the product of numbers 1 ( mod 4 ) is also 1 ( mod 4 ) , then at least one of the prime factors of 4 n + 3 , say p 1 , must be of the form 4 + 3 . Therefore, by Lemma 1, we can rewrite (2) as
z ( 4 n + 3 ) = lcm ( p 1 max { a 1 e ( p 1 ) , 0 } z ( p 1 ) , , z ( p k a k ) ) .
However, p 1 = 4 + 3 is a prime number and the previous claim implies that z ( p 1 ) is even and so is z ( 4 n + 3 ) (since z ( p 1 ) divides z ( 4 n + 3 ) ). This completes the proof. □

3.2. The Proof of Theorem 2

Proof. 
We have that
# E z ( x ) = # { n x : z ( n ) 0 ( mod 2 ) } .
Note that since z ( n ) is a multiple of z ( p ) , for all prime p in the factorization of n (by Lemma 4), then a sufficient condition for z ( n ) to be even is for n to have some prime factor of the form 4 k + 3 (since z ( 4 k + 3 ) is even, by Theorem 1). Therefore,
# E z ( x ) = # { n x : z ( n ) 2 Z } # { n x : p n with p 3 ( mod 4 ) } = x # { n x : p n p 1 ( mod 4 ) } .
Since
{ n x : p n p 1 ( mod 4 ) } = { n x : x N 1 , 4 } ,
we can apply Lemma 5 in order to get an upper bound for the cardinality of the previous set. Thus, Lemma 5 (for β = 1 / ϕ ( 4 ) = 1 / 2 ) implies the existence of a positive absolute constant C such that
# { n x : n N 1 , 4 } = x ( log x ) 1 / 2 G 1 , 4 ( 1 ) π + O C ( log x ) 1 / 4 ,
for all x e 64 6.2 × 10 27 (where we used that Γ ( 1 / 2 ) = π ). Moreover, we have that
G 1 , 4 ( s ) : = ζ ( s ) 1 / 2 p P ( 1 , 4 ) ( 1 p s ) 1
and
G 3 , 4 ( s ) : = ζ ( s ) 1 / 2 p P ( 3 , 4 ) ( 1 p s ) 1 .
Since both sets P ( 1 , 4 ) and P ( 3 , 4 ) have the same density inside the set of all primes (namely, 1 / 2 ), then G 1 , 4 ( s ) G 3 , 4 ( s ) (i.e., G 1 , 4 ( s ) G 3 , 4 ( s ) and G 3 , 4 ( s ) G 1 , 4 ( s ) ). Now, we multiply (6) and (7) to arrive at
G 1 , 4 ( s ) G 3 , 4 ( s ) = ζ ( s ) 1 p P ( 1 , 4 ) ( 1 p s ) 1 · p P ( 3 , 4 ) ( 1 p s ) 1 = ζ ( s ) 1 p P ( 1 , 2 ) ( 1 p s ) 1 = p P ( 1 p s ) · p P ( 1 , 2 ) ( 1 p s ) 1 = 1 1 2 s ,
where we used the Euler product ζ ( s ) = p P ( 1 p s ) 1 (see [16], p. 39) and that P ( 1 , 2 ) = P { 2 } . In particular, we have that G 1 , 4 ( 1 ) G 3 , 4 ( 1 ) = 1 / 2 and so G 1 , 4 ( 1 ) = O ( 1 ) . Therefore, there exists a positive constant c such that (5) becomes
# { n x : n N 1 , 4 } c x ( log x ) 1 / 2 .
Finally, we deduce by combining (3), (4), and (8) that
# E z ( x ) x c x ( log x ) 1 / 2
holds for all x > 6.3 · 10 27 > e 64 .
Additionally, since
1 δ ( E z ) = lim x # E z ( x ) x lim x 1 c ( log x ) 1 / 2 = 1 ,
we obtain that the natural density of E z is equal to 1. The proof is then complete. □

4. Further Comments

We close this paper by offering one conjecture and two questions for further research. The natural generalization of the study presented in this work is to consider a similar problem by replacing 2 with another positive integer. For example, is z ( n ) a multiple of 4, for almost all positive integers n?
To study this problem, we start by defining (for an integer m 2 ) E z ( m ) as the set of all n Z 1 for which z ( n ) is a multiple of m (i.e., E z ( m ) : = { n 1 : z ( n ) 0 ( mod m ) } ). In particular, E z ( 2 ) = E z and, by Theorem 2, one has that almost all positive integers belong to E z ( 2 ) .
Note that, by Lemma 4, z ( p ) divides z ( m ) whenever p divides m. In particular,
# E z ( z ( p ) ) ( x ) x p
yielding that the natural density of E z ( z ( p ) ) is at least 1 / p . In the case in which p = 3 , we infer that δ ( E z ( 4 ) ) 1 / 3 (observe that z ( 3 ) = 4 ). Moreover, we believe that
Conjecture 1.
The natural density of E z ( 4 ) is equal to 1.
The main difficulty in proving the last statement is to find some infinite arithmetic progression (containing infinitely many prime numbers) lying completely in E z ( 4 ) (in the spirit of Theorem 1). In general, we conclude by posing the following questions:
Question 1.
Are there infinitely many prime numbers p, for which δ ( E z ( p ) ) = 1 ?
Question 2.
Let m 1 be an integer. Is it possible to provide an explicit positive lower bound for # E z ( m ) ( x ) ?

5. Conclusions

In this paper, we study a Diophantine problem related to the arithmetic function z : Z 1 Z 1 , defined as z ( n ) : = min { k 1 : n F k } (the so-called order of appearance in the Fibonacci sequence). The problem is to establish some properties of E z : = { n 1 : 2 z ( n ) } . Indeed, we provide three infinite arithmetic progressions lying entirely within E z (which implies, in particular, that the natural density of E z is at least 2 / 3 ). Furthermore, we provide a lower bound for # ( E z [ 1 , x ] ) (for all sufficiently large x) which, in particular, confirms the expectation: the natural density of E z is 1. We close the work with some comments, questions, and a conjecture for further research. The proof combines arithmetical and analytical tools in number theory.

Funding

The research was supported by the Excellence Project PřF UHK No. 2213/2021–2022., University of Hradec Králové, Czech Republic.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author thanks the University of Hradec Králové for support.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Lucas, E. Théorie des fonctions numériques simplement périodiques. Am. J. Math. 1878, 1, 289–321. [Google Scholar] [CrossRef]
  2. Sallé, H.J.A. Maximum value for the rank of apparition of integers in recursive sequences. Fibonacci Q. 1975, 13, 159–161. [Google Scholar]
  3. Sloane, N.J.A. The On-Line Encyclopedia of Integer Sequences. Available online: http://www.research.att.com/~njas/sequences/ (accessed on 5 July 2021).
  4. Marques, D. Fixed points of the order of appearance in the Fibonacci sequence. Fibonacci Q. 2012, 50, 346–352. [Google Scholar]
  5. Somer, L.; Křížek, M. Fixed points and upper bounds for the rank of appearance in Lucas sequences. Fibonacci Q. 2013, 51, 291–306. [Google Scholar]
  6. Sun, Z.H.; Sun, Z.W. Fibonacci numbers and Fermat’s last theorem. Acta Arith. 1992, 60, 371–388. [Google Scholar] [CrossRef] [Green Version]
  7. Luca, F.; Tron, E. The distribution of self-Fibonacci divisors. In Advances in the Theory of Numbers; Fields Institute Communications; Springer: New York, NY, USA, 2015; Volume 77, pp. 149–158. [Google Scholar]
  8. Trojovský, P. On Diophantine equations related to order of appearance in Fibonacci sequence. Mathematics 2019, 7, 1073. [Google Scholar] [CrossRef] [Green Version]
  9. Trojovská, E. On the Diophantine Equation z(n) = (2 − 1/k)n Involving the Order of Appearance in the Fibonacci Sequence. Mathematics 2020, 8, 124. [Google Scholar] [CrossRef] [Green Version]
  10. Trojovský, P. Some problems related to the growth of z(n). Adv. Differ. Equ. 2020, 2020, 270. [Google Scholar] [CrossRef]
  11. Marques, D. Sharper upper bounds for the order of appearance in the Fibonacci sequence. Fibonacci Q. 2013, 51, 233–238. [Google Scholar]
  12. Fulton, J.D.; Morris, W.L. On arithmetical functions related to the Fibonacci numbers. Acta Arith. 1969, 16, 105–110. [Google Scholar] [CrossRef] [Green Version]
  13. Vinson, J. The relation of the period modulo to the rank of apparition of m in the Fibonacci sequence. Fibonacci Q. 1963, 1, 37–46. [Google Scholar]
  14. Renault, M. Properties of the Fibonacci Sequence Under Various Moduli. Master’s Thesis, Wake Forest University, Winston-Salem, NC, USA, 1996. Available online: http://webspace.ship.edu/msrenault/fibonacci/FibThesis.pdf (accessed on 3 August 2021).
  15. Chang, B.; Martin, G. The smallest invariant factor of the multiplicative group. Int. J. Number Theory 2020, 16, 1377–1405. [Google Scholar] [CrossRef] [Green Version]
  16. De Koninck, J.-M.; Luca, F. Analytic Number Theory-Exploring the Anatomy of Integers; American Mathematical Society: Providence, RI, USA, 2012. [Google Scholar]
Table 1. Values of z ( n ) for n from 1 to 50.
Table 1. Values of z ( n ) for n from 1 to 50.
nz(n)nz(n)nz(n)nz(n)nz(n)
11111021831304120
231212223032244224
34137232433204344
46142424123494430
551520252535404560
6121612262136124624
78179273637194716
861812282438184812
9121918291439284956
10152030306040305075
Table 2. Values of π ( n ) for n from 1 to 50.
Table 2. Values of π ( n ) for n from 1 to 50.
n π ( n ) n π ( n ) n π ( n ) n π ( n ) n π ( n )
111110211631304140
231224223032484248
381328234833404388
461448242434364430
520154025100358045120
6241624268436244648
7161736277237764732
8121824284838184824
92419182914395649112
1060206030120406050300
Table 3. Quantity and proportion of arguments for which z ( n ) is even.
Table 3. Quantity and proportion of arguments for which z ( n ) is even.
x10205010020050010003000
# E z ( x ) 61439811634258622632
# E z ( x ) x 0.60.70.780.810.8150.850.8620.877
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Trojovský, P. On the Parity of the Order of Appearance in the Fibonacci Sequence. Mathematics 2021, 9, 1928. https://doi.org/10.3390/math9161928

AMA Style

Trojovský P. On the Parity of the Order of Appearance in the Fibonacci Sequence. Mathematics. 2021; 9(16):1928. https://doi.org/10.3390/math9161928

Chicago/Turabian Style

Trojovský, Pavel. 2021. "On the Parity of the Order of Appearance in the Fibonacci Sequence" Mathematics 9, no. 16: 1928. https://doi.org/10.3390/math9161928

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop