线性筛判断素数

线性筛存素数的数组不能开到1e7！！！

【原始】

时间复杂度O(n*sqrt(n))

bool isprime(int n)
{
    int i;
    for(i=2; i<=sqrt(n); i++)
        if(n%i==0)
            return false;
    return true;
}

【普通筛——埃拉托斯特尼(Eratosthenes)筛法】

时间复杂度O(nloglogn)

bool number[maxn+5];
void isprime()
{
    int i,j;
    memset(number,true,sizeof(number));
    for(i=2;i<=sqrt(N);i++)
    {
        if(number[i]==true)
        {
            for(j=i*i;j<=N;j+=i)
                number[j]=false;
            //二次筛选法:i是素数，则下一个起点是i*i,把后面的所有的i*i+2*n*i筛掉
        }
    }
}

【线性筛——欧拉Euler筛】

时间复杂度为O(n)

   prime[]数组中的素数是递增的,当i能整除prime[j]，那么i*prime[j+1]这个合数肯定被prime[j]乘以某个数筛掉。
   因为i中含有prime[j],prime[j]比prime[j+1]小，即i=k*prime[j]，那么i*prime[j+1]=(k*prime[j])*prime
   [j+1]=k’*prime[j]，接下去的素数同理。所以不用筛下去了。因此，在满足i%prime[j]==0这个条件之前以及第一次
   满足改条件时,prime[j]必定是prime[j]*i的最小因子。
bool number[maxn+5];
void isprime()
{
    int prime[maxn+5];
    int i,j,c=0;
    memset(number,true,sizeof(number));
    for(i=2;i<=maxn;i++)
    {
        if(number[i])
            prime[c++]=i;
        for(j=0;j<c&&prime[j]*i<=maxn;j++)
        {
            number[prime[j]*i]=false;
            if(i%prime[j]==0) //保证每个合数只会被它的最小质因数筛去，因此每个数只会被标记一次
                break;
        }
    }
}

例题：Goldbach`s Conjecture LightOJ - 1259 Runtime Error

参考博客：判断素数（一般筛到线性筛）