素数筛法 埃氏筛 线性筛

素数筛法？
呃~，太暴力的试1~sqrt（n）时候是因数就不说了。
这里介绍埃氏筛和线性筛。

埃氏筛

埃氏筛主要原理是：
每次搜索到一个质数将他的倍数标记为合数；
本质原理就是：一个合数一定表示为一个质数*一个数；

一些tips：
1、循环时从2-sqrt（n）就好，不然小心爆掉哦~

2、内层循环从i²开始：
假设x=i*k（k

#include<iostream>
#include<cstring>
#include<cstdio>
#define maxn 10010
using namespace std;

int n,tot,p[maxn];
bool check[maxn];

void prime()
{
    memset(check,0,sizeof(check));
    for(int i=2;i*i<=n;i++)//tip1：i=1-sqrt（n）
    {
        if(!check[i])
        {   
            p[++tot]=i;
            for(int j=i*i;j<=n;j+=i)check[j]=1;//tip2：j=i^2-n
        } 
    }
}

int main()
{
    n=1000;
    prime();
    for(int i=1;i<=tot;i++) 
        printf("%d ",p[i]);
}

然而，很显然埃氏筛会重复判断许多合数；
显然，这并不优。

线性筛

线性筛较埃氏筛最大提高之处就在于它对于每一个数只会筛一次；
所以复杂度为O（n），这也是它名字的由来。

线性筛原理：
对于每一个数与质数表中的数相乘并标为合数；
然而这样还是会重复……
所以对每一个i进行筛的时候如果i是质数的倍数，
那么就跳出循环。

此做法可以保证是当i=a时a*b时筛去；
i=a*b；
所以b是i的最小质因数；
随后就会跳出当前循环。
效率由此确定。

#include<iostream>
#include<cstring>
#include<cstdio>
#define maxn 100000
using namespace std;

int n,tot,p[maxn];
bool check[maxn];

void prime()
{
    check[1]=1;
    check[0]=1;
    for(int i=2;i<=n;i++)
    {
        if(!check[i])p[++tot]=i;
        for(int j=1;j<=tot&&i*p[j]<=n;j++)
        {
            check[i*p[j]]=1;
            if(!(i%p[j]))break;//*****关键
        }
    }
}

int main()
{
    n=1000;
    prime();
    for(int i=1;i<tot;i++) 
        printf("%d ",p[i]);
}

ps：别忘了0，1，n的问题~~~~（表示被坑过）。