孪生素数的计算

无心人

呵呵

没想那么多
遇到有交流
俺就想说话
以后俺注意了

呵呵，不知道俺有腰斩的权力么？

无心人

呵呵

tprime：
我6#问的是你是完全筛掉所有可能素因子
还是筛一部分，再过米勒－罗宾测试

另外，米勒－罗宾测试对64位
我估计是大概64 / 2 * 3 * (4 + 3) * C次CPU运算
就能完成

C估计值7 <= C <= 12

gxqcn

原帖由 无心人 于 2008-7-4 19:51 发表
呵呵，不知道俺有腰斩的权力么？

版主应该都有该权力吧。

你把鼠标移到“管理选项”上，看浮动菜单上是否有“分割主题”即可得知。

无心人

有的

呵呵

知道了

tprime

我的程序没有用到素数测试方法, 将所有的因子都筛出掉
程序经过相当长的优化, 才达到目前PI2(10^16)2000小
时的结果.当然算法是最重要的,这半年都天花点时间
优化语言方面.
  btw: 打算在新的版本中试用素数测试部分, 这部分对性能
要求非常高, 希望每秒能达到10000个 10^14 - 10^17 素性
测试才能达到要求, 否则不如直接全筛

无心人

我问你，全筛每秒是淘汰多少素数？？

无心人

刚去查旧帖子
128bit的乘法最快是
$4 xx 10^-8$秒一次
测试一次素数，最多需要计算
256次乘法和同样多的除法
假设除法时间是乘法的16倍
则共需要256*17次乘法
总时间是0.0017
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
但考虑64位
则乘法应该是128位的1/4
（实际更小吧）
且最多是128个乘和同样的除
考虑64位除的时间是8个乘法时间
则总共是128*9个乘法
即0.00001152秒
＝＝＝＝＝＝＝＝＝＝＝＝
虽然计算很粗糙
但我想能达到你的要求

无心人

刚实际测试
发现64bit整数进行一次素性检测
花费是0.0001秒
但应该还可以优化

无心人

tprime有100万内的四生素数列表么？
只要第一个数就可以了

tprime

我给你10亿的吧