2,3,5,7,11,13,…是素数序列。
类似:7,37,67,97,127,157 这样完全由素数组成的等差数列,叫等差素数数列。
上边的数列公差为30,长度为6。
2004年,格林与华人陶哲轩合作证明了:存在任意长度的素数等差数列。
这是数论领域一项惊人的成果!
有这一理论为基础,请你借助手中的计算机,满怀信心地搜索:
长度为10的等差素数列,其公差最小值是多少?
注意:需要提交的是一个整数,不要填写任何多余的内容和说明文字。`
#include <iostream>
#include <set>
/* run this program using the console pauser or add your own getch, system("pause") or input loop */
using namespace std;
typedef long long LL;
const int N=5000;//素数最高到5000
const int M=1000;//公差最高到1000
set<int> sett;
bool isprime(LL n){//判断是否是素数
for(int i=2;i<n/2;i++){
if(n%i==0)return false;
}
return true;
}
int f(LL a[],int y){
for(int item=0;item<y;item++){//遍历每个素数
for(int i=1;i<M;i++){ //公差
int first=a[item];
int p;
for(p=0;p<10;p++){//循环加10次公差,如果10次相加的结果都是素数,返回该公差
first+=i;
if(sett.find(first)==sett.end())break;
}
if(p==9)return i;
}
}
}
int main(int argc, char *argv[]) {
LL a[N];
int t=5;
int y=0;//素数的个数有多少个
while(t<N){
if(isprime(t)){//将素数放进数组跟set中
a[y]=t;
sett.insert(t);
y++;
}
t++;
}
cout<<f(a,y);
return 0;
}
答案为210