1007 素数对猜想(20 分)
让我们定义dn为:dn=pn+1−pn,其中pi是第i个素数。显然有d1=1,且对于n>1有dn是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数N
(<105),请计算不超过N
的满足猜想的素数对的个数。
输入格式:
输入在一行给出正整数N
。
输出格式:
在一行中输出不超过N
的满足猜想的素数对的个数。
输入样例:
20
输出样例:
4
关键之处:
- 要会写素数判定函数。
代码:
#include <stdio.h>
#include <math.h>
#define MAXN 100000
int isPrime ( int n );
int main(void) {
int N, i, k, count;
int a[MAXN] = {0};
scanf("%d", &N);
/*将该整数前的所有素数存入数组a[]*/
k = 0;
for ( i = 2; i <= N; i++ )
if ( isPrime(i) ) a[k++] = i;
/*遍历这个全是素数的数组,相邻两项为2,则计数*/
count = 0;
for ( i = 0; i < k; i++ )
if ( a[i+1] - a[i] == 2 ) count++;
printf("%d\n", count);
return 0;
}
/*判断素数的函数*/
int isPrime ( int n ) {
int flag, i;
/*小于等于1的数一定不是素数,标记为0*/
if ( n <= 1 ) flag = 0;
else {
/*判断素数的最小范围是到根号n为止,数学上可证明*/
for ( i = 2; i <= sqrt(n); i++ )
if ( n % i == 0 ) break; //如果在根号n之前被整除了,跳出循环。
if ( i > sqrt(n) ) flag = 1; //循环正常结束,说明是素数。
else flag = 0; //循环非正常结束,说明不是素数。
}
return flag;
}