这个问题有现成的算法解决吗？我以前似乎看到过大整数加密，就是分解质因数。

没有现成的算法啊，C语言很容易实现，但一碰到fortran编程就搞不定了。多谢啊

只能穷举

判断小于此数的所有数是否是质因子

穷举的话
只要判断小于其 1/2 的整数就行了

How about the following codes？
------------------------------------------

program factt
implicit none
integer :: n,i,isprime

do 
      read*,n
      n=abs(n)
      if( n.le. 1) print*,n
100   do i=2,n
           if(mod(n,i) .eq. 0)then
               if(isprime(i)) then
                    print*,i            
                    n=n/i
                    goto 100
                endif
            endif
       enddo   
enddo
end program

integer function isprime(k)
implicit none
integer::k,i,j
i=abs(k)
if(i .le. 1)then
    isprime=0
    return
elseif(i .eq. 2)then
    isprime=1
    return
elseif(mod(i,2) .eq. 0) then
    isprime=0
    return
else
    do j=3, i-2,2
        if(mod(i,j) .eq. 0) then
            isprime=0
            return
        endif
    enddo
    isprime=1
    return
endif
end function

穷举的话

只要判断小于其[color=FF0000]平方根[/color]的整数就行了

目前还没有非常快速的算法

正因为如此

一个非常大的整数找到其所有质因子将花费相当长的时间

现代加密技术就是运用此特性

6楼的说法有问题啊

比如 数为34 质因子为 17 和 2
其平方根小于6 按你的讲法 17 是求不出来的

program primefactor
   implicit none
   integer :: n, i
   logical :: isPrime
   do 
      print *, 'please input an integer'
      read  *, n
      if(n == 0) exit
      i = 2
      do while(i<=n)
         if( mod(n,i) == 0 ) then
                 print *, i
                 n = n/i
             if(n == 1) then
              exit
             else
              i = 2
             endif
           else
           i=i+1
           endif
      enddo
   enddo
end program

这样不用isprime过程的

[quote]program primefactor
   implicit none
   integer :: n, i
   logical :: isPrime
   do 
      print *, 'please input an integer'
      read  *, n
      if(n == 0) exit
      i = 2
      do while(i<=n)
         if( mod(n,i) == 0 ) then
                 print *, i
                 n = n/i
             if(n == 1) then
              exit
             else
              i = 2
             endif
           else
           i=i+1
           endif
      enddo
   enddo
end program

这样不用isprime过程的
[/quote]


对楼上稍作修改，对于比较大的数，效率的提高还是比较明显的，不知妥否？

      program primefactor
      implicit none
      integer(kind=8) :: i,j,k,n
      
      do
       print *, 'please input an integer'
       read  *, n
       j = n
       k = 1 
       if(j == 0) exit
       i = 2
       do while(i<=int(sqrt(real(j))))
         if( mod(j,i) == 0 ) then
            print *, i
            k = k*i
            j = j/i
            if(j == 1) exit
         else
           i = i+1
         endif
      enddo
      print *, n/k

     enddo
     end program

有兴趣的网友可以验证，不妥之处，请不吝指正！