今天看的欧拉筛,好家伙没给我看吐了,看的直起生理反应,如今总算是弄懂了,赶紧记录
先上代码模板
int primes[N], cnt; // primes[]存储所有素数
bool st[N]; // st[x]存储x是否被筛掉
void get_primes(int n)
{
for (int i = 2; i <= n; i ++ )
{
if (!st[i]) primes[cnt ++ ] = i;
//如果这个数没有被标记过,代表它不是某个数的倍数,则为质数
for (int j = 0; primes[j] <= n / i; j ++ )
{//i*primes[j]需要小于要求的最大的n,也可以写成i*primes[j]<=n
st[primes[j] * i] = true;
if (i % primes[j] == 0) break;
}
}
}
引用其他博客的部分解释
对于代码中i和primes[j]的理解
还有一个打表展示的
欧拉筛法(线性筛)的学习理解
尽管上面已经解释的很详细,但我还是久久不能理解,思索半天写了几个例子才懂了,记录下我自己的部分理解。
p,i,t 的含义详见第一篇引用博客
其中i是t的最大质因数,而p不仅是t的最小质因数,
而且是i的最小质因数!!!很重要
我觉得欧拉筛比埃氏筛好就是因为这个
为什么呢?
如果看不懂博客中的解释,可以举几个例子
当t=6,p=2,12=6*2。此时满足6%2==0,即2是12和6的最小质因子,
此时如果不break
t=6,p=3,18=6*3。此时不满足6为18的最大质因数,也不满足3为18和6的最小质因数。
为什么会这样呢,因为你没break,导致本应该是当t=9,p=2,18=9*2
被剔除的18被提早剔除,这样18就会被剔除两次,欧拉筛失效。
但是没break为什么会出现这种现象呢,还是拿这个例子来说,
当没有break,说明6可以表示成一个之前本应该break的最小质因子2和
另一个数也就是3的乘积,所以6=2*3,此时2才是最小质因子,
而对应的最大质因子也由6变成了3*3=9。被筛掉的是2*9=18。
与我们想的不符啊!!我们想的是6是被筛掉的数的最大质因子,所以这不就错了。
总结,primes[j]是i和primes[j]*i的最小质因子。而枚举的i我们把它当作某一个合数的最大质因子。
而j是从小到大进行枚举,所以第一次break掉的数就是i和某个合数的最小质因子,其余的小于break的,则仅仅是应该被筛掉的数的最小质因子。