1007 素数对猜想 (20 分)
让我们定义dn为:dn=pn+1−pn,其中pi是第i个素数。显然有d1=1,且对于n>1有dn是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数N
(<105),请计算不超过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;
}