Improved Riemann Ladder Function
Dan Liu
Sichuan Normal University, Chengdu, China.
DOI: 10.4236/jamp.2021.94042   PDF    HTML   XML   233 Downloads   638 Views   Citations

Abstract

In Riemann’s prime distribution formula, there is a key Riemann ladder function. This is the basis for Riemann to study the distribution of prime numbers. But the calculation of Riemann’s ladder function is very complicated. According to the definition of Riemannian ladder function, we greatly improve the Riemannian ladder function, obtain new ladder function and improved Riemannian prime distribution formula, and prove the improved Riemannian prime distribution formula. We use the improved Riemannian prime distribution formula and merdens theorem to obtain a strong prime theorem.

Share and Cite:

Liu, D. (2021) Improved Riemann Ladder Function. Journal of Applied Mathematics and Physics, 9, 584-593. doi: 10.4236/jamp.2021.94042.

1. Introduction

The distribution of prime numbers is a difficult problem in number theory. The proof of prime theorem attracts many scholars and puzzles many wise people. Mathematicians have proved thousands of theorems on the premise of Riemann conjecture.

In 1737, the Swiss mathematician Leonhard Euler published a formula:

n = 1 n s = p ( 1 p s ) 1 , s > 1 , (1.1)

Here (1.1) is called Euler product formula. Where n is an integer, p is a prime and s is a real number.

Euler studied this formula. Let s = 1. For large x, Euler obtains an asymptotic formula

p x 1 p ~ log log x . (1.2)

From 1792 to 1793, Gauss, a German mathematician, studied the number of primes in 1000 adjacent integers around x. he found that for large x, the “average distribution density” of primes should be 1/logx, which is shown as follows [1] [2] [3] [4] [5]:

π ( x ) π ( x 1000 ) 1000 ~ 1 log x ,

where, the number of primes not exceeding x is π(x). General expression

π ( x ) π ( x r ) r ~ 1 log x , r x , (1.3)

Here (1.3) is Gauss’s guess. Represents the number of primes in r adjacent integers near x.

For example

Let x = 1238, r = 238, by (1.3) get

π ( 1238 ) π ( 1238 238 ) 238 = 203 168 238 = 0.14705 ~ 1 log 1238 = 0.14042 ,

According to (1.3), we can divide x into n parts. In this way, we can get π(x).

Let integern, x = nr, by (1.3) can get

π ( n r ) π ( n r r ) ~ r log ( n r ) , (1.4)

Set up n = 1 , 2 , 3 , , a , Substituting (1.4) get

n = 1 , π ( r ) π ( 0 ) ~ r log ( r ) , n = 2 , π ( 2 r ) π ( r ) ~ r log ( 2 r ) , n = 3 , π ( 3 r ) π ( 2 r ) ~ r log ( 3 r ) , n = a , π ( a r ) π ( a r r ) ~ r log ( a r ) ,

Add each

π ( r ) π ( 0 ) + π ( 2 r ) π ( r ) + π ( 3 r ) π ( 2 r ) + + π ( a r ) π ( a r r ) π ( a r ) = r log ( r ) + r log ( 2 r ) + r log ( 3 r ) + + r log ( a r ) ,

Let x = ar, get

π ( x ) ~ r n = 1 a 1 log ( n r ) , a = x r , (1.5)

Here (1.5) is a prime theorem without remainder. It is equal to the prime theorem π(x) ~ Lix.

For example

Let x = 64, r = 16, a = x/r = 4, by (1.6) calculation,

π ( x ) ~ 16 ( 1 log 16 + 1 log 32 + 1 log 48 + 1 log 64 ) ~ 18 ,

actual π ( 64 ) = 18 .

In this paper, we improve the Riemann ladder function and prove a strong prime number theorem by using the metens theorem [3].

π ( x ) = s ( x ) + O ( x 1 / 2 log x ) , ( x ) , (1.6)

s ( x ) = ( x + 1 ) log log x 2 n = 1 a log log ( 2 n ) , a = x / 2 ,

Before proving (1.6), we use a simple method to prove (1.3). Obviously, if we prove (1.3), we will prove (1.5).

2. The Elementary Proof of Prime Number Theorem

It is difficult to prove the prime theorem. It is especially difficult to prove the prime theorem in a simple way. We prove the prime theorem (1.3) by using the simple method according to Euler’s asymptotic formula

Prove

Let the integer r, prime p, xrpx obviously [6]

( x r ) x r p x 1 p π ( x ) π ( x r ) x x r p x 1 p , (2.1)

For example

Let x = 16, r = 12, xr = 4, 4 ≤ p ≤ 16, π(16) − π(4) = 4, by (2.1) get

4 ( 1 5 + 1 7 + 1 11 + 1 13 ) < π ( 16 ) π ( 4 ) < 16 ( 1 5 + 1 7 + 1 11 + 1 13 ) ,

By (2.1) can ge

x ( 1 r / x ) x r p x 1 p π ( x ) π ( x r ) x x r p x 1 p , (2.2)

Let r = x 1 / 2 , get

lim x ( 1 r / x ) = 1 , (2.3)

By (2.2) and (2.3) can get:

π ( x ) π ( x r ) ~ x x r p x 1 p , (2.4)

According to Euler’s asymptotic Formula (1.2), get

x r p x 1 p ~ log log x log log ( x r ) = log log x log ( x r ) = log log x log ( x ( 1 r / x ) ) = log log x log x + log ( 1 r / x ) ,

From this, we get

x r p x 1 p ~ log log x log x ( 1 x / r + 1 2 ( x / r ) 2 + 1 3 ( x / r ) 3 + + 1 k ( x / r ) k ) ,

Take main item 1 x / r , get:

x r p x 1 p ~ log log x log x 1 x / r = log 1 1 1 ( x / r ) log x = log ( 1 1 ( x / r ) log x ) 1 = 1 ( x / r ) log x + 1 2 ( x / r ) 2 ( log x ) 2 + 1 3 ( x / r ) 3 ( log x ) 3 + + 1 k ( x / r ) k ( log x ) k ,

Take the main item again 1 ( x / r ) log x , From this, we get

x r p x 1 p ~ r x log x , (2.5)

By (2.5) Substituting (2.4),

π ( x ) π ( x r ) ~ r log x , ( x ) ,

This proves (1.3).

Now, let’s look at the mean of the prime p, by (xr) ≤ px can get

p ( x r ) + x 2 = x r / 2 ,

According to (2.4), we change x into an average (xr/2) get

π ( x ) π ( x r ) ~ ( x r / 2 ) x r p x 1 p , (2.6)

For example x = 271, r = 34, 237 ≤ p ≤ 271, prime number 241, 251, 257, 263,

By (2.6) get

π ( 271 ) π ( 237 ) ~ 254 ( 1 / 241 + 1 / 251 + 1 / 257 + 1 / 263 ) = 4.02 ,

actual π(271) − π(237) = 4.

We divide x into n parts, n is a positive integer.

Let x = 2n, r = 2, can get (xr/2) = (2n − 1), by (2.6) we get

π ( 2 n ) π ( 2 n 2 ) = ( 2 n 1 ) 2 n 2 p 2 n 1 p ,

Therefore, we can improve the Riemann ladder function.

3. Improved Riemann Ladder Function

In 1859, German mathematician Riemann published a formula [2]:

π ( x ) = n = 1 k μ ( n ) J ( x 1 / n ) n , (3.1)

Here (3.1) is Riemann’s prime distribution formula [4]. Where J(x1/n) is called Riemannian ladder function [7]. The calculation is very complicated.

According to the definition of Riemann ladder function:

J ( x ) = π ( x ) + π ( x 1 / 2 ) 2 + π ( x 1 / 3 ) 3 + + π ( x 1 / k ) k , k = log x log 2 ,

A new ladder function is obtained by improving the ladder function

J ( x ) = 0 < p 2 1 p + 3 2 < p 4 1 p + 5 4 < p 6 1 p + + p 2 n 2 < p 2 n 1 p ,

Let x ≥ 2, Integer n, prime p, get [4]

π ( x ) = J ( x ) + 1 / 2 , (3.2)

J ( x ) = n = 1 a ( 2 n 1 ) 2 n 2 < p 2 n 1 p , a = x / 2 ,

Here (3.2) is an improved Riemannian prime distribution formula. It is strictly equal to π(x).

For example x = 16, a = 16/2 = 8, n = 1, 2, 3, 4, 5, 6, 7, 8.

According to (3.2), a prime number is crossed each time, and the calculation is as follows:

π ( 16 ) = 1 × 1 2 + 3 × 1 3 + 5 × 1 5 + 7 × 1 7 + 11 × 1 11 + 13 × 1 12 + 1 2 = 6 ,

In fact, π(16) = 6, it is now prove that (3.2).

Prove

If there is a prime p between 2n − 2 and 2n, it must be

p = 2 n 1 , (3.3)

Set up n = 1 , 2 , 3 , , a , by (3.2) and (3.3) get,

J ( x ) + 1 2 = 0 < p 2 1 p + 3 2 < p 4 1 p + 5 4 < p 6 1 p + + p 2 a 2 < p 2 a 1 p + 1 2 , (3.4)

According to (3.4), we can get the

J ( x ) + 1 2 = 1 × 1 2 + 3 × 1 3 + 5 × 1 5 + + p × 1 p + 1 2 = π ( x ) , (3.5)

By (3.5) get

π ( x ) = n = 1 a ( 2 n 1 ) 2 n 2 < p 2 n 1 p + 1 2 , (3.6)

Confirm (3.2) certification.

4. Merdens’ Theorem

In 1874, the mathematician Merdens proved that [7]:

lim x p x 1 p log log x = M , (4.1)

Here (4.1) is called: merdens theorem [8]. Where the mertensian constant

M = 0.2614972128476427837554268386086958 ,

Let the coefficient c (x) be obtained from (4.1)

p x 1 p log log x = c ( x ) , (4.2)

For the convenience of using (4.2), we transform (3.2).

For example x = 8, a = 8/2 = 4, by (3.2) can get

J ( 8 ) = 0 < p 2 1 p + 3 2 < p 4 1 p + 5 4 < p 6 1 p + 7 6 < p 8 1 p ,

Among

0 < p 2 1 p = p 2 1 p , 3 2 < p 4 1 p = 3 p 4 1 p 3 p 2 1 p , 5 4 < p 6 1 p = 5 p 6 1 p 5 p 4 1 p , 7 6 < p 8 1 p = 7 p 8 1 p 7 p 6 1 p ,

From the above

J ( 8 ) = p 2 1 p + 3 p 4 1 p 3 p 2 1 p + 5 p 6 1 p 5 p 4 1 p + 7 p 8 1 p 7 p 6 1 p = 7 p 8 1 p 2 p 2 1 p 2 p 4 1 p 2 p 6 1 p = 7 p 8 1 p 2 p 2 1 p 2 p 4 1 p 2 p 6 1 p 2 p 8 1 p + 2 p 8 1 p = 9 p 8 1 p 2 ( p 2 1 p + p 4 1 p + p 6 1 p + p 8 1 p ) ,

Generally speaking,

J ( x ) = ( x + 1 ) p x 1 p 2 ( p 2 1 p + p 4 1 p + p 6 1 p + + p 2 a 1 p ) , a = x / 2 , (4.3)

Replace (4.3) with (4.2), where

p x 1 p = log log x + c ( x ) ,

And

p 2 1 p + p 4 1 p + p 6 1 p + + p 2 a 1 p = log log 2 + c ( 2 ) + log log 4 + c ( 4 ) + log log 6 + c ( 6 ) + + log log x + c ( x ) = n = 1 a log log ( 2 n ) + n = 1 a c ( 2 n ) ,

It is obtained by substituting (4.3) above

J ( x ) = ( x + 1 ) ( log log x + c ( x ) ) 2 n = 1 a log log ( 2 n ) 2 n = 1 a ( 2 n ) , (4.4)

Namely

J ( x ) = ( x + 1 ) log log x 2 n = 1 a log log ( 2 n ) + ( x + 1 ) c ( x ) 2 n = 1 a ( 2 n ) , (4.5)

Let a = x/2, s(x) denote the logarithmic part, which is obtained by substituting (3.2) above

π ( x ) = s ( x ) + ( x + 1 ) c ( x ) 2 n = 1 x / 2 c ( 2 n ) + 1 2 , (4.6)

Among

s ( x ) = ( x + 1 ) log log x 2 n = 1 x / 2 log log ( 2 n ) ,

Here (4.6) is strictly equal to π(x).

For example x = 8, x/2 = 4, by (4.6) get

π ( 8 ) = 9 log log 8 2 n = 1 4 log log ( 2 n ) + 9 c ( 8 ) 2 n = 1 4 c ( 2 n ) + 1 2 = 9 log log 8 2 ( log log 2 + log log 4 + log log 6 + log log 8 ) + 9 c ( 8 ) 2 ( c ( 2 ) + c ( 4 ) + c ( 6 ) + c ( 8 ) ) + 0.5 ,

By (4.2) calculation coefficient

c ( 2 ) = p 2 ( 1 / p ) log log 2 = 0.8665129 c ( 4 ) = p 4 ( 1 / p ) ) log log 4 = 0.50669907 c ( 6 ) = p 6 ( 1 / p ) ) log log 6 = 0.45013525 c ( 8 ) = p 8 ( 1 / p ) log log 8 = 0.4440911

Get

π ( 8 ) = 9 ( 0.73209937 ) 2 ( 1.275418 ) + 9 ( 0.4440911 ) 2 ( 0.8665129 + 0.50669907 + 0.4501352 + 0.4440911 ) = 4.000005 = 4 ,

Actual π(8) = 4.

By (4.6) we can prove the prime theorem with remainder.

5. Prime Theorem with Remainder

The remainder estimation of prime number theorem is very complicated. The key is to use the mertensian constant to calculate the coefficient.

Let’s look at the remainder of (4.6). Among

2 n = 1 x / 2 c ( 2 n ) ,

Let even numbers y, x > y, according to (3.2) a = x/2, can get,

x / 2 = y / 2 + ( x y ) / 2 ,

From this we get

2 n = 1 x / 2 c ( 2 n ) = 2 n = 1 y / 2 c ( 2 n ) + 2 n = 1 ( x y ) / 2 c ( y + 2 n ) ,

Substituting (4.6),

π ( x ) = s ( x ) + ( x + 1 ) c ( x ) 2 n = 1 ( x y ) / 2 c ( y + 2 n ) + 2 n = 1 y / 2 c ( 2 n ) + 1 2 , (5.1)

where s(x) is the main term followed by the remainder. Now let’s look at the remainder of (5.1). Among

( x + 1 ) c ( x ) 2 n = 1 ( x y ) / 2 c ( y + 2 n ) , (5.2)

Here (5.2) is the key remainder. Let’s see c(x) and c(y + 2n).

In order to let c(y + 2n) approach M, let y = 2[x1/2]. According to (4.1) and (4.2), x tends to infinity, and we obtain

c ( x ) M c ( y + 2 n ) M ,

By (5.2) get

( x + 1 ) c ( x ) 2 n = 1 ( x y ) / 2 c ( y + 2 n ) ( x + 1 ) M 2 n = 1 ( x y ) / 2 M = ( x + 1 ) M ( x y ) M ,

Namely

( x + 1 ) c ( x ) 2 n = 1 ( x y ) / 2 c ( y + 2 n ) ( y + 1 ) M , (5.3)

Let’s look at a very small x.

For example x = 16, y = 4, (xy)/2 = 6, Substitute (5.3) and calculate according to (4.2)

( 16 + 1 ) c ( 16 ) 2 n = 1 6 c ( 4 + 2 n ) = ( 16 + 1 ) c ( 16 ) 2 ( c ( 6 ) + c ( 8 ) + c ( 10 ) + c ( 12 ) + c ( 14 ) + c ( 16 ) ) = 17 × 0.3242412 2 ( 0.45013 + 0.44409 + 0.342158 + 0.3568 + 0.373601 + 0.3242 ) = 5.5121 2 ( 2.2911 ) = 0.929918 ,

By (5.3) get

( y + 1 ) M = 5 × 0.2614972 1.3074861 ,

1.3074861 near 0.929918.

x The bigger, the closer (y + 1)M.

Now, let’s look at the scope of the remainder.

By (5.3) Substitute (5.1),

π ( x ) s ( x ) + ( y + 1 ) M 2 n = 1 y / 2 c ( 2 n ) + 1 2 , (5.4)

For the y of sufficient size, it is obvious that the estimation is

M < ( y + 1 ) M < y log x ,

Substitute (5.4) get

π ( x ) > s ( x ) + M 2 n = 1 y / 2 c ( 2 n ) , π ( x ) < s ( x ) + y log x 2 n = 1 y / 2 c ( 2 n ) , (5.5)

Let’s look at the coefficient,

2 n = 1 y / 2 c ( 2 n ) = 2 c ( 2 ) + 2 c ( 4 ) + 2 c ( 6 ) + + 2 c ( y ) ,

From (4.1) and (4.2) Confirm M/2 < c(2n) < 1, obviously

2 n = 1 y / 2 M 2 < 2 n = 1 y / 2 c ( 2 n ) < 2 n = 1 y / 2 1 ,

Get

y M 2 < 2 n = 1 y / 2 c ( 2 n ) < y , (5.6)

By (5.6) Substitute (5.5) get

π ( x ) > s ( x ) + M y > s ( x ) y log x , π ( x ) < s ( x ) + y log x y M 2 < s ( x ) + y log x , (5.7)

Let y = 2[x1/2], by (5.7) can get

π ( x ) = s ( x ) + O ( x 1 / 2 log x ) , (5.8)

By (5.8) Confirm, The prime theorem (1.6) is proved.

6. Conclusions

Previously, we discussed the Euler asymptotic formula

p x 1 p ~ log log x . (6.1)

From (6.2), we prove the prime theorem in a simple way

π ( x ) ~ r n = 1 x / r 1 log ( n r ) , (6.2)

Here (6.2) is a prime theorem without remainder.

We also discuss the distribution formula of prime number of Riemann

π ( x ) = n = 1 a μ ( n ) J ( x 1 / n ) n ,

We make a great improvement on the distribution formula of Riemannian prime numbers

π ( x ) = n = 1 a ( 2 n 1 ) 2 n 2 < p 2 n 1 p + 1 2 , (6.3)

This is a function strictly equal to π(x). In this way, we prove that Riemann’s conjecture is correct in principle.

According to the (6.3) theorem, we prove the strong prime number theorem

π ( x ) = s ( x ) + O ( x 1 / 2 log x ) , (6.4)

Here (6.6) is a prime theorem with remainder. Where s(x) and Lix are equal.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1] Atiyah, M. and Singer, I.M. (1963) The Index of Elliptic Operators on Compact Manifolds. Bulletin of the American Mathematical Society, 69, 422-433.
https://doi.org/10.1090/S0002-9904-1963-10957-X
[2] Ge, L.M. (2019) On the Riemann Zeta Function, I: KS-Transform. Acta Mathematica Sinica, Chinese Series, 62, 674.
[3] Lu, C.H. (2004) The Riemann Hypothesis (In Chinese). Tsinghua University Press, Beijing, Vol. 7, 9-18.
[4] Liu, J.F. (2013) Riemann Hypothesis Elementary Discussion. Progress in Applied Mathematics, 6, 74.
[5] Li, X. (1997) The Positivity of a Sequence of Numbers and the Riemann Hypothesis. Journal of Number Theory, 65, 325-333.
https://doi.org/10.1006/jnth.1997.2137
[6] Pan, C.D. and Pan, C.B. (1988) The Elementary Proof of Prime Number Theorem (In Chinese). Shanghai Science and Technology Press, Shanghai, 41-43.
[7] Rudin, W. (1987) Real and Complex Analysis, 3rd Edition, McGraw-Hill, New York.
[8] Riemann, B. (1859) Uber die Anzahl der Primzahlen unter Einer Gegebenen Grösse. Monatsber, Berlin Akad.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.