2,3,5,7,11,13,.... 是素数序列。 类似:7,37,67,97,127,1577,37,67,97,127,157 这样完全由素数组成的等差数列,叫等差素数数列。
上边的数列公差为 3030,长度为 66。
20042004 年,格林与华人陶哲轩合作证明了:存在任意长度的素数等差数列。 这是数论领域一项惊人的成果!
有这一理论为基础,请你借助手中的计算机,满怀信心地搜索:
长度为 1010 的等差素数列,其公差最小值是多少?
答案:最小公差为210
运行代码如下:
注意:在算法竞赛中,一般会指定编译器,不要使用system("pause")语句,会无法输出从而被判断错误。
#include<iostream>
#include<cmath>
using namespace std;
int isprime(int x) //判断是否为素数
{
for (int i = 2; i <= sqrt(x); i++)
{
if (x%i == 0)
return 0;
}
return 1;
}
int main()
{
int cnt = 0; //计数,cnt<10
for (int i = 2; i < 50000; i++)
{
if (isprime(i) == 1)
{
for (int j = 1; j < 1000; j++) //枚举公差
{
for (int t = 0; t < 10; t++) //找到十个数
{
if (isprime(i + j * t) == 1) //通项均为素数
{
cnt++;
}
else
{
cnt = 0;
break;
}
if (cnt == 9) //计数达到9说明已有10个数
{
cout << j << endl;
goto mark;
}
}
}
}
}
mark:
return 0;
}