C语言B组
标题:等差素数列
2,3,5,7,11,13,….是素数序列。
类似:7,37,67,97,127,157 这样完全由素数组成的等差数列,叫等差素数数列。
上边的数列公差为30,长度为6。
2004年,格林与华人陶哲轩合作证明了:存在任意长度的素数等差数列。
这是数论领域一项惊人的成果!
有这一理论为基础,请你借助手中的计算机,满怀信心地搜索:
长度为10的等差素数列,其公差最小值是多少?
注意:需要提交的是一个整数,不要填写任何多余的内容和说明文字。
用线性素数筛把10^6内的素数筛出来,然后从小到大枚举公差然后去验证。
答案:210
#include <algorithm>
#include <string.h>
#include <iostream>
#include <stdio.h>
#include <string>
#include <vector>
#include <queue>
#include <map>
#include <set>
using namespace std;
const long long N = 1000010;
int dp[N]={1,1,0};
int prim[N],tot = 0;
void init()
{
for(long long i = 2 ; i < N ; i ++)
{
if(dp[i])continue;
prim[tot++]=i;
for(long long j = i ; j * i < N ; j ++){
dp[i*j] = 1;
}
}
}
int main()
{
init();
printf("%d\n",tot);
for(int i = 1 ; i*10 < N ; i ++){
for(int j = 0 ; j < tot ; j ++){
int flag = 1,temp = prim[j];
for(int k = 1 ; k < 6 ; k ++)
{
if(temp + i >= N || dp[temp + i] == 1){
flag = 0;break;
}else{
temp = temp + i;
}
}
if(flag == 1){
printf("%d %d\n",i,prim[j]);
return 0;
}
}
}
return 0;
}
//210 199