Fermat pseudoprimes: Difference between revisions

Line 236:

end function

~~sequence~~ limits = {12000, 25000, 50000, 100000}

constant limits = {12000, 25000, 50000, 100000}, sp = repeat('-',43)

printf(1,"Base <%s first 20 %s> <=%s\n",{sp,sp,join(limits,fmt:="%6d")})

for base=1 to 20 do

integer count = 0, nlx = 1, nl = limits[1]

Line 253:

Line 254:

end if

end for

string s = join(shorten(first20,"",3,"%d"),","),

string s = join(shorten(first20,"",8,"%5d"),","),

c = join(~~columnize({limits,~~counts~~}),", "~~,fmt:="%~~d:%5d~~")

c = join(counts,fmt:="%6d")

printf(1,"~~Base~~ %2d~~, counts to~~ %s, ~~first~~ ~~20:~~ %s\n", {base, c, s})

printf(1,"%2d: %s %s\n", {base, s, c})

end for

<pre>

Base <------------------------------------------- first 20 -------------------------------------------> <= 12000 25000 50000 100000

Base 1, counts to 12000:10561, 25000:22237, 50000:44866, 100000:90407, first 20: 4,6,8,...,28,30,32

1: 4, 6, 8, 9, 10, 12, 14, 15,..., 22, 24, 25, 26, 27, 28, 30, 32 10561 22237 44866 90407

Base 2, counts to 12000: 25, 25000: 38, 50000: 55, 100000: 78, first 20: 341,561,645,...,6601,7957,8321

2: 341, 561, 645, 1105, 1387, 1729, 1905, 2047,..., 4033, 4369, 4371, 4681, 5461, 6601, 7957, 8321 25 38 55 78

Base 3, counts to 12000: 25, 25000: 38, 50000: 53, 100000: 78, first 20: 91,121,286,...,4961,5551,6601

3: 91, 121, 286, 671, 703, 949, 1105, 1541,..., 2701, 2821, 3281, 3367, 3751, 4961, 5551, 6601 25 38 53 78

Base 4, counts to 12000: 50, 25000: 75, 50000: 111, 100000: 153, first 20: 15,85,91,...,1905,2047,2071

4: 15, 85, 91, 341, 435, 451, 561, 645,..., 1387, 1581, 1695, 1729, 1891, 1905, 2047, 2071 50 75 111 153

Base 5, counts to 12000: 22, 25000: 35, 50000: 54, 100000: 73, first 20: 4,124,217,...,8029,8911,9881

5: 4, 124, 217, 561, 781, 1541, 1729, 1891,..., 5662, 5731, 6601, 7449, 7813, 8029, 8911, 9881 22 35 54 73

Base 6, counts to 12000: 31, 25000: 46, 50000: 74, 100000: 104, first 20: 35,185,217,...,4123,4495,5713

6: 35, 185, 217, 301, 481, 1105, 1111, 1261,..., 2821, 3421, 3565, 3589, 3913, 4123, 4495, 5713 31 46 74 104

Base 7, counts to 12000: 21, 25000: 32, 50000: 49, 100000: 73, first 20: 6,25,325,...,10585,10621,11041

7: 6, 25, 325, 561, 703, 817, 1105, 1825,..., 4525, 4825, 6697, 8321,10225,10585,10621,11041 21 32 49 73

Base 8, counts to 12000: 76, 25000: 110, 50000: 150, 100000: 218, first 20: 9,21,45,...,651,861,949

8: 9, 21, 45, 63, 65, 105, 117, 133,..., 481, 511, 561, 585, 645, 651, 861, 949 76 110 150 218

Base 9, counts to 12000: 55, 25000: 83, 50000: 113, 100000: 164, first 20: 4,8,28,...,1036,1105,1288

9: 4, 8, 28, 52, 91, 121, 205, 286,..., 671, 697, 703, 946, 949, 1036, 1105, 1288 55 83 113 164

Base 10, counts to 12000: 35, 25000: 53, 50000: 65, 100000: 90, first 20: 9,33,91,...,3367,4141,4187

10: 9, 33, 91, 99, 259, 451, 481, 561,..., 1729, 2409, 2821, 2981, 3333, 3367, 4141, 4187 35 53 65 90

Base 11, counts to 12000: 30, 25000: 41, 50000: 61, 100000: 89, first 20: 10,15,70,...,2821,4577,4921

11: 10, 15, 70, 133, 190, 259, 305, 481,..., 1330, 1729, 2047, 2257, 2465, 2821, 4577, 4921 30 41 61 89

Base 12, counts to 12000: 35, 25000: 60, 50000: 91, 100000: 127, first 20: 65,91,133,...,2233,2465,2701

12: 65, 91, 133, 143, 145, 247, 377, 385,..., 1649, 1729, 1885, 1891, 2041, 2233, 2465, 2701 35 60 91 127

Base 13, counts to 12000: 31, 25000: 51, 50000: 68, 100000: 91, first 20: 4,6,12,...,3605,5028,5149

13: 4, 6, 12, 21, 85, 105, 231, 244,..., 1099, 1785, 1891, 2465, 2806, 3605, 5028, 5149 31 51 68 91

Base 14, counts to 12000: 33, 25000: 51, 50000: 69, 100000: 96, first 20: 15,39,65,...,2665,2743,3277

14: 15, 39, 65, 195, 481, 561, 781, 793,..., 1541, 1891, 2257, 2465, 2561, 2665, 2743, 3277 33 51 69 96

Base 15, counts to 12000: 22, 25000: 31, 50000: 42, 100000: 70, first 20: 14,341,742,...,8701,8911,9073

15: 14, 341, 742, 946, 1477, 1541, 1687, 1729,..., 3277, 4187, 6541, 6601, 7471, 8701, 8911, 9073 22 31 42 70

Base 16, counts to 12000: 69, 25000: 102, 50000: 145, 100000: 200, first 20: 15,51,85,...,1387,1581,1687

16: 15, 51, 85, 91, 255, 341, 435, 451,..., 1105, 1247, 1261, 1271, 1285, 1387, 1581, 1687 69 102 145 200

Base 17, counts to 12000: 31, 25000: 44, 50000: 63, 100000: 83, first 20: 4,8,9,...,4005,4033,4187

17: 4, 8, 9, 16, 45, 91, 145, 261,..., 1729, 1885, 2149, 2821, 3991, 4005, 4033, 4187 31 44 63 83

Base 18, counts to 12000: 46, 25000: 69, 50000: 98, 100000: 134, first 20: 25,49,65,...,1649,1729,1921

18: 25, 49, 65, 85, 133, 221, 323, 325,..., 931, 1105, 1225, 1369, 1387, 1649, 1729, 1921 46 69 98 134

Base 19, counts to 12000: 48, 25000: 70, 50000: 93, 100000: 121, first 20: 6,9,15,...,1661,1849,1891

19: 6, 9, 15, 18, 45, 49, 153, 169,..., 905, 906, 1035, 1105, 1629, 1661, 1849, 1891 48 70 93 121

Base 20, counts to 12000: 35, 25000: 49, 50000: 66, 100000: 101, first 20: 21,57,133,...,2821,2947,3059

20: 21, 57, 133, 231, 399, 561, 671, 861,..., 1891, 2059, 2413, 2501, 2761, 2821, 2947, 3059 35 49 66 101

</pre>

=={{header|Raku}}==

<lang perl6>use List::Divvy;