让我们定义dn为:dn=pn+1−pn,其中pi是第i个素数。显然有d1=1,且对于n>1有dn是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数N
(<105),请计算不超过N
的满足猜想的素数对的个数。
输入格式:
输入在一行给出正整数N
。
输出格式:
在一行中输出不超过N
的满足猜想的素数对的个数。
输入样例:
20
输出样例:
4
这题我采用了打印素数表的方式,将<N的素数全部装入数组里面,然后进行遍历。
因为<N的数里,素数是少的,减少了计算次数,用素数来对<N的数字做除法,能整除便不是素数。(该方法参见于mooc翁恺老师的C语言课)
(注意这里适用于N不大的情况,N要是过大,这种方式会运行超时,举个例子,一般采用sqrt(N)的数学函数的方式求素数,10的平方根是3,而10里面有4个素数,所以计算次数会多余采用数学函数的方式)
//素数对猜想
#include<stdio.h>int isprime(int x, int prime[], int length); //用已知素数表的素数去做除法
int main()
{
const int number=100000;
int prime[number]={2};
int i=3,cnt=1,N=0,sum=0;
if(scanf("%d", &N)){
for(i=3;i<=N;i++){
if(isprime(i, prime, cnt)){
prime[cnt++]=i;
}
}
// for(i=0;i<cnt;i++){
// printf("%d\n", prime[i]); 检测素数表是否打印成功
// }
for(i=0;i<cnt;i++){
if(prime[i+1]-prime[i]==2) sum++;
}
printf("%d", sum);
}return 0;
}
int isprime(int x, int prime[], int numberofknownprime)
{
int ret=1;
for(int i=0;i<numberofknownprime;i++)
{
if(x%prime[i]==0){
ret=0;
break;
}
}
return ret;
}