- 算法要求
- 算法的目的是筛掉所有的合数。
- 每个合数都是至少两个(合数可以相同)质因数的乘积。
- 我们希望以O(n)时间完成任务,这就要求每个合数只被筛掉一次(以下会证明其充分性)。
- 为了保证3,欧拉筛法的做法是让每个合数只被其最小的质因数筛去。
接着,我们看欧拉筛法的核心代码:
isPrime[i*primes[j]] = false;
if(i%primes[j] == 0) break;
显然第一行是筛选合数的操作。它筛掉的合数格式是 i*primes[j],即一个质数系[2,n]内的任选数。下面我们研究为什么这样写可以使得primes[j]是最小的。
假设i=4,primes[j]=2。2*4=8(为了符合习惯,调换了顺序,与程序中相反)被筛去,4是2的倍数,因此结束循环。下一个应被筛掉的合数是3*4=12,我们研究为什么12没有在i=4时被筛去。
由
和要求3,我们希望12是被2*6筛去,而不是3*4。注意到4=2*2,说明4可以被小于3的质数2整除,因此我们可以断定:3*4不是我们要找的分解形式。这正对应着源代码第二行。而为何12会被2*6筛去,没有被6*2或4*3筛去,原理较简单,留给读者思考。
我们已经说明了所有合数都只会被筛掉一次。下面证明该算法的时间复杂度是O(n)。
算法中有两重循环,执行的总操作数=外层循环次数+内层循环总次数。外层循环共执行n-2次;源代码第一行被执行的次数为[2,n]中的偶数个数,该数显然小于n,因此内层循环的总次数小于n,总操作数=O(2n)=O(n)。