如何实现对大数的素性检测？
RT望大牛不吝赐教］
附上新手写的一段程序(被整型数限制完全算不了大数)＝＝抛砖引玉
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
问题比较高大上， 主站里面有个大数求阶乘的帖子，问题和你的类似，
http://fcode.cn/algorithm-50-1.html
希望对你有帮助，否则再讨论 大数模块需要修改一个 parameter 常量来达到你需要的位数。
（再大也不会是无穷的，要理解，计算机资源永远是有限的）
代码中的 nr_of_decimal_digits = 180 ，改为更大的值，比如 = 9999 。 就可达到 10**（9999）。
最多可达到 10**（2**31），我不确定是否能满足你的要求，但绝不限于 p=13【 2**（2*13）】
本帖最后由 306908677 于 2014-12-8 14:06 编辑
谢谢各位，改了一下算法，可以算到2203次
页:
[1]