等差素数列(蓝桥杯)

139 阅读1分钟


标题:等差素数列

2,3,5,7,11,13,....是素数序列。
类似:7,37,67,97,127,157 这样完全由素数组成的等差数列,叫等差素数数列 。
上边的数列公差为30,长度为6。

2004年,格林与华人陶哲轩合作证明了:存在任意长度的素数等差数列。
这是数论领域一项惊人的成果!

有这一理论为基础,请你借助手中的计算机,满怀信心地搜索:

长度为10的等差素数列,其公差最小值是多少?

注意:需要提交的是一个整数,不要填写任何多余的内容和说明文字

程序代码:

#include<stdio.h>
#include<math.h>
const int N=1000010;
int isprime[N],prime[N];

int main()
{
	int i,j,t=0,d,sum;
	for(i=2;i<=sqrt(N);i++)
		if(isprime[i]==0)
			for(j=2*i;j<N;j+=i)
				isprime[j]=1;
	for(i=2;i<N;i++)
		if(isprime[i]==0)
			prime[t++]=i;
			
	d=1;
	while(1)
	{
		for(i=0;i<t;i++)
		{
			sum=1;
			for(j=i+1;j<t;j++)
			{
				if(prime[j]-prime[i]==sum*d)
					sum++;
				else if(prime[j]-prime[i]>sum*d)
					break;
			}
			if(sum>=10)
				break;
		}	
		if(sum>=10)
			break;
		d++;
	}
	printf("%d\n",d);
	return 0;
}