Acwing 868. 筛质数

对于任意一个数p而言，如果它没有被删掉，就意味着从2~p-1中的任何一个数都没有把p删掉，说明p不是2~p-1中任何一个数的倍数，说明2~p-1中不存在任何一个p的约数，所以说，p是一个质数。

#include <iostream>
#include <algorithm>

using namespace std;

const int N = 1e6 + 10;

bool st[N];
int prime[N];
int n, cnt;

void get_prime(int n)
{
    for (int i = 2; i <= n; ++ i)
    {
        if (!st[i]) prime[cnt ++] = i;
        for (int j = i + i; j <= n; j += i) st[j] = true;
    }
}

int main()
{
    cin >> n;
    
    get_prime(n);
    
    cout << cnt << endl;
    
    return 0;
}

 时间复杂度分析：

 

 进一步优化得到埃氏筛法（只需筛掉质数的倍数）：

#include <iostream>
#include <algorithm>

using namespace std;

const int N = 1e6 + 10;

bool st[N];
int prime[N];
int n, cnt;

void get_prime(int n)
{
    for (int i = 2; i <= n; ++ i)
    {
        if (!st[i]) 
        {
            prime[cnt ++] = i;
            for (int j = i + i; j <= n; j += i) st[j] = true;
        }
    }
}

int main()
{
    cin >> n;
    
    get_prime(n);
    
    cout << cnt << endl;
    
    return 0;
}

下面引入线性筛法：

假设 i 是合数 t 的最大因数，t显然可能不唯一（例如 30 和 45 最大因数都是 15）。但是仔细想一想，必然有t = i * p(p为小于 i 的质数) 。
p为什么比 i 小？因为 i 是 t 的最大因数。
为什么 p 一定是质数？因为如果 p 是合数，那么 i 就一定不是 t 的最大因数，因为 p可以再拆成若干素数相乘，这些素数再与 i 相乘会使因数更大。
既然如此，我们只需要把所有小于 i 的质数 p 都挨个乘一次好了。可是，真相真的是这样的嘛？
其实不是的，一不小心就忘记了最初的条件。我们要满足 i 是 t 的最大因数。如果 p 大于 i 的最小质因数，那 i 还是 t 的最大因数嘛？显然不是，任何一个合数都能唯一分解为有限个质数的乘积，除去这其中最小的质因数，其他的都乘起来就是最大因数 i 。所以我们不能让 p 大于 i 的最小质因数。

线性筛法题解转载自：https://www.acwing.com/solution/content/2559/

1、因为prime中素数是递增的，所以如果i%prime[j]!=0代表i的最小质因数还没有找到，即i的最小质因数大于prime[j]。也就是说prime[j]就是i*prime[j]的最小质因数，于是i*prime[j]被它的最小质因数筛掉了。
2、如果当i%prime[j]==0时，代表i的最小质因数是prime[j]，那么iprime[j+k]这个合数的最小质因数就不是prime[j+k]而是prime[j]了。所以iprime[j+k]应该被prime[j]筛掉，而不是后续的prime[j+k]。于是在此时break。
3、综上所述达到了每个数仅筛一次的效果，时间复杂度O ( n ) .

#include <iostream>
#include <algorithm>

using namespace std;

const int N = 1e6 + 10;

bool st[N];
int prime[N];
int n, cnt;

void get_prime(int n)
{
    for (int i = 2; i <= n; ++ i)
    {
        if (!st[i]) prime[cnt ++] = i;
        
        for (int j = 0; prime[j] <= n / i; ++ j)
        {
            st[prime[j] * i] = true;
            if (i % prime[j] == 0) break;   // prime[j]一定是i的最小质因子
        }
    }
}

int main()
{
    cin >> n;
    
    get_prime(n);
    
    cout << cnt << endl;
    
    return 0;
}

 

问题：for(int j=0;primes[j]<=n/i;j ++)为什么不是for(int j=0;primes[j]<=n/i && j <= cnt ;j ++)

答：其实在if(i%primes[j]==0)这一步已经帮我们避免了j>cnt的情况。事实上对于每一个数我们都是在用它的最小质因数去筛的