素数筛

我们研究一下质数分布规律：大于等于5的质数一定和6的倍数相邻。例如5和7，11和13,17和19等；

神仙筛法：

我们研究一下质数分布规律：大于等于5的质数一定和6的倍数相邻。例如5和7，11和13,17和19等；

证明：令x≥1，将大于等于5的自然数表示如下：
······ 6x-1，6x，6x+1，6x+2，6x+3，6x+4，6x+5，6(x+1），6(x+1)+1 ······
可以看到，不在6的倍数两侧，即6x两侧的数为6x+2，6x+3，6x+4，由于2(3x+1)，3(2x+1)，2(3x+2)，所以它们一定不是素数，再除去6x本身，显然，素数要出现只可能出现在6x的相邻两侧。 这里要注意的一点是，在6的倍数相邻两侧并不是一定就是质数。

此时判断质数可以6个为单元快进，i++步长加大为6，加快判断速度，原因是，假如要判定的数为n，则n必定是6x-1或6x+1的形式，对于循环中6i-1，6i，6i+1,6i+2，6i+3，6i+4，其中如果n能被6i，6i+2，6i+4整除，则n至少得是一个偶数，但是6x-1或6x+1的形式明显是一个奇数，故不成立；另外，如果n能被6i+3整除，则n至少能被3整除，但是6x能被3整除，故6x-1或6x+1（即n）不可能被3整除，故不成立。综上，循环中只需要考虑6i-1和6i+1的情况，即循环的步长可以定为6，每次判断循环变量k和k+2的情况即可

判断单个素数最快：

bool isprime(int n)
{
	if(n==0 || n==1)    //先判断这些小于5的数 
		return false;
	if(n==2||n==3)    
		return true;
	if(n%6!=1 && n%6!=5)  //判断数是否是6x-1或者6x+1 
		return false;
	for(int i=5;i<=sqrt(n);i+=6)
	{
		if(n%i==0 || n%(i+2)==0) //每次只用判断是否是它两的倍数 
			return false;	
	}
	return true;		 //都满足了，即为质数 
}

欧拉筛 （O(n)，求一个区间的所有素数最快）

int seive(int n){
    int cot = 0;
    memset(a, 0, sizeof(a));
    
    for(int i = 2; i <= n; i++){   //遍历筛去所有最大因数是i的合数
        if(a[i] == 1) b[++cot] = i;
        for(int j = 1; j <= cot  && i * b[j] <= n; j++){
        	
            a[i * b[j]] = 1;
            
            if(!( i % b[j])) break;  //找到i的最小素因数，减少重复计算 
           
        }
    }
    
    return cot;
}