#include
#include
void main(void)
{
        int n,i=1,j,m,flag;
        for(n=3;1;n+=2)
        {
                flag=0;
                m=sqrt(n);
                for(j=2;j<=m;j++)
                {
                        if(n%j==0)
                                flag++;
                }
                if(flag==0) i++;
                if(i==10001) break;
        }
        printf("%d\n",n);
}

function euler7
    tic;
    n = 10001;
    primes = ones(n,1);
    primes(1:2) = [2 3];
    i = 3;
    m = 5;
    while i <= n
        idx = find(primes > sqrt(m), 1);
        if all(mod(m, primes(1:idx)) ~= 0)
            primes(i) = m;
            i = i + 1;
        end
        m = m + 2;
    end
    disp(num2str(m));
    toc
end

Program euler7
    Implicit None
    Integer, Parameter :: N = 10001
    Integer :: Primes(N)
    Integer :: I, M, IDX
    Real    :: StartTime, EndTime

    Call CPU_Time(StartTime)

    Primes(1:2) = (/2,3/)
    I = 3
    M = 5

    Do While(I <= N)
        If(All(Mod(M, Primes(1:I-1)) > 0)) Then
            Primes(I) = M
            I = I + 1
        EndIf
        M = M + 2
    EndDo
    Call CPU_Time(EndTime)
    Print *, M - 2
    Print *, EndTime-StartTime

End Program euler7

Originally posted by holmescn at 2011-05-14 13:01:05:
虽然使用了筛选法，但这个matlab程序还是需要7秒多。

[code]
function euler7
    tic;
    n = 10001;
    primes = ones(n,1);
    primes(1:2) = [2 3];
    i = 3;
    m = 5;
    while i <=  ...

Originally posted by holmescn at 2011-05-14 18:03:17:
这里的性能损失我想主要是来自一些重复（冗余）计算，比如不用计算每个数的余数，只要有一个能整除，那后面的就不用算了。这一点我不知道matlab是不是给优化掉了。

至于分段，可能是我没有太理解你的意思， ...