省赛第二题 等差素数列

标题：等差素数列

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

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

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

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

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

先用素数筛筛出素数，然后暴力

答案：210

#include<iostream> 
#include<math.h>
#include<string.h>
#define X(x) x,x,x,x,x,x,x,x,x,x//如何使数组符相同的值 
using namespace std;
int prim_tab[10010] = {X(X(X(X(1)))),X(1)};//make_tab1
int prim[10000];
int num = 0;
void solve()
{
	int i = 2;
	while(1)//第一层循环根据题意应该是公差 
	{
		int j = 0;
		for(j = 0; j < num/2;j++)//爆破不出来再把控制循环的部分变大 ，这一部分为数列的起始数 
		{
			int k;
			for(k = 0;k < 10;k++)//循环十次都不break就符合 
			{
				if(!prim_tab[prim[j] + k*i])
				{
					break;
				}			
			}
			if(k == 10)//判断 
			{
				cout << i << " " << j <<endl; //i位公差，j为prim数组中对应的起始数的下标 
				system("pause");
			}
		}
		i++;
	}
}


void make_tab1()//这个算法用来判断素数为佳，但是来生成这样的素数表，有点费时 
{
	int i,j;
	int flag = 0;
	for(i = 2;i < 10010;i++)
	{
		flag = 0;
		for(j = 2;j <= sqrt(double(i));j++)
		{
			if(i % j == 0)
			{
				flag = 1;
				break;
			}
		}
		if(flag == 0)
		{
			prim_tab[i] = 1;
			prim[num] = i;
			num++;
		}
	} 
}


void make_tab2()//看了别人的代码，发现另一个算法 
{
	int i,j,k;
	prim_tab[0] = 0;
	prim_tab[1] = 0;
	prim_tab[2] = 1;
	for(i = 2;i < 10010/2;i++)
	{
		if(prim_tab[i] == 1)
		{
			for(j = 2;i * j <10010;j++)
			{
				k = i*j;
				prim_tab[k] = 0;
			}
			prim[num++] = i;
		}
		else
		{
			continue;
		}
	} 
} 

int main()
{
	make_tab2();
	/*
	for(i = 0;i < num;i++)
	{
		cout << su[i] << " ";
	} */
	solve();

	system("pause");
	return 0;
}