简单素数问题

原文链接：https://blog.csdn.net/weixin_43318827/article/details/105615753

素数的判断

判断一个数n是否为素数，就需要判断n是否能被2、3…、n-1中的一个数整除。只有都不能整除n，n才能判定为素数。

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

这样的判断方法没问题，时间复杂度为O(n)。但判定素数只是整个算法中的一部分，这时候O(n)的复杂度其实有点大了。
可以优化一下：

假设在2~n-1中存在n的约数，设为k，即n % k == 0，那么由k * (n / k) == n可知，n / k也是n的一个约数，且k和n / k中一定满足其中一个小于等于sqrt(n)、另一个大于等于sqrt(n)。

所以我们只需要判定n能否被2,3,…,⌊sqrt(n)⌋中的一个整除，即可判定n是否为素数，时间复杂度优化到了O(sqrt(n))。

bool isPrime(int n)
{
	if (n <= 1) return false;
	int sqr = (int)sqrt(1.0 * n);
	for (int i = 2; i <= sqr; i++)
	{
		if (n % i == 0)
			return false;
	}
	return true;
}

如果n没有接近int型变量的范围上界，那就可以这样写：

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

这样写当n接近int型变量的范围上界时导致会i * i溢出(不过n在10⁹以内都是安全的)，解决办法是将i定义为long long。不过更推荐开根号写法，这样更安全。

素数表的获取

可以由上面的方法打印1 ~ n范围内素数表的方法，即从1 ~ n枚举，判断每个数是否为素数，是则列入素数表。枚举时间复杂度O(n)，判断素数复杂度O(√n)，所以总复杂度O(n√n)，这个时间复杂度对n不超过10⁵是没问题的。

const int maxn = 100;
int prime[maxn], pNum = 0;
bool p[maxn] = {0};
void Find_prime()
{
	for (int i = 0; i < maxn; i++)
	{
		if (isPrime(i) == true)
		{
			prime[pNum++] == i;
			p[i] = true;
		}
	}
}

上述方法对n不超过10⁵是没问题的，但是如果需要更大范围的素数表就不够了。现在介绍一个更高效的算法，时间复杂度有Q(nloglogn)。
埃氏筛法的关键在一个“筛”字。算法从小到大枚举所有数，对每一个素数，筛去它的所有倍数，剩下的就是素数了。
举例：求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已经被筛去了，因此10不是素数

11没有被筛，所以11是素数，筛去11的倍数，但超过了15

12已经被筛去了，因此12不是素数

13没有被筛，所以13是素数，筛去13的倍数，但超过了15

14已经被筛去了，因此14不是素数

15已经被筛去了，因此15不是素数

至此已经得到了1~15以内的所有素数2、3、5、7、11、13
由上面的例子可以发现，从小到大到达某数a的时候，如果a没有被前面的步骤筛去，那么a一定是素数。因为如果a不是素数，那么一定有小于a的素因子，这样在之前的步骤中一定会被筛掉，所以如果当枚举到a时还没被筛掉，那么a一定是素数。

const int maxn = 101;
int prime[maxn], pNum = 0; //prime存放素数
bool p[maxn] = {0}; //false表示是素数，true表示不是
void Find_Prime()
{
	for (int i = 2; i < maxn; i++)
	{
		if (p[i] == false) //如果是素数
		{
			prime[pNum++] = i; //存入prime中
			for (int j = i + i; j < maxn; j += i) //然后筛去所有i的倍数
				p[j] = true;
		}
	}
}

要注意几点：

1不是素数(也不是合数)
素数表长至少要比n大1
Find_prime()里面for循环条件必须是i < maxn而不是i <= maxn
例题
PAT乙级1013 数素数