2017年蓝桥杯C/C++程序设计B组-等差素数列

欢迎大家来到zz菜鸟的博客，近期我在备战蓝桥杯，会发布一些我对部分题目的看法。希望大家能够多多支持。

题目描述

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

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

有这一理论为基础，请你借助手中的计算机，满怀信心地搜索：
长度为10的等差素数列，其公差最小值是多少？

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

思路分析

• 进行素数的判断
• 对不同首项进行枚举
• 在相同首项下对不同公差进行枚举
• 直到满足十个数

#include<iostream>
using namespace std;
bool IsPrime(int a)   //判断是否是素数
{
	for (int i = 2; i <= a/2; i++)
	{
		if (a % i == 0) return 0;
	}
    return 1;
}
int main()
{ 
     //prime用来存储素数   isprime中为1表示该数为素数，为0则反之
	long prime[10000] = {2,3}, isprime[50000] = { 0,1,1,0 };  
	long t = 5;
	int index = 2;
	while (t<=50000)      //找出前50000个数中的素数
	{
		if (IsPrime(t))
		{
			prime[index++] = t;
			isprime[t - 1] = 1;   //注意下表减一
		}
		else
		{
			isprime[t - 1] = 0;
		}
		t++;
	}
	for (int i = 0; i < 5000; i++)    //首项循环
	{
		long m = prime[i];
		 //公差循环，注意此处公差不能设太大，若设太大isprime[m + j * k-1]会溢出
		for (int j = 1; j < 1000; j++)  
		{
			 int k=1;
			for (k; k < 10; k++)   //找10个数
			{
				if (isprime[m + j * k-1] == 0) break;   
			}
			if (k == 10) {
				cout << j;
				return 0;
			}
		}
	}
	return 0;
}