从2021.9.5日起,每天写1~2篇博客
埃氏筛法核心就是筛掉每一个素数的倍数,比方说2能筛4,6,8;但是有可能有重复,比如6可以被2,3筛掉。因此——>重点说的欧拉筛就有用了。
void eratosthenes_sieve(int n)
{
totPrimes = 0;
memset(flag, 0, sizeof(flag));
int sqrtn = sqrt(n + 0.5);
for (int i = 2; i <= sqrtn; i++) {
if (!flag[i]) {
primes[totPrimes++] = i;
for (int j = i * i; j <= n; j += i) {
flag[j] = true;
}
}
}
for (int i = sqrtn + 1; i <= n; i++) {
if (!flag[i])
primes[++totPrimes] = i;
}
}//此为埃氏筛。
(借鉴自https://blog.csdn.net/baoli1008/article/details/50788512)
欧拉筛是一种线性筛,因为每个数只被他的最小质因子筛掉一次。
void euler_sieve(int n)
{
totPrimes = 0;
memset(flag, 0, sizeof(flag));
for (int i = 2; i <= n; i++) {
if (!flag[i])
primes[totPrimes++] = i;
for (int j = 0; i * primes[j] <= n; j++) {
flag[i*primes[j]] = true;
if (i % primes[j] == 0)
break;
}
}
}//同样引用该博客。
if语句就是说看之前是否被筛掉(flag[i]是否=true)如果没有被筛掉,一定是一个素数,容易看出在i=某个合数之前,该合数一定会被它的最小质因子筛掉,因此flag[i]一定是true。
面对内层循环,判断条件是i×prime[j]<=n,j是不可能超过totprimes的为啥不会这里小证一波:由唯一分解定理,i一定可以由若干个质数的乘机表示,因此(i%primes[j])一定会被筛出的。接下来进入循环内部,举一个例子,当i==4时,它前面已经积累了2,3这两个prime了,进入内层for循环,将flag[8]标记,flag[i*primes[0]]也就是flag[4X2]=1,这个8是被它的最小质因子2筛掉,下一个prime是3,3X4=12,但是这里不能继续循环了,12不能被筛掉!!!因为4%2也就是i%这个素数=0,我们筛掉合数的原则是被它的最小质因子筛掉,12不能被3筛掉!否则就不满足这个原则了。因为该素数是他的一个因子,(2是4因子)所以4×X=2×(2X),当i=2X(6)的时候再筛掉这个数(2)就能保证被它最小的质因子筛掉,在本例子中X=3。
还有i是大于所有之前的素数的,如本例(4>2)&&(4>3),这也保证了flag[4x3]=1,flag[3x4]=1这样的重复赋值(重复筛)
因此这种线性写法可以保证了每个合数都能被其最小质因子筛掉有且仅有一次。
想到别的再补充。