Cubic Frobenius test statistics to $10^{11}$

I ran my latest prime-testing program, based on a cubic Frobenius test, on the numbers from 1 to $10^{11}$. Numbers in that range coprime to 30 are considered for testing. Given the input $n$, the program looks for the first polynomial $f_{i}$ such that $f_{i}$ is irreducible over $F_{n}$ , where $f_{i}$ is one of the 23 polynomials listed at the end of the earlier post https://dbernier.ca/2024/11/10/what-is-a-cubic-frobenius-primality-test/. Irreducibility is determined by table lookup. Occasionally, no polynomial is irreducible over $F_{n}$. Most often, in that case $n$ is composite although there are a few rare primes. Once the polynomial is selected, three congruences modulo $n$ are tested. We can represent any test result by a 3-bit binary sequence where a 1 bit means that a congruence is satisfied and a 0 bit means the congruence is not satisfied. The possible test results are thus $000, 001, 010, 011, 100, 101, 110, 111$ which can be read as an integer from 0 to 7. Then the summary statistics are as follows:

Type 0: 22548610645
Type 1: 16
Type 2: 154
Type 3: 0
Type 4: 4
Type 5: 0
Type 6: 0
Type 7: 4118054810
Grand total of tested and untested numbers between 1 and 100000000000: 26666666666
Grand total of probable primes found and untested primes, coprime to 30,  between 1 and 100000000000: 4118054810


There are 22,548,610,645 results of type 0 where none of the three conditions is satisfied. There are 4,118,054,810 results of type 7 where all three conditions are satisfied. This is 3 less than $\pi(10^{11})=4,118,054,813$ because the primes 2, 3 and 5 (not being coprime to 30) are not even considered for testing. There are 154 results of type 2, 16 of type 1 and 4 of type 4 (these three types correspond to 1 of 3 conditions being satisfied). For types 3, 5 and 6 there are no test results (these types correspond to 2 of 3 conditions being satisfied). I might modify the program to show how many numbers with a type 7 result are composite, although this might slow down the program. The same primality test has been conducted for numbers up to $10^{13}$ with no pseudoprimes found (ie type 7 results for a composite input). A recent version of the program is available on Github here: https://github.com/mariotrevi/cubicFrobenius. Below, we have the program’s output relating to composites with test results of 1 through 6. But first, let us provide the summary statistics to $10^{10}$ for an updated version of the program that counts how many numbers with a type 7 result are prime, and that shows the prime count according to Kim Walisch’s primesieve (the job took 63 minutes): (and secondly, the summary statistics to $10^{11}$, a job that took 24 hours on my PC)


Composite and prime counts by type between 1 and 10000000000:
Type 0: 2211613577
Type 1: 13
Type 2: 83
Type 3: 0
Type 4: 4
Type 5: 0
Type 6: 0
Type 7: 455052508
Type 7 primes: 455052508
Validating Type 7 primes using Primesieve:
primesieve prime count: 455052511


Composite and prime counts by type between 1 and 100000000000:
Type 0: 22548610645
Type 1: 16
Type 2: 154
Type 3: 0
Type 4: 4
Type 5: 0
Type 6: 0
Type 7: 4118054810
Type 7 primes: 4118054810
Validating Type 7 primes using Primesieve:
primesieve prime count: 4118054813



Processing range: 1 to 100000000000
1243 test_type=1
3059 test_type=2
289 test_type=2
6859 test_type=2
2197 test_type=2
37249 test_type=1
50653 test_type=2
106597 test_type=2
237169 test_type=2
465283 test_type=2
573067 test_type=2
912673 test_type=2
935089 test_type=1
2813177 test_type=4
3464161 test_type=2
4330747 test_type=2
4826809 test_type=2
5603333 test_type=2
8025641 test_type=2
10991701 test_type=2
14687681 test_type=2
16636201 test_type=1
18210401 test_type=2
24242261 test_type=2
25916687 test_type=2
34173553 test_type=1
36487547 test_type=4
39866533 test_type=2
44482901 test_type=2
47045881 test_type=2
45407909 test_type=1
57833861 test_type=2
67194401 test_type=2
69348869 test_type=1
83403041 test_type=2
86766121 test_type=2
103253333 test_type=2
109393201 test_type=2
119145439 test_type=2
136813753 test_type=2
140558671 test_type=2
140400719 test_type=1
159087769 test_type=4
176304427 test_type=4
238244041 test_type=2
240141287 test_type=2
303651041 test_type=2
303812587 test_type=2
348668321 test_type=2
467831377 test_type=2
481637381 test_type=2
500985941 test_type=2
523501693 test_type=1
624306439 test_type=2
787658257 test_type=2
818131609 test_type=1
878940833 test_type=2
999260501 test_type=2
1020221461 test_type=2
1048435417 test_type=2
1294261781 test_type=2
1284365503 test_type=2
1603209889 test_type=2
1870920577 test_type=2
2301980801 test_type=2
2462601379 test_type=1
2565726409 test_type=2
2871974561 test_type=2
2822739731 test_type=2
2973097121 test_type=2
3233558021 test_type=2
3616574081 test_type=2
3684692249 test_type=2
4196753009 test_type=2
4282136321 test_type=2
4310746133 test_type=2
4631740961 test_type=2
4683619717 test_type=1
5088036101 test_type=2
5165622721 test_type=2
5267500421 test_type=2
5427083201 test_type=2
5386796083 test_type=2
5500828481 test_type=2
5687368961 test_type=2
6399624533 test_type=2
6540055061 test_type=2
6791517697 test_type=2
6830494301 test_type=2
7115167081 test_type=2
7054034261 test_type=2
7182672721 test_type=2
7339666181 test_type=2
7343216063 test_type=2
8057220751 test_type=2
8464140833 test_type=2
8841126533 test_type=2
9495516163 test_type=1
9793709857 test_type=2
9711692033 test_type=2
10122343561 test_type=2
10658188801 test_type=2
10441867051 test_type=2
10604499373 test_type=2
12018820997 test_type=2
12175375807 test_type=2
13052405521 test_type=2
13205640533 test_type=2
13931904533 test_type=2
14436868181 test_type=2
14768399381 test_type=2
15032007649 test_type=2
15282655921 test_type=2
15058173151 test_type=2
15658898677 test_type=2
16521529753 test_type=2
16950083333 test_type=2
18164433281 test_type=2
18852287831 test_type=2
19875403241 test_type=2
21736583681 test_type=2
22247113321 test_type=2
22539300841 test_type=2
22712568341 test_type=2
23131319051 test_type=2
23508787541 test_type=2
25370292641 test_type=2
25769117921 test_type=2
27529550929 test_type=2
27108673453 test_type=1
27934820033 test_type=2
29412900833 test_type=2
30929084033 test_type=2
31391368961 test_type=2
31840678801 test_type=2
30983938133 test_type=2
32225407061 test_type=2
33378012577 test_type=2
33331696133 test_type=2
33995233301 test_type=2
34886240033 test_type=2
36930151301 test_type=2
39654644177 test_type=2
43035037781 test_type=2
43606069781 test_type=2
44430236033 test_type=2
43971413333 test_type=2
45326883931 test_type=1
46450465601 test_type=2
47718880301 test_type=2
48046912607 test_type=2
49996992997 test_type=2
50611638533 test_type=2
52775958313 test_type=2
54280972481 test_type=2
56112367501 test_type=2
60224864701 test_type=1
62098032641 test_type=2
63598662401 test_type=2
65851193633 test_type=2
66943750721 test_type=2
69254252033 test_type=2
71266006981 test_type=2
73838796571 test_type=2
81964996961 test_type=2
83664288001 test_type=2
82555158533 test_type=2
86498244001 test_type=2
86910928133 test_type=2
92361442181 test_type=2
93463810133 test_type=2
93399952367 test_type=2
97646190881 test_type=2
97743638021 test_type=2
meditationatae's avatar

By meditationatae

Canadian

Discover more from meditationatae

Subscribe now to keep reading and get access to the full archive.

Continue reading