2,3,5,7,11,13,…是素数序列。
类似:7,37,67,97,127,157 这样全由素数组成的等差数列,叫等差素数数列。
上边的数列公差为30,长度为6。
2004年,格林与华人陶哲轩合作证明了:存在任意长度的素数等差数列。
这是数论领域一项惊人的成果!
有这一理论为基础,请你借助手中的计算机,满怀信心地搜索:
长度为10的等差素数列,其公差最小值是多少?
思路:因为这题是填空题,求长度为10的等差素数列,猜想最大的数应该不会超过10000,因此先枚举2~10000每个数,判断一下是否是素数;先枚举首项a,再枚举公差d,再枚举10个数k,判断一下a + d * k是否是素数,如果是,cnt++,如果cnt ==10,就输出公差d
#include <iostream>
#include <algorithm>
using namespace std;
bool vis[10000];
bool check (int x)
{
for (int i = 2 ; i < x ; i ++)
{
if (x % i == 0) return false;
}
return true;
}
int main()
{
for (int i = 2 ; i <= 10000 ; i ++)
{
if (check(i)) vis[i] = true;
}
for (int a = 2 ; a <= 10000 ; a ++)
{
for (int d = 1 ; d < 1000 ; d ++)
{
int cnt = 0;
for (int k = 0 ; k < 10 ; k ++)
{
if (vis[a + k * d]) cnt ++;
}
if (cnt ==10) cout << d << endl;
}
}
return 0;
}