2,3,5,7,11,13,....是素数序列。
类似:7,37,67,97,127,157 这样完全由素数组成的等差数列,叫等差素数数列。
上边的数列公差为30,长度为6。
2004年,格林与华人陶哲轩合作证明了:存在任意长度的素数等差数列。
这是数论领域一项惊人的成果!
有这一理论为基础,请你借助手中的计算机,满怀信心地搜索:
长度为10的等差素数列,其公差最小值是多少?
注意:需要提交的是一个整数,不要填写任何多余的内容和说明文字。
思路分析:利用埃氏法求素数表,然后再暴力查找,以公差做为最外层循环,第二层循环确定数列首项,第三层循环查找在该等差下是否能够找到连续的n项。
#include <stdio.h>
#include <stdlib.h>
const int N = 10000;//N大小无法确定,可以逐渐增大,直至找到解
int PRIME[N]={0};//素数表
int solve(int n);
void getPrimetable()//埃氏法求素数表
{
int i,j;
for(i=2;i<N;i++)
{
if(!PRIME[i])
{
for(j=i+i;j<N;j=j+i)
PRIME[j]=1;
}
}
}
int main( )
{
int i,j,n;
getPrimetable();
scanf("%d",&n);//等差数列项数
printf("%d",solve(n));
return 0;
}
int solve(int n)
{
int count,i,j,k;
for(i=1;i<N;i++)//公差
{
for(j=2;j<N-n*i;j++)//数列起始位置 ,连续n个,所以j<N-n*i
{
if(!PRIME[j])//出现质数
{
for(k=0;k<n;k++)//连续n个
{
if(PRIME[j+i*k])//出现非质数
break;//结束循环
}
if(k==n)//达到指定个数
return i;//返回公差
}
}
}
return 0;
}