PAT1007. 素数对猜想

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

现给定任意正整数N (< 10⁵)，请计算不超过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;
}