《算法笔记》 5.4.2 节 —— 素数表的获取

《算法笔记》 5.4.2 节 —— 素数表的获取

一般方法的分析：

根据一般的判断素数方法，我们已经可以从 $1$ ~ $n$ 中打印素数表，即从$1$ ~ $n$ 进行枚举，判断每个数

是否为素数，如果是素数就加入素数表，这种方法的时间复杂度是 $O(n)$ 而判断素数的时间复杂度为

$O(\sqrt{n})$ ，因此总时间复杂度为 $O(n\sqrt{n})$ 这个复杂度对于 n 不超过 $10^5$的大小是没有问题的。

例（求解100以内所有的素数）：

#include <cstdio>
#include <cmath>
int IsPrime(int n)
{
    if (n <= 1) return 0;
    for (int i = 2; i < sqrt(n); ++i)
    {
        if (n % i == 0) return 0;
    }
    return 1;
}
int prime[101], pNum = 0;
bool p[101] = { 0 };
void Find_Prime()
{
    for (int i = 1; i < 101; ++i)
    {
        if (IsPrime(i))
        {
            prime[pNum++] = i;
            p[i] = true;
        }
    }
}
int main()
{
    Find_Prime();
    for (int i = 0; i < pNum; ++i)
    {
        printf("%d ", prime[i]);
    }
}

---

改进方法：

上面的算法如果在出现 n > $10^5$的情况时将显得力不从心。

这时我们可以考虑比较简单和容易理解的**“埃氏算法”**进行求解，其时间复杂度为 $O(n\log \log n)$

即Eratosthenes筛法。

算法从小到大枚举所有数，对每一个素数，筛去它的所有倍数，剩下的就都是素数了。

例：求 $1$~ $15$ 之间所有的素数，思路如下

①2是素数（唯一一个需要事先确定的素数），因此筛去 2 的倍数：

​ 2 3 4 5 6 7 8 9 10 11 12 13 14 15

②3 没有被前面的步骤筛去，因此 3 是素数，筛去所有 3 的倍数：

​ 2 3 4 5 6 7 8 9 10 11 12 13 14 15

③ 4 已经在①中被筛去，因此 4 不是素数。

④ 5 没有被前面的步骤筛去，因此 5 是素数，筛去所有 5 的倍数：

​ 2 3 4 5 6 7 8 9 10 11 12 13 14 15

⑤ 6 已经在①中被筛去，因此 6 不是素数。

⑥ 7 没有被前面的步骤筛去，因此 7 是素数，筛去所有 7 的倍数：

​ 2 3 4 5 6 7 8 9 10 11 12 13 14 15

⑦ 8 已经在①中被筛去，因此 8 不是素数。

⑧ 9 已经在②中被筛去，因此 9 不是素数。

⑨ 10 已经在①中被筛去，因此 9 不是素数。

⑩ 11 没有被前面的步骤筛去，因此 11 是素数，筛去所有 11 的倍数：

​ 2 3 4 5 6 7 8 9 10 11 12 13 14 15

⑪ 12 已经在①中被筛去，因此 12 不是素数。

⑫ 13 没有被前面的步骤筛去，因此 13 是素数，筛去所有 13 的倍数：

​ 2 3 4 5 6 7 8 9 10 11 12 13 14 15

⑬ 14 已经在①中被筛去，因此 14 不是素数。

⑭ 15 已经在②中被筛去，因此 15 不是素数。

至此 $1$~ $15$ 之间所有的素数已全部得到，即 $2,3,5,7,11,13$。

---

改进代码：

由上面的例子可以知道：当大小达到某数 a 时，如果 a 没有被前面的步骤消去，那么 a 一定是素数。

而“筛”这个动作，可以由一个布尔型数组来标记，若 a 被筛掉，则 p[a] 为 true ，否则 p[a] 为 false。

在程序开始，可以把数组的元素初始化为 false。

以下是求100以内所有素数的程序：

#include <cstdio>
const int maxn = 101;
int prime[maxn], pNum = 0;
bool p[maxn] = { 0 };
void Find_Prime()
{
	for (int i = 2; i < maxn; ++i)
	{
		if (p[i] == false)
		{
			prime[pNum++] = i;
			for (int j = i + i; j < maxn; j += i)
			{
				p[j] = true;
			}
		}
	}
}
int main()
{
	Find_Prime();
	for (int i = 0; i < pNum; ++i)
	{
		printf("%d ",prime[i]);
	}
}