让我们定义 dn 为:dn = pn+1 - pn,其中 pi 是第i个素数。显然有 d1=1 且对于n>1有 dn 是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数N (< 105),请计算不超过N的满足猜想的素数对的个数。
输入格式:每个测试输入包含1个测试用例,给出正整数N。
输出格式:每个测试用例的输出占一行,不超过N的满足猜想的素数对的个数。
输入样例:20输出样例:
4
思路看起来很简单,但是实现的过程中遇到两个问题
1:运行超时,第一次的代码两个for循环中外层写的n,内层写的n/2+1,,改进后外层只考虑奇数,内层只考虑到根号j即可
2:出现段错误问题,,整了半天,发现是vector定义的时候出问题,于是指定了他的空间2000;段错误:
可能是数组越界,堆栈溢出(比如,递归调用层数太多)等情况引起
代码如下:
#include<iostream>
#include"math.h"
#include<vector>
#include<iterator>
using namespace std;
void findPrimNum(vector<int>& src,int n)
{
for (int j = 3;j <= n;j+=2)
{
int flag = 0;
for (int i = 2;i < sqrt(j)+1;i++)
{
if (j%i == 0)
{
flag = 1;
break;
}
}
if (flag == 0)
src.push_back(j);
}
}
void findResult(vector<int> src)
{
int counter = 0;
vector<int>::iterator it1=src.begin(), it2=src.begin();
for (it1 = src.begin(),it2=it1+1;it2 != src.end();it1++,it2++ )
{
if ((*it2 - *it1) == 2)
counter++;
}
cout << counter;
}
int main()
{
int n;
vector<int> src(2000);
cin >> n;
findPrimNum(src,n);
findResult(src);
int a;
cin >> a;
return 0;
}
补充知识:
好点的算法:
我们知道,我们的算法如果写成线性算法,也就是O(n),已经算是不错了,但是最好的是O(Log(n))的算法,这是一个对数级的算法,著名的二分取中(Binary Search)正是O(Log(n))的算法。通常来说,O(Log(n))的算法都是以排除法做为手段的。所以,找质数的算法完全可以采用排除法的方式。如下所示,这种算法的复杂度是O(n(log(logn)))。
示例:打印30以内的质数
一、初始化如下列表。
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
二、把第一个数(2)取出来,去掉所有可以被2整除的数。
2 3 5 7 9 11 13 15 17 19 21 23 25 27 29
三、取第二个数(3),去掉所有可以被 3整除的数。
2 3 5 7 11 13 17 19 23 25 29
四、取第三个数(5),因为4已经被去除了,再去掉所有可以被5整除的数。
2 3 5 7 11 13 17 19 23 29
接下来的数是7,但是7的平方是49,其大于了30,所以我们可以停止计算了。剩下的数就是所有的质数了。
拓展:
http://coolshell.cn/articles/3738.html
http://coolshell.cn/articles/2704.html
http://www.cnblogs.com/luluping/archive/2010/03/03/1677552.html