【数论】素数（二）：素数筛法（埃式筛、欧拉筛、区间筛）

我的数论-素数部分博客共5part：
基本概念、性质、猜想、定理
素数筛法（埃式筛、欧拉筛、区间筛）
素数判断法（朴素法、模6法、Rabin-Miller及改进）
数的分解（Pollard-rho）
梅森素数（Lucas_Lehmer判定法）

素数筛法

埃式筛法

$O(nloglogn)$

bool b[N];
bool prime[N];
void solve(int x)
{
    for (int i = 0; i <= x; i++)
        prime[i] = true;
    b[0] = b[1] = 0;
    for (int i = 2; i <= x; i++)
        if (b[i])
            for (int j = i << 1; j <= x; j += i)
                prime[j] = false;
}

欧拉筛法

• $O(n)$

• 因为$i$ 可以写成$prime[j]\times k$，如果要继续筛$i\times prime[j+1]$ 就可以写成$prime[j]\times k\times prime[j+1]$ ，是$prime[j]$ 的整数倍，因为$prime[j+1]>prime[j]$ 在$i$ 继续增的过程中一定会增到$k\times prime[j+1]$ ，那个时候会先把$i\times prime[j+1]$ 筛去。

  因此跳出循环时$prime[j]$ 正是$prime[j]\times i$ 的最小质因子

typedef int ll; //原来用long long爆空间了
const int MAX_L = (1e8) + 2;
ll a, b;
bool check[MAX_L]; //检查某个数是否是素数，初始化置false
ll prime[MAX_L];
ll tot = 0;
void sieve(ll x)
{
    for (int i = 2; i <= x; i++)
    {
        if (!check[i])
            prime[tot++] = i; //i是素数，添加到prime中
        for (int j = 0; j < tot; j++)
        {
            if (prime[j] * i > x)
                break;
            check[prime[j] * i] = true;
            if (i % prime[j] == 0)
                break; //已被筛过，可以不往下筛了，简单说明见下：
        }
    }
}

区间筛素数

• $[a,b)$

typedef long long ll;
const int MAX_L = (1e8) + 2;
const int MAX_SQRTB = (1e4) + 2;
ll a, b;
bool isprime[MAX_L];           //[a,b)区间，用isprime[i-a]=true表示i为素数
bool isprime_small[MAX_SQRTB]; //[2,sqrt(b))区间
void segment_sieve(ll x, ll y)
{
    for (ll i = 0; (ll)i * i < y; i++)
        isprime_small[i] = true;
    for (ll i = 0; i < y - x; i++)
        isprime[i] = true;
    for (int i = 2; (ll)i * i < y; i++)
        if (isprime_small[i])
        {
            for (ll j = 2 * i; (ll)j * j < y; j += i)
                isprime_small[j] = false; //从2开始筛[2,sqrt(b))区间
            for (ll j = max((ll)2, (x + i - 1) / i) * i; j < y; j += i)
                //max找离x最近的i的倍数，如果没有就是2倍i
                isprime[j - a] = false;
        }
}