素数筛选法
素数筛选法的目的就是筛选出在某一区间[m,n)内的所有素数,常见的办法有如下几种。
1.朴素的筛选法
朴素的筛选法思路很简单,先写一个判断是否是素数的函数isPrime(),然后从2到n分别调用isPrime()函数来检查。检查是否是质数的算法是核心,使用从2到n的开根的数作为除数(为什么是开根数呢?因为如果当前筛选到p,则说明1*p,2*p.....(p-1)*p已经在筛选1,2,3……p-1的时候被处理过了,现在我们只需要接着处理p*p即可,即要求p*p <=n)。复杂度大概是O(n*log(n)).虽然复杂度看起来比较和谐,但是显然直觉告诉我们调用函数次数过多。
bool isPrime(int n)
{
//if we process p turn,then ,2*p,3*p,4*p,....(p-1)*p was processed
//so we begin with p*p < n;
int i,sqr = sqrt((double)n);
for(i=2; i<=sqr; ++i)
{
if ( (n%i) == 0 )
return false;
}
return true;
}
接下来如果需要找到[m,n)的素数,只需要在[m,n)上调用isPrime()函数n-m次即可。
2.比朴素法更高效的筛选法
我们知道要求某一区间内的素数,只需要将这一区间内的合数去除,即筛除即可,这种办法利用的就是这种思想。
例如:2,3,4,5,6,7,8,9,10,11,12,13,14,15.......
第一遍筛掉2的倍数:剩下2,3,5,7,9,11,13,15.....
第二遍筛掉3的倍数:剩下2,3,5,7,11,13,....
第三遍筛掉5的倍数:(为什么是5而不是4,因为4已经被2的倍数筛掉了,再筛已经无意义)