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

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次