![]() |
|
|
#34 |
|
Romulan Interpreter
"name field"
Jun 2011
Thailand
240638 Posts |
With a good eps viewer (I use acdSee) this pari liner will generate in just 4 seconds a very beautiful plot which can be scrolled, panned, and zoomed in. To have an idea how this function looks. It will "go under" an infinite number of times, but how much it will stay under when the primes get really HUGE, nobody knows. In the ranges it was tested, it is mostly positive, as already mentioned.
Code:
gp > default(realprecision,3); maxprime=5*10^6; v=primepi(maxprime); x=vector(v); y=vector(v); for(i=1,v,x[i]=i); i=0; c1=0; c5=0; forprime(p=5,maxprime,if(p%6==1,c1++,c5++);i++; y[i]=c5-c1); plothraw(x,y) Edit: of course the "implementation" is far away of being efficient, a sieve would do wonders, hihi, but pari stores a large number of small primes in a vector. You may need "allocatemem()" few times if you want to go larger. But yet, this is only didactic, it is too slow for "large study". Last fiddled with by LaurV on 2021-08-27 at 05:36 |
|
|
|
|
|
#35 |
|
"ม้าไฟ"
May 2018
2×3×89 Posts |
|
|
|
|
|
|
#36 | |
|
Apr 2020
100001000112 Posts |
Quote:
may have something to do with this.In this case, OEIS is a good place to look: an OEIS search for the first crossover point, 608981813029, might find the sequence of crossover points. Indeed, it gives A007352 as the first result. The b-file tells us that there are over 9000 further crossovers up to 610968213803, at which point 3n+2 takes the lead until 6148171711663 (note this has one more digit than the previous two numbers!). A Google search for 6148171711663 turns up this presentation. They verify this value, and presumably the first crossover too, although they claim that there were errors in the search that found it. I assume this means some of the terms in A007352 beyond 6148171711663 are wrong, unless the b-file has been corrected. |
|
|
|
|
|
|
#37 |
|
"ม้าไฟ"
May 2018
2·3·89 Posts |
I will reach the crossover point in days with two Raspberry Pi 4B devices running the same task. Of course, I could do that much faster with my gaming PC, but I am also testing a custom solar-powered configuration independent of the power grid. One Raspberry Pi 4B device is connected to the power grid and another one to a battery charged by a solar panel. After the testing phase is over and hopefully both devices produce the same result, I will be able to arrange an array of Edge computing devices for prime testing for which there will be no need to pay electricity bills.
|
|
|
|
|
|
#38 |
|
"ม้าไฟ"
May 2018
10268 Posts |
The computational experiment was able to reach the first crossover point x = 608,981,813,017 (π(x) = 23,338,590,791) where π6,5(x) - π6,1(x) = -1.
This is preceded by three equilibrium points x = 608,981,812,891, 608,981,812,951, and 608,981,812,993 (π(x) = 23,338,590,786, 23,338,590,788, and 23,338,590,790) for which π6,5(x) - π6,1(x) = 0. It appears that there is an existing OEIS A096629 (https://oeis.org/A096449) description of the equilibrium points (and a list of 85,508 equilibrium points up to x = 6,156,051,951,809 (π(x) = 216,682,882,516)) with the following description: "Values of n for which {p_3, p_4, ..., p_n} (mod 3) contains equal numbers of 1's and 2's." I was able to find it after a Google search with the π(x) values of equilibrium points obtained with the Raspberry Pi devices. The current version of the Wolfram code generates consecutive lists of 10,000,000 consecutive primes for which a prime count of π6,5(x) and π6,1(x) is performed. |
|
|
|
|
|
#39 |
|
"ม้าไฟ"
May 2018
21616 Posts |
The prime number races draw an analogy with the Big Bang of the universe and the existence of more matter than antimatter in the current universe.
What if the unknown fabric of the universe creates new universes with Big Bangs after consecutive equilibrium zero points and in some universes the antimatter dominates over matter? Last fiddled with by Dobri on 2021-09-18 at 14:43 |
|
|
|
|
|
#40 |
|
"ม้าไฟ"
May 2018
2×3×89 Posts |
The maximum difference observed in the entire interval before the first crossover point, π(x) < 23,338,590,791, is Max[π6,5(x) - π6,1(x)] = 47,050 = 2×52×941.
Correction to post #38: The link to A096629 is https://oeis.org/A096629. Last fiddled with by Dobri on 2021-09-18 at 15:59 |
|
|
|
|
|
#41 |
|
"ม้าไฟ"
May 2018
2×3×89 Posts |
The maximum difference Max[π6,5(x) - π6,1(x)] = 47,050 = 2×52×941 occurs at x = 457,861,654,499 (π(x) = 17,741,340,390).
The primes 2, 5, 941, and 457,861,654,499 are congruent to 2 (mod 3) thus 5, 941, and 457,861,654,499 are of the type 6k-1 = 3(2k-1)+2. Therefore, said primes are also Eisenstein primes with zero imaginary part (see https://mathworld.wolfram.com/EisensteinPrime.html). Last fiddled with by Dobri on 2021-09-18 at 19:05 |
|
|
|
|
|
#42 |
|
"ม้าไฟ"
May 2018
10268 Posts |
As the computational experiment continues beyond the first crossover point (and multiple subsequent crossovers are encountered), the observed minimum Min[π6,5(x) - π6,1(x)] is equal to -1,539 = -34×19 so far but its value may change with the future updates on this prime number race.
|
|
|
|
|
|
#43 |
|
"ม้าไฟ"
May 2018
2×3×89 Posts |
The following local maxima were observed:
π6,5(x) - π6,1(x) = 47,716 = 22×79×151, x = 683,008,329,317, π(x) = 26,060,805,816; π6,5(x) - π6,1(x) = 48,542 = 2×13×1,867, x = 684,347,039,021, π(x) = 26,109,930,072; π6,5(x) - π6,1(x) = 48,591 = 32×5,399, x = 684,349,485,899, π(x) = 26,110,020,031; and π6,5(x) - π6,1(x) = 48,910 = 2×5×67×73, x = 684,706,312,967, π(x) = 26,123,111,894. The primes 2, 5, and 5,399 are congruent to 2 (mod 3), then 13, 67, 73, 79, 151, and 1, 867 are congruent to 1 (mod 3), and 3 is congruent to 0 (mod 3). The observed difference 48,910 is the largest one since the computational experiment was initiated. |
|
|
|
|
|
#44 |
|
"Serge"
Mar 2008
San Diego, Calif.
240358 Posts |
Well, isn't it fun racing the race that already finished in 1976 and reporting results as if they were happening right now, with great urgency, with minute-by-minute updates and all?
Hunter S. Thompson you are not, Sir. This is no Las Vegas, ...only fear and loathing. P.S. Your "blogging" belongs in Blogorrhea section (where the blogger is the only one reading). Blog as much as you want. |
|
|
|
![]() |
| Thread Tools | |
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| Number of primes 2^n-1 within certain range of n | bur | Homework Help | 8 | 2021-02-18 20:54 |
| Number of relative primes for P-1 | petrw1 | Factoring | 4 | 2018-11-07 20:20 |
| The expected number of primes | Batalov | Computer Science & Computational Number Theory | 5 | 2016-08-11 01:17 |
| Estimating the number of primes in a partially-factored number | CRGreathouse | Probability & Probabilistic Number Theory | 15 | 2014-08-13 18:46 |
| Guess a Number | davar55 | Puzzles | 11 | 2008-03-19 21:33 |