浅谈C++如何求等差素数列

题目

标题：等差素数列

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

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

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

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

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

题解

絮絮叨叨（骂骂咧咧

一开始看到这道题还是有点懵的，毕竟我个数学小白，对素数什么的最发怵了。

然后找了好多大佬的题解都没看明白，甚至有一个大佬的代码看的我晕头转向~

然后终于被我找到一份能看懂并且觉得非常正确的代码，思路如下：

思路

两层循环，一层循环用于循环公差，一层循环用于循环起始素数。

需要注意的是，内层循环起始素数的时候，不能无边界循环下去，要设置一个上限，否则外层循环永远无法走到下一个公差（自己写的时候自以为是犯的错

内层循环走的时候，只需要判断：
①这个数是不是素数（作为起始素数最基本的条件）
②判断从这个素数开始，以cha为公差能否存在连续10个等差的素数。【用ok函数来判断的】

如果以上两个条件都满足，则这就是我们要找的长度为10的等差素数列，其公差的最小值

因为我们是从小到大找的，那我们找到的满足条件的第一个就是答案~

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll maxn=1e6+50;
ll a[maxn];
bool ok(ll n,ll cha)
{
	for(ll i=0;i<10;i++)
	{
		if(!a[n+i*cha])return 0;
	}
	return 1;
}

int main()
{
	a[1]=0;
	a[2]=1;
	a[3]=1;
	for(ll i=4;i<=1000000;i++)
	{
		bool flag=0;
		for(ll j=2;j*j<=i;j++)
		{
			if(i%j==0)
			{
				flag=1;
				break;
			}
		}
		if(flag)a[i]=0;
		else a[i]=1;
	}
	
	for(ll cha=1;;cha++)
	{
		for(ll i=2;i<1000000;i++)
		{
			if(a[i]&&ok(i,cha))
			{
				printf("%lld\n",cha);
				return 0;
			}
		}
	}
}

后记

其实我对素数一直都怀有敬畏之心，希望能找个时间把素数的相关算法摸摸透，把板子整理齐全~（先给自己挖个坑

要是整理好了，我就把链接更新上来！（咕咕咕~