1,746
Views
0
CrossRef citations to date
0
Altmetric
ORIGINAL ARTICLE

On the Eight Levels theorem and applications towards Lucas-Lehmer primality test for Mersenne primes, I

ORCID Icon
Pages 267-284 | Received 30 Dec 2022, Accepted 15 Apr 2023, Published online: 11 May 2023

Abstract

Lucas-Lehmer test is the current standard algorithm used for testing the primality of Mersenne numbers, but it may have limitations in terms of its efficiency and accuracy. Developing new algorithms or improving upon existing ones could potentially improve the search for Mersenne primes and the understanding of the distribution of Mersenne primes and composites. The development of new versions of the primality test for Mersenne numbers could help to speed up the search for new Mersenne primes by improving the efficiency of the algorithm. This could potentially lead to the discovery of new Mersenne primes that were previously beyond the reach of current computational resources. The current paper proves what the author called the Eight Levels Theorem and then highlights and proves three new different versions for Lucas-Lehmer primality test for Mersenne primes and also gives a new criterion for Mersenne compositeness.

1. Introduction

Primes of special form have been of perennial interest (Guy, Citation2004). Among these, the primes of the form 2p1 which are called Mersenne prime. It is outstanding in their simplicity.

  • Clearly, if 2p1 is prime then p is prime. The number 2p1 is Mersenne composite if p is prime but 2p1 is not prime. For example, the prime number 271=127 is called Mersenne prime. However, the number 2111=2047=23×89 is called Mersenne composite (see Cambraia, Knapp, Lemos, Moriya, & Rodrigues, Citation2022; Ferradi & Xagawa, Citation2020; Jo & Kim, Citation2022; Kharaghani & Suda, Citation2020; Skula, Citation2019; Zheng & Yang, Citation1985).

  • Mathematics is kept alive by the appearance of challenging unsolved problems. The current paper gives new expansions related to the following two major open questions in number theory: Are there infinitely many Mersenne primes? Are there infinitely many Mersenne composite? Many mathematicians believe that there are infinitely many Mersenne primes but a proof of this is still one of the major open problems in number theory (see Deza, Citation2021; Gallier, Citation2011; Guy, Citation2004; Effinger & Mullen, Citation2021; Kraft & Washington, Citation2018; Kundu & Mazumder, Citation2022).

  • Mersenne primes have a close connection to perfect numbers, which are numbers that are equal to the sum of their proper divisors. It is known that Euclid and Euler proved that a number N is even perfect number if and only if N=2p1(2p1) for some prime p, and 2p1 is prime. Euclid proved only that this statement was sufficient. Euler, 2000 years later, proved that all even perfect numbers are of the form 2p1(2p1) where 2p1 is a Mersenne prime (see Gallier, Citation2011; Wells, Citation1987). Thus the theorems of Euclid and Euler characterize all even perfect numbers, reducing their existence to that of Mersenne primes.

  • The odd perfect numbers are quite a different story (see Dickson, Citation1919; Ochem & Rao, Citation2012; Wells, Citation1987). An odd perfect number is a hypothetical number that is both an odd integer and a perfect number, meaning that it is equal to the sum of its proper divisors. No odd perfect number has been discovered, and it is not known whether any exist. The search for odd perfect numbers has been ongoing for centuries and has involved many of the greatest mathematicians in history. It is still unknown whether there is any odd perfect number. Recently (Ochem & Rao, Citation2012), showed that odd perfect numbers, if exist, are greater than 101500

  • There are practical reasons for seeking out bigger and bigger primes. Very big primes are crucial to the most widely used systems for encrypting data such as those that underpin all online banking and shopping. Based on the properties of Mersenne primes, instead of polynomial rings, Aggarwal, Joux, Prakash, and Santha described an elegant public-key encryption (see Aggarwal, Joux, Prakash, & Santha, Citation2018).

  • The discovery of Mersenne primes has been an active area of research for centuries, and the search for larger and more complex prime numbers continues to this day. The discovery of new Mersenne primes is a significant achievement in the field of mathematics, as these numbers have interesting mathematical properties and are very useful in various applications, such as cryptography and computer science. Using the Lucas-Lehmer test, the computation was carried out at the UCLA Computing Facility, it was reported that 242531 and 244231 are primes (Hurwitz, Citation1962). Discovering Mersenne primes is a challenging and time-consuming task, and the discovery of new Mersenne primes often represents a significant achievement in the field of mathematics. The Great Internet Mersenne Prime Search (GIMPS), which is a collaborative project aimed at discovering new Mersenne primes, has discovered the largest known prime number

282,589,9331, having 24,862,048 digits (see Deza, Citation2021; GIMPS, Citation2016). The GIMPS project has been instrumental in discovering the majority of the known Mersenne primes, and it continues to search for new ones. GIMPS is a fascinating example of how the power of distributed computing can be used to advance mathematical knowledge and uncover new discoveries. Each new discovery sets a new record for the largest known prime number. Despite its success in discovering large prime numbers, the Great Internet Mersenne Prime Search (GIMPS) has some limitations. It’s worth noting that Mersenne primes are rare and become increasingly difficult to find as their size increases. Surprisingly, as of March 2023, only 51 Mersenne primes are known. While 51 Mersenne primes are currently known, it is not known whether there are any more beyond those that have been discovered. The question of whether there are infinitely many Mersenne primes is one of the major unsolved problems in mathematics. Whether the list of known Mersenne primes is finite or infinite is an open question in mathematics, and it remains an area of active research (see Deza, Citation2021; GIMPS, Citation2016; Kundu & Mazumder, Citation2022).

1.1. Lucas-Lehmer primality test

Mersenne primes have some special properties that make them useful in certain types of computations, particularly in the field of cryptography and computing. One of these properties is that they are easy to test for primality using the Lucas-Lehmer primality test. This makes it easier to find large prime numbers, which are important for secure cryptographic systems. The Lucas-Lehmer primality test is a primality test specifically designed for testing the primality of Mersenne numbers (i.e. numbers of the form 2p1). It is a fast and efficient test that can quickly determine whether a given Mersenne number is prime or composite. It is well-known the following theorems (see Deza, Citation2021; Gallier, Citation2011; Kraft & Washington, Citation2018):

Theorem 1.

(Lucas-Lehmer primality test)

2p1 is Mersenne prime if and only if (1) 2n1|(1+3)n+(13)n(1) where n:=2p1.

1.2. Euclid-Euler-Lucas-Lehmer association

Perfect numbers have fascinated mathematicians for centuries, and many properties and patterns have been discovered about them. One of the most famous results about perfect numbers is that every even perfect number can be written in the form 2p1(2p1), where 2p1 is a Mersenne prime. Moreover, every Mersenne prime of the form 2p1 gives rise to an even perfect number through this formula. The following theorem tells us that there is a strong association between perfect numbers and Mersenne primes.

Theorem 2.

(Euclid-Euler-Lucas-Lehmer)

A number N is even perfect number if and only if N=2p1(2p1) for some prime p, and (2) 2n1|(1+3)n+(13)n(2) where n:=2p1.

For example, the first few even perfect numbers are: (3) 6=221(221),28=231(231),496=251(251),8128=271(271),  (3) and the corresponding Mersenne primes are: (4) 221=3,231=7,251=31,271=127.(4)

1.3. Notations

For a natural number n, we define δ(n)=n(mod 2). For an arbitrary real number x, x2 is the highest integer less than or equal x2. We also need the following notation λ=01(n2(4λ)2):=1.

For simplicity, we sometimes write λn2(4λ)2 for λ[n2(4λ)2].

2. The purpose of the current paper

The aim of this paper is to study some arithmetical properties of the coefficients of the expansion (5) xn+yn=k=0n2Ψk(n)(xy)n2k(x2+y2)k,(5) for n0,2,4,6(mod 8). Then we study the expansion (6) xn+ynx+y=k=0n2Ψk(n)(xy)n2k(x2+y2)k,(6) for n1,3,5,7(mod 8).

2.1. Reasons for why the expansions of 5 and 6 are interesting

We show in the current paper that the numbers of the form (7) λ(n2(4λ)2)(7) arise up naturally in the coefficients of the expansions (5, 6), and enjoy some unexpected new interesting arithmetical properties which could be helpful in the study of Mersenne primes and Mersenne composites.

2.2. New results

We state and prove the following new results in order:

  1. The Eight Levels theorem

  2. The first new version for Lucas-Lehmer primality test

  3. The second new version for Lucas-Lehmer primality test

  4. The third new version for Lucas-Lehmer primality test

  5. New version for Euclid-Euler-Lucas-Lehmer association

  6. New criteria for compositeness of Mersenne numbers

  7. New combinatorial identities

We end the paper with a discussion about further potential investigations of how this new versions for Lucas-Lehmer primality test may provide a better theoretical understanding of the two major questions about Mersenne numbers; whether there are infinitely many Mersenne primes or Mersenne composites.

3. Summary for the main new results of the paper

This paper proves the following new theorems.

3.1. The first new version for Lucas-Lehmer primality test

The following result proposes a new double-indexed recurrence relation for the Lucas-Lehmer test.

Theorem 3.

Given prime p5. 2p1 is prime if and only if (8) 2n1|k evenk=0,n2ϕk(n)(8) where n:=2p1,ϕk(n) are defined by the double index recurrence relation ϕk(m)=4 ϕk1(m2)ϕk(m4) and the initial boundary values satisfy (9) ϕ0(m)={+2m0(mod 8)0m±2(mod 8)2m4(mod 8),ϕ1(m)={+2 mm2(mod 8)0m0, 4(mod 8)2 mm6(mod 8).  (9)

3.2. The second new version for Lucas-Lehmer primality test

The following result proposes a new explicit sum of products of difference of squares for the Lucas-Lehmer test.

Theorem 4.

Given prime p5,n:=2p1. The number 2p1 is prime if and only if (10) 2n1|k evenk=0,n2(1)k2 λ=0k21n2(4λ)2k!.(10)

3.3. The third new version for Lucas-Lehmer primality test

The following result proposes a new nonlinear recurrence relation for the Lucas-Lehmer test.

Theorem 5.

Given prime p5. The number 2p1 is prime if and only if (11) 2n1|kevenk=0,n2ϕk(n)(11) where n:=2p1,ϕk(n) are generated by the double index recurrence relation ϕk(n)ϕk2(n)=(2k4)2n2k(k1), and we can choose either of the following initial values to generate ϕk(n) from the starting term ϕ0(n) or the last term ϕn2(n): (12) ϕ0(n)=+2,ϕn2(n)=2n.(12)

3.4. New version for Euclid-Euler-Lucas-Lehmer association

Theorem 6.

A number N is even perfect number if and only if N=2p1(2p1) for some prime p, and (13) 2n1|kevenk=0,n2(1)k2λ=0k21[n2(4λ)2]k!(13) where n:=2p1.

3.5. New criteria for compositeness of Mersenne numbers

Theorem 7.

(Criteria for compositeness of Mersenne numbers)

Given prime p5. The number 2n1=2p1 is Mersenne composite number if (14) 2n-1kevenk=0,n2λ=0k21[(4λ)2-41]k!.(14)

3.6. New combinatorial identities

Theorem 8.

(combinatorial identities)

For any natural number n, the following combinatorial identities are correct n0(mod 4)¯4n2(n2)!=2λ=0n44[n2-(4λ)2]n1(mod 4)¯4n2(n2)!=λ=1n14[(n+1)2-(4λ2)2]n2(mod 4)¯4n2(n2)!=2nλ=1n24[n2-(4λ2)2]n3(mod 4)¯4n2(n2)!=(n+1)λ=1n34[(n+1)2-(4λ)2]

4. Discussion for the proposed method’s theoretical analysis

4.1. Algebraically independent polynomials

Algebraically independent polynomials are important in several areas of mathematics and its applications, including algebraic geometry, commutative algebra, and number theory. Algebraic independent polynomials are important because they provide a way to study and understand the relationships between different algebraic objects, such as numbers, functions, and algebraic structures (Perron, Citation1932). Two polynomials in two variables, say h(x, y) and g(x, y), are said to be algebraically independent over the field of rational numbers if there exist a polynomial f with rational coefficients such that f(h(x,y),g(x,y))=0, then this entails that f = 0. The polynomials xy and x2+y2 are algebraically independent over the field of rational numbers. Therefore, if we have the following identity, with integer coefficients μk(n), (15) 0=k=0n2μk(n)(xy)n2k(x2+y2)k,(15) then this entails that all the coefficients vanish; which means μk(n)=0.

4.2. Symmetric polynomials

Symmetric polynomials are polynomials that remain unchanged under the permutation of their variables. For example, the polynomial x2+y2+z2 is symmetric because it remains the same if we swap the variables x, y, and z. For any natural number n, xn+yn is symmetric polynomial. Then we have integer coefficients fk(n), from the fundamental theorem of symmetric polynomials (Perron, Citation1932), that satisfy (16) xn+yn=k=0n2fk(n)(x+y)n2k(x2+y2)k.(16)

Dividing by (x+y)δ(n), we get (17) xn+yn(x+y)δ(n)=k=0n2fk(n)(2xy+x2+y2)n2k(x2+y2)k.(17)

Hence, we get the integer sequence Ψk(n) that satisfy (18) xn+yn(x+y)δ(n)=k=0n2Ψk(n)(xy)n2k(x2+y2)k.(18)

5. The Eight Levels theorem

Now we prove what we call the Eight Levels Theorem for the expansion of the polynomial xn+yn(x+y)δ(n) in terms of the symmetric polynomials xy and x2+y2. Then we investigate some properties for the coefficients of the expansion (19) xn+yn(x+y)δ(n)=k=0n2Ψk(n)(xy)n2k(x2+y2)k.(19)

Then we study some applications and prove the results of the summary one by one.

5.1. The statement of the Eight Levels theorem

Theorem 9.

(The Eight Levels Theorem)

For any complex numbers x, y, any non negative integers n, k, the coefficients Ψk(n) of the expansion (20) xn+yn(x+y)δ(n)=k=0n2Ψk(n)(xy)n2k(x2+y2)k(20) are integers and (21) Ψ0(n)={+2n0(mod 8)+1n±1(mod 8)0n±2(mod 8)1n±3(mod 8)2n±4(mod 8)(21) and, for each 1kn2, the coefficients satisfy the following statements

  • For n0,2,4,6(mod 8):

    n0(mod 8)¯ Ψk(n)={0for k odd2(1)k2λ=0k21[n2-(4λ)2]4kk!for k even

n2(mod 8)¯ Ψk(n)={0for k even2(1)k2nλ=1k2[n2-(4λ2)2]4kk!for k odd

n4(mod 8)¯ Ψk(n)={0for k odd2(1)k2+1λ=0k21[n2-(4λ)2]4kk!for k even

n6(mod 8)¯ Ψk(n)={0for k even2(1)k2+1nλ=1k2[n2-(4λ2)2]4kk!for k odd

  • For n1,3,5,7(mod 8):

    n1(mod 8)¯

Ψk(n)=(1)k2(n+12k)δ(k)λ=1k2[(n+1)2-(4λ2)2]4kk!

n3(mod 8)¯ Ψk(n)=(1)k2+δ(k1)(n+1)(n+12k)δ(k1)λ=1k12[(n+1)2(4λ)2]4kk!

n5(mod 8)¯ Ψk(n)=(1)k2+1(n+12k)δ(k)λ=1k2[(n+1)2-(4λ2)2]4kk!

n7(mod 8)¯ Ψk(n)=(1)k2+δ(k)(n+1)(n+12k)δ(k1)λ=1k12[(n+1)2(4λ)2]4kk!  

Proof.

Put (x,y)=(1,1) in Equation(20), we immediately get Equation(21). To prove the statements of Theorem Equation(9), we need to prove the following lemmas.

Lemma 10.

For each natural number n, 1kn2, the coefficients Ψk(n) of the expansion of Equation(20) are integers, unique and satisfy (22) Ψk(n)=Ψk1(n2)Ψk(n4).(22)

Proof.

From the fundamental theorem on symmetric polynomials (Ibrahim, Citation2004; Perron, Citation1932), we have a sequence of integers Ψk(n) satisfy the expansion of Equation(20). From the algebraic independence of xy,x2+y2, the coefficients Ψk(n) of the expansion of Equation(20) are unique. Now, multiply Equation(20) by xy(x2+y2), and noting (x+y)δ(n+2)=(x+y)δ(n)=(x+y)δ(n2) and n+22=n2+1 and n22=n21, we get (23) k=1n2+1Ψk1(n)(xy)n2+2k(x2+y2)k=k=0n2+1Ψk(n+2)(xy)n2+2k(x2+y2)k+k=0n21Ψk(n2)(xy)n2+2k(x2+y2)k(23)

Again from the algebraic independence of xy,x2+y2, and from Equation(23), we get the following identity for any natural number n: (24) Ψk(n+2)=Ψk1(n)Ψk(n2).(24)

Replace n by n − 2 in Equation(24), we get Equation(22). □

Lemma 11.

For every even natural number n, 0kn2, the following statements are true for each case: (25) n0or4(mod 8)¯Ψk(n)=0forkodd,n2 or 6(mod 8)¯Ψk(n)=0forkeven.(25)

Proof.

Consider n even natural number, and replace δ(n)=0 in Equation(20) to get (26) xn+yn=k=0n2Ψk(n)(xy)n2k(x2+y2)k(26)

Then replace x by—x in Equation(26), we get (27) xn+yn=k=0n2Ψk(n)(xy)n2k(x2+y2)k(27)

From the algebraic independence of xy,x2+y2, and from (26, 27) we get the proof. □

From Equation(20), we get the following initial values for Ψk(n).

Lemma 12.

Ψ0(0)=+2,Ψ1(0)=0,Ψ2(0)=0,Ψ3(0)=0,Ψ0(2)=0,Ψ1(2)=+1,Ψ2(2)=0,Ψ3(2)=0,Ψ0(4)=2,Ψ1(4)=0,Ψ2(4)=+1,Ψ3(4)=0,Ψ0(6)=0,Ψ1(6)=3,Ψ2(6)=0,Ψ3(6)=+1.  

Lemma 13.

For any odd natural number n, the coefficients Ψk(n) of the expansion of Equation(20) satisfy the following property (28) Ψk(n)=2(k+1)(n+1)Ψk+1(n+1)+(n+12k)(n+1)Ψk(n+1).(28)

Proof.

For n odd, n + 1 is even. Then from the expansion of Equation(20) we get the following (29) xn+1+yn+1=k=0n+12Ψk(n+1)(xy)n+12k(x2+y2)k.(29)

Now acting the differential operator (x+y) on Equation(29) and noting that (x+y)(xn+1+yn+1)=(n+1)(xn+yn),(x+y)xy=x+y,(x+y)(x2+y2)=2(x+y), and equating the coefficients, we get the proof. □

5.2. The proof of Theorem (9) for n0,2,4,6(mod 8)

Now in this section we prove Theorem Equation(9) for n even. The values of Ψk(n) that comes from the formulas of Theorem Equation(9) for n=0,,2,4,8 are identical with the correct values that come from Lemma Equation(12). So, Theorem Equation(9) is correct for the initial vales n=0,2,4,6. Now we assume the validity of Theorem Equation(9) for each 0m<n, with m0,2,4,6(mod8) and need to prove that Theorem Equation(9) for n. Lemma Equation(11) proves the validity of Theorem Equation(9) for n0 or 4(mod 8) if k odd, and for n2 or 6(mod 8) if k even. Therefore it remains to prove the validity of the following cases:

  • n0(mod 8),keven

  • n2(mod 8),kodd

  • n4(mod 8),keven

  • n6(mod 8),kodd

We prove these cases one by one as following.

Consider n0(mod 8)andkeven

Proof.

In this case, n26(mod 8) and n44(mod 8) and from Lemma Equation(10), we get Ψk(n)=Ψk1(n2)Ψk(n4)=2(1)k12+1(n2)λ=1k12(n2)2(4λ2)24k1(k1)!2(1)k2+1λ=0k21(n4)2(4λ)24kk!=2(1)k2(n2)λ=1k21(n4λ)λ=1k21(n+4λ4)4k1(k1)!+2(1)k2λ=0k21(n4λ4)λ=0k21(n+4λ4)4k1(k1)!=2(1)k24kk!λ=1k21(n4λ)λ=0k22(n+4λ)[4k(n2)+(n4k2)(n4)]=2(1)k24kk!λ=1k21(n4λ)λ=0k22(n+4λ)n(n+4(k21))=2(1)k24kk!λ=0k21(n4λ)λ=0k21(n+4λ)=2(1)k2λ=0k21n2(4λ)24kk!.  

Consider n2(mod 8)andkodd

Proof.

In this case, n20(mod 8) and n46(mod 8) and from Lemma Equation(10), we get Ψk(n)=Ψk1(n2)Ψk(n4)=2(1)k12λ=0k121(n2)2(4λ)24k1(k1)!2(1)k2+1(n4)λ=1k12(n4)2(4λ2)24kk!=2(1)k2λ=0k21(n4λ2)λ=0k21(n+4λ2)4k1(k1)!+2(1)k2(n4)λ=1k2(n4λ2)λ=1k2(n+4λ6)4kk!=2(1)k24kk!λ=1k21(n4λ2)λ=0k21(n+4λ2)[4k(n2)+(n4)(n4(k2)2)].  

As 4k(n2)+(n4)(n4(k2)2)=n(n+2k4), we get the following Ψk(n)=2(1)k24kk!nλ=1k2(n4λ+2)λ=1k2(n+4λ2)=2n(1)k2λ=1k2n2(4λ2)24kk!.  

Consider n4(mod 8)andkeven

Proof.

In this case, n22(mod 8) and n40(mod 8) and from Lemma Equation(10), we get Ψk(n)=Ψk1(n2)Ψk(n4)=2(1)k12(n2)λ=1k22(n2)2(4λ2)24k1(k1)!2(1)k2λ=0k21(n4)2(4λ)24kk!=2(1)k2+1(n2)λ=1k21(n4λ)λ=1k21(n+4λ4)4k1(k1)!+2(1)k2+1λ=0k21(n4λ4)λ=0k21(n+4λ4)4kk!  

Hence Ψk(n)=2(1)k2+14kk!λ=1k21(n4λ)λ=0k22(n+4λ)n(n+4(k21))=2(1)k2+14kk!λ=0k21(n4λ)λ=0k21(n+4λ)=2(1)k2+1λ=0k21n2(4λ)24kk!.  

Consider n6(mod 8)andkodd

Proof.

In this case, n24(mod 8) and n42(mod 8) and from Lemma Equation(10), we get Ψk(n)=Ψk1(n2)Ψk(n4)=2(1)k12+1λ=0k121(n2)2(4λ)24k1(k1)!2(1)k2(n4)λ=1k12(n4)2(4λ2)24kk!=2(1)k2+1λ=0k21(n4λ2)λ=0k21(n+4λ2)4k1(k1)!+2(1)k2+1(n4)λ=1k2(n4λ2)λ=1k2(n+4λ6)4kk!  

Hence Ψk(n)=2(1)k2+14kk!nλ=1k2(n4λ+2)λ=1k2(n+4λ2)=2n(1)k2+1λ=1k2n2(4λ2)24kk!.  

5.3. The proof of Theorem (9) for n1,3,5,7(mod 8)

With the help of Lemma Equation(13) and Theorem Equation(9) for the even case that we already proved, together with Lemma Equation(11), we prove Theorem Equation(9) for the odd cases n1,3,5,7(mod 8), one by one, for each parity for k.

Consider n1(mod 8)andkodd

Proof.

In this case, n+12(mod 8) and from Lemma Equation(11), we get Ψk+1(n+1)=0. Hence, from Lemma Equation(13), and from Theorem Equation(9), we get the following relation Ψk(n)=(n+12k)(n+1)Ψk(n+1)=(1)k2(n+12k)λ=1k2(n+1)2(4λ2)24kk!.  

Consider n1(mod 8)andkeven

Proof.

In this case, n+12(mod 8) and from Lemma Equation(11), we get Ψk(n+1)=0. From Lemma Equation(13), and from Theorem Equation(9), we get the following relation Ψk(n)=2(k+1)(n+1)Ψk+1(n+1)=(1)k2λ=1k2(n+1)2(4λ2)24kk!.  

Consider n3(mod 8)andkodd

Proof.

In this case, from Lemma Equation(11), we get Ψk(n+1)=0. Hence, from Lemma Equation(13), and from Theorem Equation(9), we get the following relation Ψk(n)=2(k+1)(n+1)Ψk+1(n+1)=(1)k2(n+1)λ=1k2(n+1)2(4λ)24kk!.  

Consider n3(mod 8)andkeven

Proof.

In this case, Ψk+1(n+1)=0. From Lemma Equation(13), and from Theorem Equation(9), we get the following relation Ψk(n)=(n+12k)(n+1)Ψk(n+1)=(1)k2+1(n+1)(n+12k)λ=1k21(n+1)2(4λ)24kk!.  

Consider n5(mod 8)andkodd

Proof.

In this case, from Lemma Equation(11), we get Ψk+1(n+1)=0. Hence, from Lemma Equation(13), and from Theorem Equation(9), we get the following relation Ψk(n)=(n+12k)(n+1)Ψk(n+1)=(1)k2+1(n+12k)λ=1k2(n+1)2(4λ2)24kk!.  

Consider n5(mod 8)andkeven

Proof.

In this case, from Lemma Equation(11), we get Ψk(n+1)=0. From Lemma Equation(13), and from Theorem Equation(9), we get the following relation Ψk(n)=2(k+1)(n+1)Ψk+1(n+1)=(1)k2+1λ=1k2(n+1)2(4λ2)24kk!.  

Consider n7(mod 8)andkodd

Proof.

In this case, from Lemma Equation(11), we get Ψk(n+1)=0. Hence, from Lemma Equation(13), and from Theorem Equation(9), we get the following relation Ψk(n)=2(k+1)(n+1)Ψk+1(n+1)=(1)k2+1(n+1)λ=1k2(n+1)2(4λ)24kk!.  

Consider n7(mod 8)andkeven

Proof.

In this case, Ψk+1(n+1)=0. From Lemma Equation(13), and from Theorem Equation(9), we get the following relation Ψk(n)=(n+12k)(n+1)Ψk(n+1)=(1)k2(n+1)(n+12k)λ=1k21(n+1)2(4λ)24kk!.  

This completes the proof of Theorem Equation(9). □

6. General characteristics for Ψ sequence

6.1. Examples for Ψ sequence

From Theorem Equation(9), we list some examples that show the splendour of the natural factorization of Ψk(n) for k=0,1,2,3,4,5,6,7: (30) Ψ0(n)={+2n0(mod 8)+1n1(mod 8)0n2(mod 8)1n3(mod 8)2n4(mod 8)1n5(mod 8)0n6(mod 8)+1n7(mod 8),Ψ1(n)={0n0(mod 8)+(n1)4n1(mod 8)+2n4n2(mod 8)+(n+1)4n3(mod 8)0n4(mod 8)(n1)4n5(mod 8)2n4n6(mod 8)(n+1)4n7(mod 8),  (30) (31) Ψ2(n)={2n2422!n0(mod 8)(n1)(n+3)422!n1(mod 8)0n2(mod 8)+(n3)(n+1)422!n3(mod 8)+2n2422!n4(mod 8)+(n1)(n+3)422!n5(mod 8)0n6(mod 8)(n3)(n+1)422!n7(mod 8),Ψ3(n)={0n0(mod 8)(n5)(n1)(n+3)433!n1(mod 8)2(n2)n(n+2)433!n2(mod 8)(n3)(n+1)(n+5)433!n3(mod 8)0n4(mod 8)+(n5)(n1)(n+3)433!n5(mod 8)+2(n2)n(n+2)433!n6(mod 8)+(n3)(n+1)(n+5)433!n7(mod 8),  (31) (32) Ψ4(n)={+2(n4)n2(n+4)444!n0(mod 8)+(n5)(n1)(n+3)(n+7)444!n1(mod 8)0n2(mod 8)(n7)(n3)(n+1)(n+5)444!n3(mod 8)2(n4)n2(n+4)444!n4(mod 8)(n5)(n1)(n+3)(n+7)444!n5(mod 8)0n6(mod 8)+(n7)(n3)(n+1)(n+5)444!n7(mod 8),(32) and (33) Ψ5(n)={0n0(mod 8)+(n9)(n5)(n1)(n+3)(n+7)455!n1(mod 8)+2(n6)(n2)n(n+2)(n+6)455!n2(mod 8)+(n7)(n3)(n+1)(n+5)(n+9)455!n3(mod 8)0n4(mod 8)(n9)(n5)(n1)(n+3)(n+7)455!n5(mod 8)2(n6)(n2)n(n+2)(n+6)455!n6(mod 8)(n7)(n3)(n+1)(n+5)(n+9)455!n7(mod 8),(33) and (34) Ψ6(n)={2(n8)(n4)n2(n+4)(n+8)466!n0(mod 8)(n9)(n5)(n1)(n+3)(n+7)(n+11)466!n1(mod 8)0n2(mod 8)+(n11)(n7)(n3)(n+1)(n+5)(n+9)466!n3(mod 8)+2(n8)(n4)n2(n+4)(n+8)466!n4(mod 8)+(n9)(n5)(n1)(n+3)(n+7)(n+11)466!n5(mod 8)0n6(mod 8)(n11)(n7)(n3)(n+1)(n+5)(n+9)466!n7(mod 8),,(34) and the following example should give us a better vision about that sequence: (35) Ψ7(n)={0n0(mod 8)(n13)(n9)(n5)(n1)(n+3)(n+7)(n+11)477!n1(mod 8)2(n10)(n6)(n2)n(n+2)(n+6)(n+10)477!n2(mod 8)(n11)(n7)(n3)(n+1)(n+5)(n+9)(n+13)477!n3(mod 8)0n4(mod 8)+(n13)(n9)(n5)(n1)(n+3)(n+7)(n+11)477!n5(mod 8)+2(n10)(n6)(n2)n(n+2)(n+6)(n+10)477!n6(mod 8)+(n11)(n7)(n3)(n+1)(n+5)(n+9)(n+13)477!n7(mod 8).  (35)

As we see, sometimes the numerator of the sequence Ψk(n) has a centre factor that all the other factors in the numerators get around it. For n0,2,4,6(mod 8), it is symmetric in this sense; meaning if (na) is a factor of the numerator then (n+a) would be a factor of that numerator and vice versa. However, it gives a different story for n1,3,5,7(mod 8), and a natural phenomenon for these numbers arises up here and needs a closer attention.

6.2. The right tendency concept for Ψk(n)

For n±1(mod 8), from the data above, and from the formulas of Ψk(n), we can observe that the factors of the numerators for k=1,2,3,4,5,6, always fill the right part first then go around the centre factor to fill the left part and so on. For example, for k = 1 the centre factor is (n1). Then for k = 2, the next factor (n+3) is located on the right of (n1). We call this natural behaviour by the right tendency which is explained in table Equation(36). (36) Right Tendency For Ψk(n)kThe Center1(n1)2(n1) (n + 3)3(n5)(n1)(n+3)4(n5)(n1)(n+3)(n+7)5(n9)(n5)(n1)(n+3)(n+7)6(n9)(n5)(n1)(n+3)(n+7)(n+11)(36)

6.3. The left tendency concept for Ψk(n)

However, for n±3(mod 8), and from the data above, and the formulas of Ψk(n), we can also observe that the factors of the numerators for k=1,2,3,4,5,6, always fill the left part first then go around the centre factor to fill the right part and so on. For example, for k = 1 the centre factor is (n+1). Then for k = 2, the next factor (n3) is located on the left of (n+1). We call this natural behaviour by the left tendency which is explained in table Equation(37). (37) Left Tendency For Ψk(n)kThe Center1(n+1)2 (n-3)(n+1)3(n3)(n+1)(n+5)4(n7)(n3)(n+1)(n+5)5(n7)(n3)(n+1)(n+5)(n+9)6(n11)(n7)(n3)(n+1)(n+5)(n+9)(37)

6.4. The signs phenomena of Ψ integers

For k0,1,2,3,4,5,6,7(mod 8), we get the following data (38) k(mod 8)(1)k2(1)k2+1(1)k2+δ(k)(1)k2+δ(k1)0+11+111+111+121+11+131+1+114+11+115+111+161+11+171+1+11(38)

Therefore, the signs of the Ψ sequence get back periodically every 8 steps. Table Equation(39) shows that three plus are followed by zero then it followed by three minus then followed by zero and so on. For n0,1,2,3,4,5,6,7(mod 8) and for k0,1,2,3,4,5,6,7(mod 8), the following data is useful (39) Signs of Ψk(n)n(mod 8)Ψ0(n)Ψ1(n)Ψ2(n)Ψ3(n)Ψ4(n)Ψ5(n)Ψ6(n)Ψ7(n)0+00+001++++20+00+03++++40+00+05++++600+00+7++++(39)

6.5. General divisibility relationship for the Ψ integers

Whether n0,1,2,3,4,5,6,7(mod 8), it is rather surprising that the quantity Ψk(n)÷Ψk2(n) always gives n2(2k4)216k(k1)or(n±1)2(2k2)216k(k1)or0÷0.

Therefore, computing these ratios, at the eight levels n0,1,2,3,4,5,6,7(mod 8), immediately proves the following desirable theorem that gives a fraction with a difference of two squares divided by a multiplication of the factors 16, k, and k − 1.

Theorem 14.

(Generating The Ψintegers From The Previous Term)

If Ψk(n) not identically zero, then Ψk(n)Ψk2(n)=[n+(1)n2+kδ(n)]2[2k22δ(n1)]216k(k1), and we can choose either of the following initial values to generate Ψk(n) from the starting term Ψ0(n) or the last term Ψn2(n): Ψ0(n)={+2n0(mod 8)+1n±1(mod 8)0n±2(mod 8)1n±3(mod 8)2n±4(mod 8),Ψn2(n)=1.  

7. The emergence of ϕ sequence

Another natural sequence that emerges naturally from Ψk(n) is the integer sequence ϕk(n) which is defined as follows

Definition 15.

ϕk(n):=4kΨk(n).

7.1. Recurrence relation of order 4 to generate ϕ sequence

From Equation(10), we get the following recurrence relation

Lemma 16.

For each natural number n, 0kn2, the integers ϕk(n) satisfy the following property (40) ϕk(n)=4ϕk1(n2)ϕk(n4).(40)

From Equation(30), we easily get the following initial values

Lemma 17.

For each natural number n (41) ϕ0(n)={+2n0(mod 8)+1n±1(mod 8)0n±2(mod 8)1n±3(mod 8)2n±4(mod 8),ϕ1(n)={0n0(mod 8)+(n1)n1(mod 8)+2nn2(mod 8)+(n+1)n3(mod 8)0n4(mod 8)(n1)n5(mod 8)2nn6(mod 8)(n+1)n7(mod 8).  (41)

7.2. Explicit formulas for ϕ sequence

Now, from Theorem Equation(9), we get the following explicit formulas for the integer sequence ϕk(n).

Lemma 18.

For any non negative integers n, k, the sequences ϕk(n) satisfy the following statements

n0(mod 8)¯ ϕk(n)={0for k odd2(1)k2λ=0k21[n2(4λ)2]k!for k even    n1(mod 8)¯ ϕk(n)=(1)k2(n+12k)δ(k)λ=1k2[(n+1)2(4λ2)2]k!   n2(mod 8)¯ ϕk(n)={0for k even2(1)k2nλ=1k2[n2(4λ2)2]k!for k odd n3(mod 8)¯ ϕk(n)=(1)k2+δ(k1)(n+1)(n+12k)δ(k1)λ=1k2δ(k1)[(n+1)2(4λ)2]k!   n4(mod 8)¯ ϕk(n)={0for k odd2(1)k2+1λ=0k21[n2(4λ)2]k!for k even    n5(mod 8)¯ ϕk(n)=(1)k2+1(n+12k)δ(k)λ=1k2[(n+1)2(4λ2)2]k!   n6(mod 8)¯ ϕk(n)={0for k even2(1)k2+1nλ=1k2[n2(4λ2)2]k!for k odd n7(mod 8)¯ ϕk(n)=(1)k2+δ(k)(n+1)(n+12k)δ(k1)λ=1k2δ(k1)[(n+1)2(4λ)2]k!

7.3. Nonlinear recurrence relation to generate ϕ sequence

To study the arithmetic properties of ϕ Sequence, we need to generate ϕk(n) from the previous one, ϕk2(n), or generate ϕk(n) from the next one, ϕk+2(n). As Ψk(n)Ψk2(n)=ϕk(n)42ϕk2(n), and noting that δ(n)=0,δ(n1)=1, for n even, we immediately get, from Theorem Equation(14), the following desirable theorem.

Theorem 19.

(Generating The ϕintegers From The Previous Term)

If ϕk(n) not identically zero, n even, then ϕk(n)ϕk2(n)=n2(2k4)2k(k1), and depending on the parity of k we can choose either of the following initial values to generate ϕk(n) from the starting terms ϕ0(n) or ϕ1(n): (42) ϕ0(n)={+2n0(mod 8)2n4(mod 8),ϕ1(n)={+2 nn2(mod 8)2 nn6(mod 8). (42)

8. Three new versions for Lucas-Lehmer primality tests for Mersenne numbers

From Lemmas (16–18), we are ready now to prove the following new results.

8.1. The proof of the first new version

Theorem 20.

Given prime p5. 2p1 is prime if and only if (43) 2p1|k evenk=0,n2ϕk(n)(43) where n:=2p1,ϕk(n) are defined by the double index recurrence relation ϕk(m)=4 ϕk1(m2)ϕk(m4) and the initial boundary values satisfy (44) ϕ0(m)={+2m0(mod 8)0m±2(mod 8)2m4(mod 8),ϕ1(m)={+2mm2(mod 8)0m0,4(mod 8)2mm6(mod 8).  (44)

Proof.

Given prime p5, let n:=2p1. From Lucas-Lehmer-Test, see Gallier(Citation2011) and Deza (Citation2021), we have 2p1is prime2p1|(1+3)n+(13)n.

Hence, as n0(mod 8), replace x=1+3,y=13 in Theorem Equation(9), we get the following equivalent statements: 2p1is prime2p1|kevenk=0,n2Ψk(n) (2)n2k(8)k2p1|2n2kevenk=0,n2Ψk(n)4k2p1|kevenk=0,n2ϕk(n).

Lemmas (16, 17) already proved the rest of Theorem Equation(20). □

8.2. The proof of the second new version

From Lemmas (16, 17), we should observe that the recurrence relation of ϕk(n) is always even integer for n0(mod 8). Hence from Lemma Equation(18), we get the following theorem.

Theorem 21.

Given prime p5,n:=2p1. The number 2p1 is prime if and only if (45) 2 n1|kevenk=0,n2(1)k2λ=0k21n2(4λ)2k!.(45)

8.3. The proof of the third new version

When we generate the ϕ integers needed for Mersenne numbers, we should notice that for n:=2p1,p5, we have n0(mod 8).

From the proof of Theorem Equation(20), we know that 2p1is prime2p1|kevenk=0,n2ϕk(n).

It is plain that, from EquationEq. (20), Ψn2(n)=1. Consequently, ϕn2(n)=4n2=2n. Therefore, from Theorem Equation(19), we get the following result.

Theorem 22.

Given prime p5, n:=2p1. The number 2p1 is prime if and only if (46) 2n1|kevenk=0,n2ϕk(n)(46) where ϕk(n) are defined and generated by the double index recurrence relation ϕk(n)ϕk2(n)=(2k4)2n2k(k1), and we can choose either of the following initial values to generate ϕk(n) from the starting term ϕ0(n) or the last term ϕn2(n): ϕ0(n)=+2,ϕn2(n)=2n. 

9. Criteria for compositeness of Mersenne numbers

The following theorem is an immediate consequence of Theorem Equation(21).

Theorem 23.

(Criteria for compositeness of Mersenne numbers)

Given prime p5. The number 2n1=2p1 is Mersenne composite number if (47) 2n-1kevenk=0,n2λ=0k21[(4λ)2-41]k!.(47)

10. The proof of the combinatorial identities

Choosing k=n2 in the Eight Levels Theorem Equation(9), and noting Ψn2(n)=1, we surprisingly get the following eight combinatorial identities which reflect some unexpected facts about the nature of numbers. n0(mod 8)¯4n2(n2)!=2λ=0n41n2(4λ)2n1(mod 8)¯4n2(n2)!=λ=1n14(n+1)2(4λ2)2n2(mod 8)¯4n2(n2)!=2nλ=1n24n2(4λ2)2n3(mod 8)¯4n2(n2)!=(n+1)λ=1n34(n+1)2(4λ)2n4(mod 8)¯4n2(n2)!=2λ=0n44n2(4λ)2n5(mod 8)¯4n2(n2)!=λ=1n14(n+1)2(4λ2)2n6(mod 8)¯4n2(n2)!=2nλ=1n24n2(4λ2)2n7(mod 8)¯4n2(n2)!=(n+1)λ=1n34(n+1)2(4λ)2

Writing only the different identities, which are 4, we get the following theorem which gives formulas to factor any factorial in terms of a product of difference of squares.

Theorem 24.

(combinatorial identities) For any natural number n, the following combinatorial identities are correct n0(mod 4)¯4n2(n2)!=2λ=0n44[n2-(4λ)2]n1(mod 4)¯4n2(n2)!=λ=1n14[(n+1)2-(4λ2)2]n2(mod 4)¯4n2(n2)!=2nλ=1n24[n2-(4λ2)2]n3(mod 4)¯4n2(n2)!=(n+1)λ=1n34[(n+1)2-(4λ)2]

11. Formulas

Now we compute the summation (48) kevenk=0,n2ϕk(n)(48) for the first few terms from the starting and from the end. From Theorem Equation(22), we get the following explicit formulas: ϕ0(n)=+2.

Now, compute ϕ2(n) as following ϕ2(n)=n2(44)22(21)ϕ0(n)=n2.

Proceeding this way, we get ϕ4(n)=n2(84)24(41)ϕ2(n)=+n2(n242)12.

Similarly ϕ6(n)=n2(n242)(n282)360, ϕ8(n)=+n2(n242)(n282)(n2122)20160,etc

Now, compute ϕk(n) from the end, and from Theorem Equation(22), we get ϕk2(n)=k(k1)n2(2k4)2ϕk(n).

Initially ϕn2(n)=2n.

Hence ϕn22(n)=n2n5, and ϕn24(n)=+n(n6)2n11.

Similarly ϕn26(n)=n(n8)(n10)32n16, ϕn28(n)=+n(n10)(n12)(n14)32n23,etc.

The author feels that we need a clever way to evaluate the sum Equation(48). We may like to add the terms in a way reflect some elegant arithmetic. Remember that we do not need to compute the sum Equation(48) exactly; but we just need to find the sum modulo 2n1. According to the following theorem, and working modulo 2n1, the last term always gives the value of the first term.

Theorem 25.

For p5 prime, and n:=2p1, (49) ϕn2(n)2(mod 2n1).(49)

Proof.

We should observe that if p5 prime, then n=2p1=1+ζp, for some positive integer ζ. Then ϕn2(n)=2n=21+ζp=21(2p)ζ2(mod 2n1).

Hence, this encourages one to compute the partial sums of (50) kevenk=0,n2ϕk(n)(mod 2n1),(50) in the following order ϕn2(n)+ϕ0(n)+ϕ2(n)+ϕ4(n)++ϕn21(n).

12. The 5 scenario

For example, take p = 5, then n=24. Hence (51) ϕ8(16)+2(mod 31),ϕ0(16)+2(mod 31),ϕ2(16)23(mod 31),ϕ4(16)+22+1(mod 31),ϕ6(16)1(mod 31).(51)

Hence, we get the partial sums (52) ϕ8(16)+2,ϕ8(16)+ϕ0(16)+22,ϕ8(16)+ϕ0(16)+ϕ2(16)22,ϕ8(16)+ϕ0(16)+ϕ2(16)+ϕ4(16)+1,ϕ8(16)+ϕ0(16)+ϕ2(16)+ϕ4(16)+ϕ6(16)0.(52)

As we ended up with zero, this shows that 251=31 is Mersenne prime.

13. Further research investigations

Now, consider p prime greater than 3, and n:=2p1. The previous sections give various explicit formulas and techniques to compute and generate all of the terms ϕk(n) needed for checking the primality of the Mersenne number 2p1.

13.1. Searching for a hypothetical pattern

This previous particular example, for p=5=22+1, is illuminating and should motivate us for further theoretical investigations for other similar scenarios. We need to find a general pattern for n similar to this special case, for p = 5, such that, working modulo 2n1, the partial sums give: (53) ϕn2(n)+2,ϕn2(n)+ϕ0(n)+22,ϕn2(n)+ϕ0(n)+ϕ2(n)22+ϵ1,ϕn2(n)+ϕ0(n)+ϕ2(n)+ϕ4(n)+23+ϵ2,ϕn2(n)+ϕ0(n)+ϕ2(n)+ϕ4(n)+ϕ6(n)23+ϵ3,ϕn2(n)+ϕ0(n)+ϕ2(n)+ϕ4(n)+ϕ6(n)+ϕ8(n)+24+ϵ4,ϕn2(n)+ϕ0(n)+ϕ2(n)+ϕ4(n)+ϕ6(n)+ϕ8(n)+ϕ10(n)24+ϵ5,(53) where ϵk{0,1,+1} for each k. We need to investigate certain types of prime p. We should try to investigate the primes p of the form 2a±1 or 2a±2b±1 to make the sum (54) kevenk=0,n2ϕk(n)(mod 2n1),(54) easily to understand and help us identify a general pattern for p for which this sum gives zero modulo 2n1 infinitely many times (in this case we could prove that Mersenne primes are infinite), or gives nonzero modulo 2n1 infinitely many times (in this case we could prove that Mersenne composites are infinite).

13.2. New identities to help understand the sums of (48) and (54)

To motivate the readers about this future vital research and investigations, and in the spirit of the previous results, I feel compelled to mention, even though succinctly, the following new identities that I found recently for ϕk(n):

For any natural number p greater than 3, and n:=2p1, we get the identities (55)  kevenk=0,n2ϕk(n) 2k=2, kevenk=0,n2ϕk(n) 22k=1, kevenk=0,n2ϕk(n) 22k 3k=L(n),(55) where L(n) is Lucas sequence defined by L(m+1)=L(m)+L(m1),L(1)=1,L(0)=2.

13.3. The nature of the prime factors of the sum

We use SAGE (SAGE Mathematical Software, Version 2.6), with double-checking provided by Mathematica, to carry out all the numerical computations of the current paper. For p=5, n=16, the following sum (56) k evenk=0,n2ϕk(n)= 2×31×607.(56)

This shows that the number 31 is Mersenne prime because it is one of the factors of this sum. Surprisingly, although the prime number 607 is not a Mersenne prime, the number 607 is the exponent of the Mersenne prime 26071. Moreover, the prime 607 is an irregular prime since it divides the numerator of the Bernoulli number B592. So, for a given prime p, n:=2p1, we should also investigate the arithmetical nature of all the other prime factors for the sum (57) k evenk=0,n2ϕk(n).(57)

14. Conclusions

The discovery of new Mersenne primes is significant in the field of mathematics because they are relatively rare and difficult to find. Moreover, they have important applications in areas such as cryptography, number theory, and computer science. As of April 2023, there are currently only 51 known Mersenne primes. The largest known Mersenne prime as of this date is 282,589,9331, which has 24,862,048 digits. The discovery of Mersenne primes is a very challenging, computationally intensive process, and time-consuming task, and the search for new Mersenne primes via GIMPS is an ongoing effort by many mathematicians and computer scientists around the world. While the Great Internet Mersenne Prime Search (GIMPS) has been successful in discovering 51 Mersenne primes to date, there are several limitations to the search. The search for Mersenne primes is a probabilistic process, and it is not guaranteed that a new Mersenne prime will be found. Even if a new Mersenne prime exists, there is a chance that it may not be discovered in a given search due to the limitations of the algorithm or computational resources. The search becomes increasingly difficult as the size of the potential primes increases, and it may become impractical or impossible to search for larger Mersenne primes without significant advances in theoretical understanding of the nature of Mersenne primes. Therefore, the current paper developed three new versions of primality tests for Mersenne numbers which could potentially help in the search for new Mersenne primes or provide insights into the open questions surrounding the infinitude of Mersenne primes or Mersenne composites.

Disclosure statement

The author declares that he has no conflict of interest.

References

  • Aggarwal, D., Joux, A., Prakash, A., & Santha, M. (2018). A new public-key cryptosystem via Mersenne numbers. In CRYPTO (Vol. 10993, pp. 459–482). Cham: Springer.
  • Cambraia, A., Knapp, P., Lemos, A., Moriya, B. K., & Rodrigues, P. H. A. (2022). On prime factors of Mersenne numbers. Palestine Journal of Mathematics, 11(2), 449–456.
  • Deza, E. (2021). Mersenne numbers and Fermat numbers. Hackensack, NJ: World Scientific Publishing Co. Pte. Ltd.
  • Dickson, L. E. (1919). History of the theory of numbers (Vol. I). Washington: Carnegie Institution of Washington.
  • Effinger, G., & Mullen, G. L. (2021). Elementary number theory. New York: Chapman and Hall/CRC.
  • Ferradi, H., & Xagawa, K. (2020). Post-quantum provably-secure authentication and MAC from Mersenne primes. In Topics in cryptology–CT-RSA 2020 (Lecture notes in computer science, Vol. 12006, pp. 469–495). Cham: Springer.
  • Gallier, J. (2011). Discrete mathematics. New York: Universitext, Springer-Verlag.
  • GIMPS. (2016). Great internet mersenne prime search, project [Online]. Retrieved September 16, from http://www.mersenne.org/default.php
  • Guy, R. (2004). Unsolved problems in number theory. New York: Springer.
  • Hurwitz, A. (1962). New Mersenne primes. Mathematics of Computation, 16(78), 249–251. doi:10.1090/S0025-5718-1962-0146162-X
  • Ibrahim, M. (2004). A new approach to polynomial identities. The Ramanujan Journal, 8, 423–457.
  • Jo, G., & Kim, D. (2022). Mersenne prime factor and sum of binomial coefficients. Journal of Applied Mathematics and Informatics, 40(1–2), 61–68.
  • Kharaghani, H., & Suda, S. (2020). Commutative association schemes obtained from twin prime powers, Fermat primes, Mersenne primes. Finite Fields and Their Applications, 63, 101631. doi:10.1016/j.ffa.2019.101631
  • Kraft, J. S., & Washington, L. C. (2018). An introduction to number theory with cryptography (2nd ed.). New York: Chapman and Hall/CRC.
  • Kundu, S., & Mazumder, S. (2022). Number Theory and its Applications. London: CRC Press.
  • Ochem, P., & Rao, M. (2012). Odd perfect numbers are greater than 101500. Mathematics of Computation, 81(279), 1869–1877. doi:10.1090/S0025-5718-2012-02563-4
  • Perron, O. (1932). Algebra I, Die Grundlagen. Berlin: Walter de Gruyter.
  • SAGE Mathematical Software, Version 2.6. http://www.sagemath.org.
  • Skula, L. (2019). Prime power divisors of Mersenne numbers and Wieferich primes of higher order. Integers, 19, 4.
  • Wells, D. (1987). The penguin dictionary of curious and interesting numbers. New York: Penguin Books.
  • Witno, A. (2021). Hypothetical elite primes for Mersenne numbers and repunits. Journal of Integer Sequences, 24(1), 21.1.7.
  • Youssri, Y. H., Abd-Elhameed, W. M., & Atta, A. G. (2022). Spectral Galerkin treatment of linear one-dimensional telegraph type problem via the generalized Lucas polynomials. Arabian Journal of Mathematics, 11(3), 601–615. doi:10.1007/s40065-022-00374-0
  • Zheng, D. X., & Yang, Q. (1985). On the complex Mersenne transforms. Sichuan Daxue Xuebao, 1, 6–9.