线性筛法求素数

和朴素的筛法求素数不同，线性筛法可以省去许多重复的删除

 

例如12：它可以被2,6；3,4筛掉，被筛了2次；

所以新的算法是在处理数i时，枚举i之前的所有质数并用i*质数进行删除，直到第一个质数k能整除i就在这之后停止，因为接下来的所有质数所产生的需要被筛的数num=p*i都可以在后面的枚举到(i/k*p)时被筛掉：(i/k*p)*k（显然可以保证在k之前没有一个数能整除(i/k*p)），所以这样就能保证所有合数只被筛掉了一遍

void init(int k)
{
    memset(bl,1,sizeof(bl));
    bl[1]=0;
    for(int i=2;i<=k;i++)
    {
        if(bl[i]) prime[++tot]=i;
        for(int j=1;j<=tot;j++)
        {
            if(i*prime[j]>k) break;
            bl[i*prime[j]]=0;
            if(i%prime[j]==0) break;
        }
    }
}