让我们定义 dn 为:dn = pn+1 - pn,其中 pi 是第i个素数。显然有 d1=1 且对于n>1有 dn 是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数N (< 105),请计算不超过N的满足猜想的素数对的个数。
输入格式:每个测试输入包含1个测试用例,给出正整数N。
输出格式:每个测试用例的输出占一行,不超过N的满足猜想的素数对的个数。
题目来自:http://pat.zju.edu.cn/contests/pat-b-practise/1007
#include<iostream>
#include<cmath>
#include<vector>
using namespace std;
bool prime_check(int number) //检查数字是否为素数
{
int sqrt_number=sqrt(number);
for(int i=2;i<=sqrt_number;++i)
if(!(number%i)) return false;
return true;
}
int prime_pair(int number)
{
vector<int> prime_number;
for(int i=1;i<=number;++i) //记录素数
{
if(prime_check(i))
{
prime_number.push_back(i);
}
}
int number_prime=prime_number.size();
int number_prime_pair=0;
for(int j=0,k=1;k<number_prime;++j,++k) //判断是否满足素数对要求
{
if(prime_number[k]-prime_number[j]==2)
++number_prime_pair;
}
return number_prime_pair;
}
int main()
{
unsigned long int max_number;
cin>>max_number;
if(max_number<1 || max_number>100000)
return -1;
int prime_pair_number=prime_pair(max_number);
cout<<prime_pair_number<<endl;
return 0;
}
上面的代码中,使用下标操作符[]进行vector元素的存取,在我改成迭代器之后,开始代码并不能正常工作,究其原因,是因为vector迭代器不支持前向自增运算(++a),在我改成后向自增(a++)之后,代码便工作正常。
#include<iostream>
#include<cmath>
#include<vector>
using namespace std;
bool prime_check(int number) //检查数字是否为素数
{
int sqrt_number=sqrt(number);
for(int i=2;i<=sqrt_number;++i)
if(!(number%i)) return false;
return true;
}
int prime_pair(int number)
{
vector<int> prime_number;
for(int i=1;i<=number;++i) //记录素数
{
if(prime_check(i))
{
prime_number.push_back(i);
}
}
vector<int>::iterator first=prime_number.begin();
vector<int>::iterator second=first;
int number_prime_pair=0;
for(second++;second!=prime_number.end();first++,second++) //判断是否满足素数对要求
{
if((*second)-(*first)==2)
++number_prime_pair;
}
return number_prime_pair;
}
int main()
{
unsigned long int max_number;
cin>>max_number;
if(max_number<1 || max_number>100000)
return -1;
int prime_pair_number=prime_pair(max_number);
cout<<prime_pair_number<<endl;
return 0;
}