1007. 素数对猜想 (20)

如题：http://www.patest.cn/contests/pat-b-practise/1007

让我们定义 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的满足猜想的素数对的个数。

输入样例：

20

输出样例：

4

 

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

int prime[100001];
int vis[100001];

int isp(int x)
{
 int i;
 for(i=2;i*i<=x;i++)
  if(x%i==0)
   return 0;
 return 1;
}
int main()
{
 int n;
 cin>>n;
 int i,j;
 memset(prime,1,sizeof(prime));
 prime[0]=prime[1]=0;
 vis[0]=vis[1]=0;
 for(i=2;i<=n;i++)
 { 
  if(!vis[i])
  {
   if(isp(i))
   {
    for(j=2*i;j<=n;j+=i)
    {
     prime[j]=0;
     vis[j]=1;
    }
   }
  }
 }
 int cnt=0;
 for(i=2;i<=n-2;i++)
 {
  if(prime[i]&&prime[i+2])
   cnt++;
 }
 cout<<cnt<<endl;
}