素数筛法详解

在大于1的自然数中，除了1和它本身以外不再有其他因数的自然数称为素数，也叫质数。

注意

1.算术基本定理：任何一个大于1的自然数 n,如果n不为质数，那么n可以唯一分解成有限个质数的乘积。
2. 若一个数 n 可以进行因数分解，则得到的两个数一定是有一个>=sqrt(n),另一个<=sqrt(n)。

以下介绍3种常见的素数筛法。

朴素筛法

此方法一般用于判断单个数是否为素数。时间复杂度为O(√n) 。

现设要判断n是否为素数，

1. 枚举从2到√n的所有自然数，若这些数都不是n的因数，则n为质数。

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>

int visit[100000],n,sum;//visit[i]=1说明i不是素数
int prime_permutation[100000];//数组prime_permutation存储得到的素数

//朴素筛法
int simple_sieve(int n)
{
    int i,m;
    m=(int)sqrt(n)+1;
    for(i=2;i<=m;i++)
    {
        if(n%i==0)
        {
            return 0;//n不是质数
        }
    }
    return 1;//n为质数
}

埃拉托斯特尼筛法

时间复杂度：O(nloglogn)。此法一般用于质数的区间筛选，即选出区间（2，n）中的质数。

1. 由当前已经找到的素数，从后面的数中筛去当前素数的倍数。

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>

int visit[100000],n,sum;//visit[i]=1说明i不是素数
int prime_permutation[100000];//数组prime_permutation存储得到的素数

//埃氏筛法
int eratosthenes_sieve(int n)//n为右界
{
    int i,j;
    sum=0;
    memset(visit,0,sizeof(visit));//先将标记数组清零
    visit[0]=1;
    visit[1]=1;//0和1不是素数
    for(i=2;i<=n;i++)
    {
        if(!visit[i])//若i为质数
        {
            prime_permutation[sum]=i;//将i存入数组
            sum++;
            for(j=2*i;j<=n;j+=i)//质数i的整数倍一定不是质数
            {
                visit[j]=1;
            }
        }
    }
    return 0;
}

缺点

某些合数会被多次筛除，从而导致运算重复。

例如对于合数6，（设右界n=7）

• 当 i =2时，在循环for(j=2*i;j<=n;j+=i)中，会执行一次visit[6]=1;。
• 当 i =3时，在循环for(j=2*i;j<=n;j+=i)中，又会执行一次visit[6]=1;。

欧拉筛法

时间复杂度：O（n）。此法一般用于质数的区间筛选，即选出区间（2，n）中的质数。

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>

int visit[100000],n,sum;//visit[i]=1说明i不是素数
int prime_permutation[100000];//数组prime_permutation存储得到的素数

//欧拉筛法
int euler_sieve(int n)//n为右界
{
    int i,j;
    sum=0;
    memset(visit,0,sizeof(visit));//先将标记数组清零
    visit[0]=1;
    visit[1]=1;//0和1不是素数
    for(i=2;i<=n;i++)
    {
        if(!visit[i])//若i为质数
        {
            prime_permutation[sum]=i;//将素数i存入数组
            sum++;
        }
        for(j=0;j<sum&&(i*prime_permutation[j]<=n);j++)
        {
            visit[i*prime_permutation[j]]=1;//先将合数 i*prime_permutation[j] 筛除掉
            if(i%prime_permutation[j]==0)//若prime_permutation[j]为i的因数,可设i=m*prime_permutation[j]
            {//则合数i*prime_permutation[j+1]=m*prime_permutation[j]*prime_permutation[j+1]
                break;//让该合数在i=m*prime_permutation[j+1]时被筛掉即可,所以此时跳出循环
            }
        }
    }
    return 0;
}

欧拉筛法是在埃氏筛法的基础上进行了优化。它的核心思想是保证每个合数仅被它最小的质因子筛除一次。这里的质因子是指存储在数组prime_permutation中的质数。

着重注意下面这段代码中的注释：

        for(j=0;j<sum&&(i*prime_permutation[j]<=n);j++)
        {
            visit[i*prime_permutation[j]]=1;//先将合数 i*prime_permutation[j] 筛除掉
            if(i%prime_permutation[j]==0)//若prime_permutation[j]为i的因数,可设i=m*prime_permutation[j]
            {//则合数i*prime_permutation[j+1]=m*prime_permutation[j]*prime_permutation[j+1]
                break;//让该合数在i=m*prime_permutation[j+1]时被筛掉即可,所以此时跳出循环
            }
        }

欧拉筛法的难点就在于如何理解if(i%prime_permutation[j]==0)这条判断语句的作用。

• 若i%prime_permutation[j]==0成立，则说明prime_permutation[j]为 i 的质因数，可设 i =m×prime_permutation[ j ]，则合数 i ×prime_permutation[ j+1]=m×prime_permutation[ j ]×prime_permutation[ j+1]。
• 又因为数组prime_permutation为递增数列，即▽j>0，prime_permutation[ j +1]>prime_permutation[ j ]。所以一旦i%prime_permutation[j]==0成立，即可说明prime_permutation[j]为 i 的最小质因数，则prime_permutation[ j ]也是合数 i ×prime_permutation[ j+1]的最小质因数。
• 根据每个合数仅被它最小的质因子筛除一次的原则，合数 i ×prime_permutation[ j+1]只能在i==m*prime_permutation[j+1]时，在循环for(j=0;j<sum&&(i*prime_permutation[j]<=n);j++)中，通过执行visit[i*prime_permutation[j]]=1语句被筛除。
• 因为m*prime_permutation[j+1]>m*prime_permutation[j]，所以合数 i ×prime_permutation[ j+1]不能在i==m×prime_permutation[j]时被筛除。
• 例如，讨论以下i=6时的情况。当i=6时，此时的数组prime_permutation已经存储的质数有2、3、5。显然，在进行for(j=0;j<sum&&(i*prime_permutation[j]<=n);j++)循环时，会有6%prime_permutation[0]即6%2==0。若此时不跳出该for循环，则在下一步会有visit[6*3]=1,也就是说会筛除合数18。18=2*9，它的最小质因数为2，也就是说，18应该在j==0即prime_permutation[j]==2时被筛除，此时应有i==9。问题是如前所述，若不及时跳出for循环，在i=6时18已被筛除，矛盾。

集成代码

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>

int visit[100000],n,sum;//visit[i]=1说明i不是素数
int prime_permutation[100000];//数组prime_permutation存储得到的素数

//朴素筛法
int simple_sieve(int n)
{
    int i,m;
    m=(int)sqrt(n)+1;
    for(i=2;i<=m;i++)
    {
        if(n%i==0)
        {
            return 0;
        }
    }
    return 1;
}

//埃氏筛法
int eratosthenes_sieve(int n)//n为右界
{
    int i,j;
    sum=0;
    memset(visit,0,sizeof(visit));//先将标记数组清零
    visit[0]=1;
    visit[1]=1;//0和1不是素数
    for(i=2;i<=n;i++)
    {
        if(!visit[i])//若i为质数
        {
            prime_permutation[sum]=i;//将i存入数组
            sum++;
            for(j=2*i;j<=n;j+=i)//质数i的整数倍一定不是质数
            {
                visit[j]=1;
            }
        }
    }
    return 0;
}

//欧拉筛法
int euler_sieve(int n)//n为右界
{
    int i,j;
    sum=0;
    memset(visit,0,sizeof(visit));//先将标记数组清零
    visit[0]=1;
    visit[1]=1;//0和1不是素数
    for(i=2;i<=n;i++)
    {
        if(!visit[i])//若i为质数
        {
            prime_permutation[sum]=i;
            sum++;
        }
        for(j=0;j<sum&&(i*prime_permutation[j]<=n);j++)
        {
            visit[i*prime_permutation[j]]=1;//先将合数 i*prime_permutation[j] 筛除掉
            if(i%prime_permutation[j]==0)//若prime_permutation[j]为i的因数,可设i=m*prime_permutation[j]
            {//则合数i*prime_permutation[j+1]=m*prime_permutation[j]*prime_permutation[j+1]
                break;//让该合数在i=m*prime_permutation[j+1]时被筛掉即可
            }
        }
    }
    return 0;
}

int main()
{
    int i;
    //输入数据
    scanf("%d",&n);
    //使用埃氏筛法筛选素数并输出结果
    eratosthenes_sieve(n);
    for(i=0;i<sum;i++)
    {
        printf("%d ",prime_permutation[i]);
    }
    memset(prime_permutation,0,sizeof(prime_permutation));
    printf("\n");
    //使用欧拉筛法筛选素数并输出结果
    euler_sieve(n);
    for(i=0;i<sum;i++)
    {
        printf("%d ",prime_permutation[i]);
    }
    return 0;
}