素数线性筛 之（欧拉筛）

素数线性筛 之（欧拉筛）

最近做项目很久么得刷题了，昨天偶然用到筛素数，突然想整理一下欧拉筛
我们都知道普通筛和埃氏筛会造成一些不必要的浪费，欧拉筛可以改善这一问题。。
还是借蓝桥的一道题说明吧

二、看一下2017年8届蓝桥的2t吧

标题：等差素数列
2,3,5,7,11,13,…是素数序列。

类似：7,37,67,97,127,157 这样完全由素数组成的等差数列，叫等差素数数列。

上边的数列公差为30，长度为6。
2004年，格林与华人陶哲轩合作证明了：存在任意长度的素数等差数列。

这是数论领域一项惊人的成果！
有这一理论为基础，请你借助手中的计算机，满怀信心地搜索：
长度为10的等差素数列，其公差最小值是多少？
注意：需要提交的是一个整数，不要填写任何多余的内容和说明文字。
话不多说，直接在代码中解释吧

#include<bits/stdc++.h>
using namespace std;
const int n=100000;//数组尽量开的大一些，方便查询  
bool jude[n];//判断是否是素数 
int  yes[n]={0};//存放素数 
int len=0;//数组yes的长度 （即有几个素数） 
void init()//重点！！！线性素数筛来了！！！！！ 
{
 
 jude[0]=jude[1]=true;//素数从2开始，0和1显然不是。p掉。 
 for(int i=2;i<=n;i++)
 {
  if(!jude[i])
  yes[len++]=i;//如果是素数，存入 
  for(int j=0;j<=len&&i*yes[j]<=n;j++)
  {
   jude[i*yes[j]]=true;//将数的倍数标记为非素数 
   if(i%yes[j]==0)
   break; //!这里是比较难理解的一点，网上比较普遍的解释是防止重复标记。但究竟是怎么个防止法，就很耐人寻味了~，~ 
//我个人是这样理解的，因为倍数最小是从2而不是1开始的（与质数相乘的话，该质数一定是乘积的最大因子），所以这个是别的数无法代替的(标记上)，但也到头了，
//这可能是一种潜规则吧，整个规律可能就跟这个有关，再深我也悟不透了(*/ω＼*) 
//而后面再遇到的数的乘积随着i的增大，肯定还会再次遇到（因为i是不断递增的可以变成任何数，而只要遇到自己的质因子了，再大的话就可能由别的数组成） 
//所以可以将整个质数循环看作是一个最小质因子循环，如果这个最小质因子能被整除了，
//再往下乘就会重复标记了,因为后面可能在遇到  ####### 下面会配图说明 
  }
 }
 return ;
}

————————————以下是图片说明————————————————————————————————————————————————————————————————

————————————————————————————————————————————————————————————————————————————————————————————————

上面就是本片的灵魂部分了，以下将主函数奉上，可解出本题答案

int main()
{
init();//先将所有素数筛出来 
for(int i=1;i<=n;i++)
{//这个循环的是差值的循环，从小到大慢慢找 
 for(int j=0;j<=len;j++)
 {//这个循环是列举所有的素数个数 
  int f=0;//突破口，到10停止
  for(int k=yes[j];k<=n;k=k+i)
  {//上面列举的是素数个数，这里列举的是素数 
   if(jude[k])//要是这个数被标记过不是素数了，说明这不是等差素数列，重新从下一个循环找 
   break;
   f++;
   if(f>=10)
   { 
    cout<<i<<endl;
    exit(0);
    //成功找到后输出差值是几，并且强制结束 
    //exit(0)是强制退出程序 
   }
   } 
  } 
}
return 0;
}