标题:等差素数列
2,3,5,7,11,13,…是素数序列。
类似:7,37,67,97,127,157 这样完全由素数组成的等差数列,叫等差素数数列。
上边的数列公差为30,长度为6。
2004年,格林与华人陶哲轩合作证明了:存在任意长度的素数等差数列。
这是数论领域一项惊人的成果!
有这一理论为基础,请你借助手中的计算机,满怀信心地搜索:
长度为10的等差素数列,其公差最小值是多少?
注意:需要提交的是一个整数,不要填写任何多余的内容和说明文字。
——————————————————————————————————————————————
- 素数打表,列一个纯素数数组,根据素数规律从30开始取公q为6的倍数,即q+=6;
答案:210
#include<iostream>
#include<cstdio>
using namespace std;
const int N=1e6+5;
int a[N];
int ct=0;
void prime()
{
bool p[N];
fill(p,p+N,1);
p[1]=0;
for(int i=2;i<N;i++)
{
if(p[i])
for(int j=i+i;j<N;j+=i)
p[j]=0;
}
for(int i=2;i<N;i++)
{
if(p[i])a[ct++]=i;
}
}
int main()
{
prime();
int q=30;
int ctt;
while(1)
{
q+=6;//
for(int s=3;s<100;s++)//等差数列第一个素数,这里从7开始
{
int h=s;
ctt=1;
for(int j=h+1;j<ct;j++)//后面搜索符合q的素数加入以a[s]为头等差数列
{
if(a[j]-a[h]>q)break;//对于公差q,以a[s]开头最长处
if(a[j]-a[h]==q)
{
cout<<ctt<<' '<<a[j]<<endl;//
h=j;
ctt++;
if(ctt==10)
{
printf("%d\n",q);
return 0;
}
}
}
}
}
return 0;
}