算法笔记（四）——欧拉筛法求素数

 先附代码 （来自http://www.cnblogs.com/A-S-KirigiriKyoko/articles/6034572.html）

/*求小于等于n的素数的个数*/
#include<stdio.h>
#include<string.h>
using namespace std;
int main()
{
    int n, cnt = 0;
    int prime[100001];//存素数 
    bool vis[100001];//保证不做素数的倍数 
    scanf("%d", &n);
    memset(vis, false, sizeof(vis));//初始化 
    memset(prime, 0, sizeof(prime));
    for(int i = 2; i <= n; i++)
    {
        if(!vis[i])//不是目前找到的素数的倍数 
        prime[cnt++] = i;//找到素数
        for(int j = 0; j<cnt && i*prime[j]<=n; j++)
        {
            vis[i*prime[j]] = true;//找到的素数的倍数不访问 
            if(i % prime[j] == 0) break;//关键！！！！ 
        }
    }
    printf("%d\n", cnt);
    return 0;
}

 其中最难理解的一步是 if(i%prime[j]==0) break; 以下是我看一些大牛的理解后的小结（有错请指出^_^）。

 如果i%prime[j]==0，那么i就可以看成prime[i]乘以一个数（记为n）；因为我们存的prime是从小到大存的，所以prime[j+1]>prime[j]；那么i*prime[j+1]就可以看成prime[j]*n*prime[j+1];那么i*prime[j+1]的最小质因子是prime[j]。加上if(i%prime[j]==0) break后我们就可以保证每一个数只被它的最小质因子给筛去,那么每一个数只会被筛一次，那么时间复杂度为O(n)。举个例子，6=2*3；当i=2时i%2=0，跳出循环，那它就不会被质因子3给筛去；当i=3时，6就被质因子2给筛去了。下图是n=20的模拟。（p[j]是代码中的prime[j]，tot是cnt）（本图来自https://www.cnblogs.com/grubbyskyer/p/3852421.html ）