線性篩法求素數

 

這個是經典的Eraosthenes篩法：

for (int i = 2; i * i < N; i++)
{
    if (tag[i]) continue;
    for (int j = i; j * j < N; j++)
        tag[i*j] = 1;
}
for (int i = 2; i < N; i++)
    if (!tag[i])
        prime[tol++] = i;

 

但是Eraosthenes篩法的速度並不快，原因在於對於一個合數，這種方法會重復的標記。一種線性篩素數的方法有效的解決了這一點，代碼如下：

 

void get_prime()
{
    int cnt = 0;
    for (int i = 2; i < N; i++)
    {
        if (!tag[i])    p[cnt++] = i;
        for (int j = 0; j < cnt && p[j] * i < N; j++)
        {
            tag[i*p[j]] = 1;
            if (i % p[j] == 0)
                break;
        }
    }
}