欧拉筛,也被称为线性筛,是一种高效的筛选质数的方法。相比传统的埃拉托斯特尼筛法(简称埃氏筛),欧拉筛在筛选质数的过程中具有更高的效率,主要体现在它能保证每个合数只被它的最小质因子筛除,从而降低了重复操作的次数。
欧拉筛的基本原理
-
初始化数组:定义一个数组来记录每个数是否为质数,一般初始状态为“是质数”。
-
筛选质数:从最小的质数开始,对每个质数的倍数进行标记,标记它们为“非质数”。不同于埃氏筛,欧拉筛在标记时有特殊的规则。
-
特殊的标记规则:对于每个合数,只用其最小质因子进行筛除。这意味着对于合数
x
,如果它有一个最小质因子p
,则只在p
的筛选轮次中将x
标记为非质数。 -
线性复杂度:由于每个合数只被其最小质因子筛除,每个数只被标记一次,因此整个筛选过程的时间复杂度为O(n)。
代码演示
#include<stdio.h>
#include<stdbool.h>
int main()
{
bool a[100001] = {1,1};//bool默认为0,这样设置确保0,1不参与运算
int prime[100001];//素数存到这个数组里面
int tmp = 0;
int i, j;
long long n;
scanf_s("%lld", &n);
for (i = 2; i <= n; i++)//从2开始,0,1不参与运算
{
if (a[i] == 0)//如果布尔值为0,就是素数
prime[++tmp] = i;
for (j = 1; j <= tmp; j++)
{
if (i * prime[j] > n)//如果大于输入的n,就跳出循环
break;
a[i * prime[j]] = 1;//素数的倍数都不是素数,将其布尔值设为1
if (i % prime[j] == 0)//如果 i 可以被 Prime[j] 整除,则停止内部循环,因为 i * Prime[j] 将会有更小的质因数。
break;
}
}
return 0;
}
欧拉筛的重点
if(i % prime[j] == 0)break;
当 i % Prime[j] == 0
(即 i
能被 Prime[j]
整除),i
的最小质因子肯定不大于 Prime[j]
。此时,根据算法的逻辑,i
的所有倍数都已经在之前的迭代中被 Prime[j]
(或更小的质数)筛选过了。因为 Prime[j]
是 i
的因子,所以 i * Prime[j]
实际上是 Prime[j]
的某个倍数。这个倍数一定已经在处理 Prime[j]
的那一轮筛选中被标记过了,因为我们是按照顺序一个接一个地标记 Prime[j]
的倍数。因此,没有必要再次对 i * Prime[j]
进行标记,这也意味着内部循环可以停止。
这样做可以避免冗余的标记操作,因为任何数的倍数不会被它的因子的倍数多次标记。这是欧拉筛算法保持线性时间复杂度的关键步骤之一。