题目一:
公差最小的等差素数列
Problem Description
2,3,5,7,11,13,....是素数序列。
类似:7,37,67,97,127,157 这样完全由素数组成的等差数列,叫等差素数数列。
上边的数列公差为30,长度为6。
2004年,格林与华人陶哲轩合作证明了:存在任意长度的素数等差数列。
这是数论领域一项惊人的成果!
有这一理论为基础,请你借助手中的计算机,满怀信心地搜索:
长度为10的等差素数列,其公差最小值是多少?
注意:需要提交的是一个整数,不要填写任何多余的内容和说明文字。
Input
无
Output
一个整数,即程度为10的等差素数列的最小公差
思路:1.搞定两个数组,prime[]数组放从小到大的素数,num[]数组标记哪些数为素数。
2.先枚举公差,再枚举起始素数,然后判断起始素数+k倍公差(k从1~9)所得到的数字是否也是素数(num[]判断),然后再判断起始素数+10倍公差是否是素数(因为要求长度为10,长度为11的等差素数列也满足上述条件,但不符合题意)。
3.标记flag,遇到解即跳出循环。
#include<iostream>
using namespace std;
#define N 1000009
int prime[N],tot=0;
int num[N];
int check(int x)
{
if(x==2) return 1;
if(x%2==0) return 0;
else{
for(int i=2;i*i<=x;i++)
if(x%i==0)
return 0;
return 1;
}
}
void init()
{
for(int i=2;i<N;i++)
if(check(i))
{
prime[tot++]=i; //把素数筛出来从小到大放到prime数组中
num[i]=1; //整数数组里标记。
}
}
int main()
{
init();
int flag;
for(int i=1;i*10<N;i++){ //枚举公差
for(int j=0;j<N;j++){ //枚举起始素数
flag=0;
int tmp=prime[j];
for(int k=1;k<=9;k++){ //枚举长度
int ne_tmp=tmp+k*i;
if(num[ne_tmp]!=1)
break;
if(k==9&&num[tmp+10*i]!=1)//确保长度为10,而不是11,12等
{
flag=1;
cout<<i;
}
}
if(flag) break;
}
if(flag) break;
}
return 0;
}