原始筛:
就是根据素数的定义:素数的因子只有1和它本身。
所以就只需要从2开始遍历到根号n,看定义是否成立。
若符合定义,则是素数,反之,就不是素数。代码如下:
void isprime1(int n) {
for(int i=2; i<=sqrt(n); i++)
if(n%i==0) return false;
return true;
}
埃氏筛:
埃氏筛又称为普通筛,其原理就是从2开始打表,若2为质数,则2的倍数必然不是质数,因此便让表中的4 6 8 10。。。。全部删除,最后剩下的就是所要的素数。其运行速度要远高于普通筛。
bool is_prime[MAX+5];//标记是否为素数
int prime[MAX/2];//存放素数
int cnt=0;//已存放的素数个数
void isprime2() {
memset(is_prime,true,sizeof is_prime);
for(int i=2; i<=MAX; i++) {
if(is_prime[i]) {
prime[cnt++]=i;
for(int j=i+i; j<=MAX; j+=i) //二次筛选法:i是素数,则i的倍数必不为素数,筛掉
is_prime[j]=false;
}
}
}
欧拉筛:
欧拉筛的原理就是利用:最小质因数 × 最大因数(非自己) = 这个合数。
只需利用最小质因子即可。原理和埃氏筛较为类似,但并不是单纯的一个数一个数的找倍数,其效率又远远高于埃氏筛。我认为其中也利用打表的思路。
#include <cstdio>
#include <cstring>
bool isPrime[100000010];
//isPrime[i] == 1表示:i是素数
int Prime[6000010], cnt = 0;
//Prime存质数
void GetPrime(int n)//筛到n
{
memset(isPrime, 1, sizeof(isPrime));
//以“每个数都是素数”为初始状态,逐个删去
isPrime[1] = 0;//1不是素数
for(int i = 2; i <= n; i++)
{
if(isPrime[i])//没筛掉
Prime[++cnt] = i; //i成为下一个素数
for(int j = 1; j <= cnt && i*Prime[j] <= n/*不超上限*/; j++)//注意j要<=cnt, 至少要把自身的平方给筛掉把
{
//从Prime[1],即最小质数2开始,逐个枚举已知的质数,并期望Prime[j]是(i*Prime[j])的最小质因数
//当然,i肯定比Prime[j]大,因为Prime[j]是在i之前得出的
isPrime[i*Prime[j]] = 0;
if(i % Prime[j] == 0)//i中也含有Prime[j]这个因子
break; //重要步骤。见原理
//若满足i%Prime[j]==0条件,就可以保证Prime[j]是i的因子,而欧拉筛是利用最小质因子,因而就不需要i了,所以break即可。
}
}
}
int main()
{
int n, q;
scanf("%d %d", &n, &q);
GetPrime(n);
while (q--)
{
int k;
scanf("%d", &k);
printf("%d\n", Prime[k]);
}
return 0;
}
最近看了一篇大神的题解,似乎发现了更高效简洁的筛法:
void init(int n){
ntp[1]=1;
for(int i=2;i<=n;++i){
if(!ntp[i]){
for(int j=i+i;j<=n;j+=i)ntp[j]=1;
}
}
}
ntp数组的值若为1则说明该数不是素数,否则是