天天看点

素数筛法打表

素数打表比较快的方法是利用素数的倍数不是素数的性质来找的,算法将近为线性时间O(n)

2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30

第一遍先删去二的倍数

3 5 7 9 11 13 15 17 19 21 23 25 27 29

再删去3的倍数

2 3 5 7 11 13 17 19 23 29

再删5的倍数

2 3 5 7 11 13 17 19 23 29

因为7*7=49>30了所以我们到7既可以了

那么剩下的都是30以内的素数了

const int maxn=1e6+5;
bool vis[101000000];
int prime[maxn];
void init()
{
    memset(vis,true,sizeof(vis));
    int num=0;
    for(int i=2;i<maxn;i++)
    {
        if(vis[i]==true)
        {
            prime[num++]=i;
        }
        for(int j=0;(j<=num)&&(i*prime[j]<maxn);j++)
        {
            vis[i*prime[j]]=false;
            if(i%prime[j]==0) break;
        }
    }
}