[交流] Euler 工程 第五题：能被1到20所有的数都整除的最小正数 已有11人参与

        PROGRAM MAIN
        IMPLICIT NONE
        INTEGER N,I,J,K,M,FLAG
        INTEGER(8) SUM
        DIMENSION K(10)
        DO I=1,10
        K(I)=1
        END DO
        READ(*,*)N
        J=1
        DO I=2,N
                CALL ZS(I,FLAG)
                IF(FLAG.EQ.1)THEN
                        K(J)=I
                        J=J+1
                END IF
        END DO       
        DO I=1,J-1
                CALL KN(K(I),N)
        END DO
        SUM=1
        DO I=1,J-1
                SUM=SUM*K(I)
        END DO
        WRITE(*,*)SUM
        END

C        判断质数，FLAG=1为质数，否则不为质数。
        INTEGER        FUNCTION ZS(II,FLAG)
        INTEGER II,I,J,K,FLAG
        K=INT(II/2)
        FLAG=1
        DO I=2,K
        IF(MOD(II,I).EQ.0)THEN
                FLAG=0
                GOTO 10
        END IF
        END DO
10        RETURN
        END

C        判断某个质数的n次方在规定范围以内
        INTEGER FUNCTION KN(X,N)
        IMPLICIT NONE
        INTEGER X,Y,N,SUM
        Y=1
        SUM=X
        DO WHILE(SUM.LT.N)
                SUM=SUM*X
                Y=Y+1
C        WRITE(*,*)X,N
        END DO
        X=SUM/X
        END

%% evenly divided by 1:20
% Elapsed time is 173.211097 seconds.
% ans =
%    232792560
function result = euler5()
tic;
flag = 0;
result = 2520; % 能被1-20整除,肯定比能被1-10整除的2520大
while flag==0
    result = result+10; % 能被10整除,所以每次增加10
    flag = ~any(mod(result,2:20)); % 检测2-20,如果全部可以整除,改变flag结束循环
end
toc;
end

Originally posted by libralibra at 2011-05-12 16:36:07:
[code] %% evenly divided by 1:20
% Elapsed time is 173.211097 seconds.
% ans =
%    232792560
function result = euler5()
tic;
flag = 0;
result = 2520; % 能被1-20整除,肯定比能被1-10整除的2520 ...

#include
#include
int x(int i);
void main(void)
{
        int i;
        float dif;
        time_t start,end;

        time(&start);
        for(i=1;1;i++)
                if(x(i))
                {
                        printf("%d\n",i);
                        break;
                }
        time(&end);
        dif=difftime(end,start);
        printf("运算时间：%.1f秒\n",dif);
}
int x(int i)
{
        int j,k=0;
        for(j=2;j<21;j++)
                if(i%j==0)
                        k++;
        if(k==19)
                return 1;
        else return 0;
}