素数的判定

素数判定的演化过程：

判定一个数字p是不是质数：拿1 ~p 的所有数字来试除。

bool divisibiity_test(int n)
{
    for (int i=2; i<n; i++)
        if (n % i == 0)
            return false;   // 不是质数
    return true;    // 是质数

}

改进

bool divisibiity_test(int n)
{
//一个数至少有一个比sprt(n)小的质因数。
      for (int i=2; i<=sqrt(n); i++)
        if (n % i == 0)
            return false;
    return true;
}
时间复杂度：O(sqrtN) 

开一个数组来标记范围内的书是否是素数

bool prime[20000000];
void eratosthenes()
{
    memset(prime,1,sizeof(prime));
    prime[0] = 0;  prime[1] = 0;  // 0 和 1 不是质数
    for (int i=2; i<20000000; i++)
        if (prime[i]==1) // 删掉质数i的倍数 
           for (int j=i+i; j<20000000; j+=i)
                prime[j] = 0;
}

注意：删掉质数 i 的倍数时，早已删掉1 倍~ i-1 倍之间的合数了，所以直接从i倍开始删除。

代码可优化如下：
void eratosthenes(){
    memset(prime,1,sizeof(prime));
    prime[0] = 0;prime[1] = 0;
    for (int i=2; i<20000000; i++)
        if (prime[i]==1)
      //删掉i的倍数从i倍开始。
            for (int j=i*i; j<20000000; j+=i)
                prime[j] = 0;
}