标题:等差素数列
2,3,5,7,11,13,....是素数序列。
类似:7,37,67,97,127,157 这样完全由素数组成的等差数列,叫等差素数数列。
上边的数列公差为30,长度为6。
2004年,格林与华人陶哲轩合作证明了:存在任意长度的素数等差数列。
这是数论领域一项惊人的成果!
有这一理论为基础,请你借助手中的计算机,满怀信心地搜索:
长度为10的等差素数列,其公差最小值是多少?
注意:需要提交的是一个整数,不要填写任何多余的内容和说明文字。
解题思路:可以先把一定范围之内的素数存起来,枚举这些作为首项,枚举公差,看看是不是10个都是素数,一开始提前存储素数可以用欧拉筛法
#include<iostream>
#include<cstdio>
using namespace std;
const int N = 100005;
bool is_prime[N];
int prime[N/2], cnt=0;
bool check(int a, int d){
for(int i=1; i<=9; ++i){
int num = a+i*d;
if(is_prime[num])return false;
}
return true;
}
int main()
{
is_prime[1]=true;
for(int i=2; i<N; ++i){
if(!is_prime[i]){
prime[cnt++] = i;
}
for(int j=0; i*prime[j]<N && j<cnt; ++j){
is_prime[i*prime[j]]=true;
if(i%prime[j]==0)break;
}
}
for(int i=1; i<10000; ++i){
for(int j=0; j<cnt; ++j){
if(prime[j]+9*i>=N)break;
if(check(prime[j], i)){
printf("%d", i);
return 0;
}
}
}
return 0;
}