如何实现对大数的素性检测？

306908677

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

306908677

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

fcode

大数模块需要修改一个 parameter 常量来达到你需要的位数。

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

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

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

珊瑚虫

问题比较高大上， 主站里面有个大数求阶乘的帖子，问题和你的类似，
http://fcode.cn/algorithm-50-1.html
希望对你有帮助，否则再讨论