参考自:蓝桥杯2017第二题。
标题:等差素数列
2,3,5,7,11,13,....是素数序列。
类似:7,37,67,97,127,157 这样完全由素数组成的等差数列,叫等差素数数列。
上边的数列公差为30,长度为6。
2004年,格林与华人陶哲轩合作证明了:存在任意长度的素数等差数列。
这是数论领域一项惊人的成果!
有这一理论为基础,请你借助手中的计算机,满怀信心地搜索:
长度为10的等差素数列,其公差最小值是多少?
注意:需要提交的是一个整数,不要填写任何多余的内容和说明文字。
用了很长时间自己写了一个算法,发现最多可以计算长度为10的等差素数数列。也就是说代码不够健全。正在寻找解决方案。
#include<stdio.h>
#include<stdbool.h>
#include<math.h>
#define N 15000
//算法:长度为10的等差素数列,其公差最小值
int count;
bool ismin;
bool isPrime(int n){//判断是否为素数,这个算法性能不佳,可以使用线性素数筛
int i;
for(i=2;i<=sqrt(n);i++){
if(n%i==0) return false;
}
return true;
} //一个数的因数一定一个大于等于sqrt(n),一个小于等于sqrt(n)。
void test(int n,int temp){
if(count==10){
printf("temp=%ld n=%ld\n",temp,n-10*temp);
ismin = true;
return;
}
else if(isPrime(n+temp)&&count!=10){
++count;
test(n+temp,temp);//递归
}
else {
count = 0;
return;
}
}
int main(){
int n,temp;
for(temp=1;temp<=N;temp++){
for(n=2;n+temp<N;n++){ //temp指的是公差,n是首项。基本思路是运用
//等差数列的公式an=a1+(n-1)d,即知道了首
//项和公差就可以确定一个等差数列。
if(isPrime(n)){
test(n,temp);
if(ismin) return 0;
}
}
}
}
顺便补充一下数学上的一些基本概念:
素数:大于1的自然数中,除了1和它本身以外不再有其他因数。
2是唯一的偶素数,其他素数都是奇素数。
比1大但不是素数的书称为合数。任何一个合数都可以表示成一个质数和一个数的乘积。