无心人 发表于 2008-10-11 20:32:41

上面贴的论文有如何计算伪素数的算法
用1.2G CPU计算到10^15花费8.6天
我想如果精心编码加更快的CPU
有希望在60天内计算出10^16的结果
呵呵

medie2005 发表于 2008-10-11 20:51:43

现在,顺带求一下下面的问题:
分别求出第一次误判的数大于10^8,10^9,10^10,10^11,10^12,10^13,10^14,10^15 的测试基数目最少的测试基组合.

无心人 发表于 2008-10-11 20:58:14

刚去测试你提供的基的
(2, 3,102199) 通过5262
(2, 3, 696811) 通过5398
都是很多的
所以目前我想并没有很快捷的方法来筛选基组合

无心人 发表于 2008-10-11 21:00:26

10^13应该能找到4元基组合的

无心人 发表于 2008-10-11 21:06:02

被你杀掉的基首数字就不错
(2,3,7,61,24251)
669094855201

无心人 发表于 2008-10-11 21:09:01

*Primes> let t8 = filter (<100000000) t23
*Primes> t8
[1373653,1530787,1987021,2284453,3116107,5173601,6787327,11541307,13694761,15978
007,16070429,16879501,25326001,27509653,27664033,28527049,54029741,61832377,6609
6253,74927161,80375707]
*Primes> length t8
21

10^8内仅21组通过2,3的

medie2005 发表于 2008-10-11 21:09:52

原帖由 无心人 于 2008-10-11 20:58 发表 http://bbs.emath.ac.cn/images/common/back.gif
刚去测试你提供的基的
(2, 3,102199) 通过5262
(2, 3, 696811) 通过5398
都是很多的
所以目前我想并没有很快捷的方法来筛选基组合
我现在觉得,小素数的筛选功能更强一些,这个可能也是(2, 3,102199) 和(2, 3, 696811) 误判数较多,以及(2,3,5),(2,3,13),(2,3,61)误判数较少的原因。

无心人 发表于 2008-10-11 21:13:55

10^8的得到了
(2, 3, 23)

medie2005 发表于 2008-10-11 21:16:47

原帖由 无心人 于 2008-10-11 21:06 发表 http://bbs.emath.ac.cn/images/common/back.gif
被你杀掉的基首数字就不错
(2,3,7,61,24251)
66,90948,55201
呵呵,这个(2,3,7,61,24251)还是一般般.实质上我在17#就找到了一个比它好的了:
(2,3,5,7,3439829)第一次误判发生在9905,03702,60227.

mathe 发表于 2008-10-11 21:18:51

原帖由 无心人 于 2008-10-11 20:21 发表 http://bbs.emath.ac.cn/images/common/back.gif
http://oldweb.cecm.sfu.ca/pseudoprime/psp-search-slides.pdf

541
这个论文好像还没有用到我142#的结论。所以应该还可以改善
页: 7 8 9 10 11 12 13 14 15 16 [17] 18 19 20 21 22
查看完整版本: 能通过2,3,5,7的检验的合数