Fortran Coder

标题: 如何实现对大数的素性检测？ [打印本页]

作者: 306908677 时间: 2014-12-2 21:56
标题: 如何实现对大数的素性检测？
RT
望大牛不吝赐教］
附上新手写的一段程序(被整型数限制完全算不了大数)＝＝抛砖引玉

[Fortran] 纯文本查看 复制代码



program main



implicit none



integer(kind=8),Mp,Ln,q

integer(kind=4) ,parameter  ::  p=19

integer ,parameter :: Lo=4

integer :: n,i

Mp=1

do i=1,p

Mp=Mp*2

end do



Mp=Mp-1

Ln=Lo

n=p-2

loop:do i=1,n

Ln=Ln**2-2

q=mod(Ln,Mp)



if(q==0)then

call print_big(Mp)

else 

cycle loop

end if

end do loop



end program main

作者: 珊瑚虫 时间: 2014-12-2 22:21
问题比较高大上，主站里面有个大数求阶乘的帖子，问题和你的类似，
http://fcode.cn/algorithm-50-1.html
希望对你有帮助，否则再讨论

作者: fcode 时间: 2014-12-3 08:30
大数模块需要修改一个 parameter 常量来达到你需要的位数。

（再大也不会是无穷的，要理解，计算机资源永远是有限的）

代码中的 nr_of_decimal_digits = 180 ，改为更大的值，比如 = 9999 。就可达到 10**（9999）。

最多可达到 10**（2**31），我不确定是否能满足你的要求，但绝不限于 p=13【 2**（2*13）】

作者: 306908677 时间: 2014-12-4 11:30
本帖最后由 306908677 于 2014-12-8 14:06 编辑

谢谢各位，改了一下算法，可以算到2203次

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)