你必学的素数筛选法！PAT 乙级 1007（C++）

1007 素数对猜想 （20 分）

让我们定义d​n​​为：d​n​​=p​n+1​​−p​n​​，其中p​i​​是第i个素数。显然有d​1​​=1，且对于n>1有d​n​​是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。

现给定任意正整数N(<10​5​​)，请计算不超过N的满足猜想的素数对的个数。

输入格式:

输入在一行给出正整数N。

输出格式:

在一行中输出不超过N的满足猜想的素数对的个数。

输入样例:

20

输出样例:

4

---

 思路分析：

用数组记录下N之内的所有素数，再遍历数组，若数组中相邻数字之差为二，则计数器加一。

素数筛选法：

欧拉素数筛选法是必学算法之一，也称线性筛选法。其思路是：用布尔数组筛选素数，默认全部是素数。逐个遍历过去，若是素数则存入素数数组中，并将其倍数筛选出来，如果超过了n就退出；而每次循环就会将仍为false的数字加入素数数组，因为合数的因子一定比自己小，遍历因子时必然作为倍数出现过，所以当前布尔数组仍是false则说明是素数。我觉得这里的思想是计算机这门学科非常通用的思想，在动态规划、KMP中都有体现。我现在还讲不太清，大家结合代码体会吧嘻嘻嘻。

代码：

#include <iostream>
#include <string>
#include <cmath>
#include <stack>
#include <algorithm>
using namespace std;
bool vis[100000];//被筛出的数字为true；默认全是false
int Prime[100000];//存储素数
int cnt = 0;
void Shift(int n)//素数筛选函数
{
	int i, j;
	vis[1] = 1;//1既不是素数也不是合数
	for (i = 2;i <= n;i++)
	{
		if (vis[i] == false) Prime[++cnt] = i;//遍历因子时会把合数都筛去
		for (j = 1;j <= cnt&&Prime[j] * i <= n;j++)
		{
			vis[Prime[j] * i] = 1;//筛出
			if (i%Prime[j] == 0) break;
			//如果i是合数，则i及其所有倍数一定已经被筛去了
		}
	}
}
int main()
{
	int N, i, j,c=0;
	cin >> N;
	Shift(N);
	for (i = 2;i <= N;i++)
	{
		if (Prime[i] - Prime[i - 1] == 2)
			c++;
	}
	cout << c << endl;
	system("pause");
	return 0;

}