2,3,5,7,11,13,…是素数序列。
类似:7,37,67,97,127,157 这样完全由素数组成的等差数列,叫等差素数数列。
上边的数列公差为30,长度为6。
2004年,格林与华人陶哲轩合作证明了:存在任意长度的素数等差数列。
这是数论领域一项惊人的成果!
有这一理论为基础,请你借助手中的计算机,满怀信心地搜索:
长度为10的等差素数列,其公差最小值是多少?
注意:需要提交的是一个整数,不要填写任何多余的内容和说明文字。
思路:把所有素数放在一个数组里,方便查找枚举。
用枚举法寻找,先枚举首个元素,再枚举公差,再枚举其它九个元素。直到其它的元素都符合条件
用set.find()方法找起来比较快,省的再次判读元素。
#include<bits/stdc++.h>
using namespace std;
set<int>all;
bool isprime(long long a){
for(int i = 2; i < a/2; i++){
if(a%i==0) return false;
}
return true;
}
//枚举法寻找
int f(long long a[],int n){
int ans=2,j=0;
for(int i = 0; i < n; i++){//枚举首项
long long first = a[i];
for(int j = 1; j < a[n-1]-first; j++){//枚举公差
int m=first;
for(int z = 1; z < 10; z++){//枚举个数
m += j;
if(all.find(m)==all.end())break;
if(m > a[n - 1])break;
if(z == 9)return j;
}
}
}
}
int main(){
int i=1,max=10,t=5;
long long a[505];
a[0] = 2;
all.insert(2);
all.insert(3);
while(i < 500){
if(isprime(t)){
a[i++] = t;
all.insert(t);
}
t++;
}
cout<<f(a,500);
return 0;
}