[½»Á÷] ¡¾ÇóÖú¡¿³ÌÐòÎÞ´íÎ󣬵«ÊÇÔËÐеÄÊý¾Ý²»¶Ô£¬ÊÇΪʲô£¿¡¾ÒÑÍê½á¡¿ ÒÑÓÐ1È˲ÎÓë

SUBROUTINE Power(m,last,term,length)
IMPLICIT NONE
INTEGER::m,last,length,i,j,k
INTEGER(2),DIMENSION(1:last)::term
term=0;term(last)=1
doi: DO i=1,m
      DO j=1,last
           IF(TERM(j)/=0)THEN
            DO k=last,j,-1
                  term(k)=term(k)**2
                END DO
                CYCLE
           END IF
          END DO
        DO k=last,j-1,-1
           term(k-1)=term(k-1)+term(k)/10
           term(k)=mod(term(k),10)
        END DO
   ENDDO doi
   DO k=1,last
     IF(term(k)/=0)THEN
           length=k;exit
         END IF
        END DO
        END SUBROUTINE Power
!
   PROGRAM Mersenne_Number
   IMPLICIT NONE
   INTEGER::m,n,k,i
   INTEGER,PARAMETER::max=20
   INTEGER(2),DIMENSION(3)::p=(/19,31,61/)
   INTEGER(2),DIMENSION(1:max)::term
   OPEN(1,FILE='MYF1.OUT')
   DO i=1,3
     m=p(i);CALL Power(m,max,term,n)
         term(max)=term(max)-1
         WRITE(*,'("P=",I3,"Mersenne=",20I1)')p(i),(term(k),k=n,max)
         WRITE(1,'("P=",I3,"Mersenne=",20I1)')p(i),(term(k),k=n,max)
        END DO
        END PROGRAM Mersenne_Number

Originally posted by azengjuan at 2010-09-24 19:29:56:
ÄãÓÐfortran¼¸Ð´µÄ°¢