[通用算法] 如何实现对大数的素性检测？

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