[模板]素数筛选

这套素数模板包含两个内容：

1、使用notprime[]数组构成一张布尔类型的表，用O(1)的时间复杂度判断一个数是不是素数。

2、素数筛选，筛选出小于MAXN的素数。 prime[0] 存的是素数的个数。

两种素数判断法，在MAXN 小于1000w时都能轻松搞定。

1、

<pre name="code" class="cpp">/* 素数筛选
 * 判断小于MAXN的数是不是素数。
 * notprime是一张表，为false的是素数，
 * 为true的不是素数。
 */<pre name="code" class="cpp">const int MAXN=100010;
bool notprime[MAXN];
void init()
{
    memset(notprime,false,sizeof(notprime));
    notprime[0] = notprime[1] = true;
    for(int i = 2; i<MAXN; i++) if(!notprime[i])
        {
            if (i>MAXN/i) continue;
            for(int j = i*i; j<MAXN; j+=i)
            {
                notprime[j] = true;
            }
        }
}

2、

const int MAXN=100010;
//素数筛选，筛选出小于MAXN的素数。 prime[0] 存的是素数的个数。
int prime[MAXN+1];
void getprime()
{
    memset(prime,0,sizeof(prime));
    for(int i = 2; i < MAXN;i++)
    {
        if (!prime[i]) prime[++prime[0]] = i;
        for(int j = 1; j<= prime[0] &&prime[j] <= MAXN/i;j++)
        {
            prime[prime[j]*i] = 1;
            if (i%prime[j] == 0) break;
        }
    }
=
}