质数筛选:
1.朴素筛法:
直接把2~n-1中质数和合数的倍数都筛一遍,其代码如下所示:
int primes[N],cnt=0;
bool st[N]={false};
void get_prime(int n)
{
for(int i=2;i<=n;i++)
{
if(!st[i])
primes[cnt++]=n;
for(int j=i+i;j<=n;j+=i)
st[j]=true;
}
}
时间复杂度为O(nlogn);
2.埃拉托色尼筛法(the Sieve of Eratosthenes):
首先说一下它的概念:
埃拉托斯特尼筛法,简称埃氏筛或爱氏筛,是一种由希腊数学家埃拉托斯特尼所提出的一种简单检定素数的算法。要得到自然数n以内的全部素数,必须把不大于根号n的所有素数的倍数剔除,剩下的就是质数;
要得到自然数n以内的全部素数,必须把不大于的所有素数的倍数剔除,剩下的就是素数。给出要筛数值的范围n,找出以内的素数。先用2去筛,即把2留下,把2的倍数剔除掉;再用下一个质数,也就是3筛,把3留下,把3的倍数剔除掉;接下去用下一个质数5筛,把5留下,把5的倍数剔除掉
其代码如下所示:
int primes[N],cnt=0;
bool st[N]={false};
void get_prime(int n)
{
for(int i=2;i<=n;i++)
{
if(!st[i])
{
primes[cnt++]=n;
for(int j=i+i;j<=n;j+=i)
st[j]=true;
}
}
}
埃氏筛法是由朴素筛法改进的,时间复杂度为O(nloglogn)。
3.线性筛法(欧拉筛法):
欧拉筛法的基本思想 :在埃氏筛法的基础上,让每个合数只被它的最小质因子筛选一次,以达到不重复的目的。其代码如下所示:
int primes[N],cnt=0;
bool st[N]={false};
//算法核心:x仅会被其最小质因子筛去
void get_prime(int x) {
for(int i = 2; i <= x; i++) {
if(!st[i]) prime[cnt++] = i;
for(int j = 0; prime[j] <= x / i; j++) {
//对于任意一个合数x,假设pj为x最小质因子,当i<x/pj时,一定会被筛掉
st[prime[j]*i] = true;
if(i % prime[j] == 0) break;
/*pj是prime[j]的缩写
1.i%pj == 0, pj定为i最小质因子,pj也定为pj*i最小质因子
2.i%pj != 0, pj定小于i的所有质因子,所以pj也为pj*i最小质因子
*/
}
}
}
它的时间复杂度是O(n)。
一般来说用线性筛会比较好;