判断素数优化

判断素数优化

第一版

#include <iostream>
using namespace std;
int is_prime(int x)
{
	if (x < 2 || (x > 2 && x % 2 == 0)) return 0;
	for (int i = 3; i * i <= x; i += 2)
		if (x % i == 0) return 0;
	return 1;
}
int main(void)
{
	int n;
	cin >> n;
	if (is_prime(n)) cout << "Yes";
	else cout << "No";
	return 0;
}

第二版

#include <iostream>
using namespace std;
int is_prime(int x)
{
	if (x == 1) return 0;
	for (int i = 2; i * i <= x; ++i)
		if (x % i == 0) return 0;
	return 1;
}
int main(void)
{
	int n;
	cin >> n;
	if (is_prime(n)) cout << "Yes";
	else cout << "No";
	return 0;
}

第三版（最优算法）

#include <iostream>
using namespace std;
bool is_prime(int x)
{
	if (x == 1) return false;
	if (x == 2 || x == 3) return true;
	if (x % 6 != 1 && x % 6 != 5) return false;
	for (int i = 5; i <= x / i; i += 6)
		if (x % i == 0 || x % (i + 2) == 0) return false;
	return true;
} 
int main(void)
{
	int n;
	cin >> n; 
	if (is_prime(n)) cout << "Yes";
	else cout << "No";
	return 0;
}

证明

首先素数有这样一个规律：大于等于5的质数一定和6的倍数相邻。例如5和7，11和13等。

证明过程如下：

令x≥1；则大于等于5的自然数表示如下：…，6x-1，6x，6x+1，6x+2，6x+3，6x+4，…可以看到，不在6的倍数两侧的数：6x，2(3x+1)，3(2x+1)，2(3x+2)都不是素数。可能为素数的就只有：6x-1和6x+1