素数筛算法与素数打表[数学]

素数，即质数，是很常见的问题，如何判断一个数是不是质数，如何求出一个区间内的所有质数，是很常见的问题。

因数只有1和它本身的数是质数。其中1不属于质数和合数。根据定义，我们就可以直接判断一个数字n是否是素数。从2开始，迭代到sqrt(n)。

这样得到求一个数是不是质数的时间复杂度是O(n^1/2)，得到一个1到n区间内所有质数的复杂度就是O(n^3/2)。

在求解问题时，为了通过预运算，简化时间复杂度，我们常通过打表的方式，将数值预先计算出来，然后直接存储到数组，素数表就是常用的打表方法。

所以我们关心如何求得n范围内所有的质数。

素数筛：

先把n个自然数按次序排列起来。1不是质数，也不是合数，要划去。第二个数2是质数留下来，而把2后面所有能被2整除的数都划去。2后面第一个没划去的数是3，把3留下，再把3后面所有能被3整除的数都划去。3后面第一个没划去的数是5，把5留下，再把5后面所有能被5整除的数都划去。这样一直做下去，就会把不超过N的全部合数都筛掉，留下的就是不超过N的全部质数。因为希腊人是把数写在涂腊的板上，每要划去一个数，就在上面记以小点，寻求质数的工作完毕后，这许多小点就像一个筛子，所以就把埃拉托斯特尼的方法叫做“埃拉托斯特尼筛法”，简称“筛法”。

其中存在一些点：

#define Max 1000000  
bool prime[Max];  
void IsPrime(){  
     prime[0]=prime[1]=0;prime[2]=1;  
     for(int i=3;i<max;i++)  
        prime[i]=i%2==0?0:1;  //首先将2的倍数全部筛掉
     int t=(int)sqrt(Max);  
     for(int i=3;i<=t;i++)  
       if(prime[i])  //没有被筛掉,即i是质数
         for(int j=i*i;j<Max;j+=2*i)//一个数j如果能表示成一个数的倍数，则一定不是质数，j+=2*i保证了不是偶数，省去了重复计算。  
            prime[j]=0;  
}  

相比上面，有更快速的实现方法：（避免了重复筛同一个数）

//prime[]存储了小于n的素数,下标从1开始  
//visit[]数组标记了是否为素数，bool类型  
//返回值num表示小于n的素数的个数  
ll prime[maxn];  
bool visit[maxn];  
ll init_prim(ll n)  
{  
    memset(visit, true, sizeof(visit));  
    ll num = 0;  
    for (ll i = 2; i <= n; ++i)  
    {  
        if (visit[i] == true)  
        {  
        num++;  
        prime[num] = i;  
        }  
        for (ll j = 1; ((j <= num) && (i * prime[j] <= n)); ++j)  
        {  
            visit[i * prime[j]] = false;  
            if (i % prime[j] == 0) break; //点睛之笔  
        }  
    }  
    return num;//返回质数个数  
}  

通过整体一次遍历，利用已经判断出来的素数进行筛除，通过合数判断进行了剪枝重复。实现了几乎线性的素数筛法。