欧拉筛法的算法说明

• 算法要求

1.  算法的目的是筛掉所有的合数。
2.  每个合数都是至少两个（合数可以相同）质因数的乘积。
3.  我们希望以O(n)时间完成任务，这就要求每个合数只被筛掉一次（以下会证明其充分性）。
4.  为了保证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)。