题目
让我们定义dn为:dn=pn+1−pn,其中pi是第i个素数。显然有d1=1,且对于n>1有dn是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数N(<10^5),请计算不超过N的满足猜想的素数对的个数。
输入格式:
输入在一行给出正整数N。
输出格式:
在一行中输出不超过N的满足猜想的素数对的个数。
输入样例:
20
输出样例:
4
解题思路
首先理解题意,分析后得出,这道题目需要两个函数,素数判断+dn求解;
1、素数求解:
涉及的两个概念:
素数:只有1和自身是它的除数,0、1均不为素数
除数:24/8–>24被除数,8除数;
函数求解思路:输入数字n,循环判断它的除数是不是只有自身和1
其中除数从1到n循环当然可以,但是有一种更加优化的方法,也是这里我选择采用的方法,根据数学分析(这里不再赘述具体的分析过程)可得除数最大只可能到√n,所以只需循环到√n即可。
2、dn求解
通过dn计算相邻的两个素数之差是否为2,如果是则素数对加一,反正不是。
代码如下:
#include <stdio.h>
#include <math.h>
int prime (int n){//判断n是否为素数
int i,k;k=sqrt(n);
if(n==1||n==0)
return 0; //1 0不是素数,需单独考虑
for(i=2;i<=k;i++)
if(n%i==0) return 0;
return 1;
}
int dn(int m,int n,int* x){//求dn
if(m>=n){
*x= m-n; //注意返回值要使用指针
}
*x=n-m;
}
int main(){
int N,i,j,k,x,m;//输入正整数N
int a[10000];
int b[10000];
scanf("%d",&N);
m=N;
for(j=N-1;j>=0;j--){//将整数N拆分开放在数组a中
a[j]=N;
N--;
}
for(i=0,k=0;i<m;i++){//判断a[i]是否为素数
if(prime(a[i])){
b[k]=a[i];//b中存放所有小于N的素数
k++;
}
}
j=0;i=0;
while(j<k-1){//求素数的差
dn(b[j],b[j+1],&x);//返回函数的值
if(x==2){
i++;
}
j++;
}
printf("\n%d",i);
return 0;
}
总结
这次的代码不是很复杂,一定要捋清楚代码参数之间的关系,而且注意=和==之间的区别,好几次都是因为二者没有分清楚,导致浪费了很多时间。