判断素数的高效的方法

      以前判断素数的方法是：将每个值传入一个函数，这个函数里面用 for循环判读是否为素数。这种方法的效率比较低，所以经过在网上查找资料和大佬的整理，出现里另外一种效率较高的判断素数的方法——孪生素数法。

      首先我们要知道什么是孪生素数：素数对（p, p + 2）称为孪生素数。

      孪生素数中有一个定理：大于3的素数只分布在6n-1和6n+1两数列中。（n非0自然数）。

        这个定理说明了，大于3的素数一定是6的倍数相邻，即6n-1和6n+1,例如11和13，17和19等等。

综上说明得： 大于3的素数 ==》分布在6n-1和6n+1这两个数列中。

这是一个逻辑关系，即“分布在6n-1和6n+1”是“大于3的素数”的必要条件。由逻辑关系得到，它的逆否定理得到 “这个是不分布在6n-1和6n+1数列中”==》“这个数不是大于3的素数。” 即：不再这两个数列中的数都不是素数（2、3单独讨论，这里除外）。再有逻辑关系得：如果一个数在6n-1和6n+1数列中，这个数不一定为素数，需要另外判断。

 

    综上，我们可以得到一个关于素数的规律：大于3的素数一定和6的倍数相邻。

 

   由此，我们可以得出一个判断素数的算法：

  1、  当被判断数为1时，直接判断为不是素数

 if(num<2)
        return NO;

2、当被判断素数为2或者3时，直接判断为素数       

   

if(num==2||num==3)
        return OK;

  3、当被判断数不再6n-1和6n+1两个数列任何一个数列是，直接判断为不是素数。

    if(num%6!=1&&num%6!=5)//被判断数即不再6n-1数列也不再6n+1数列才进入
        return NO;

4、剩余的数都在6n-1和6n+1数列，判断是否为素数

      经过上面3重筛选，剩余的数全部都是在6n-1和6n+1两数列之中的数，有上面的逻辑关系得到这个数不一定为素数，所以我们还需要两外判断。 

    int temp=sqrt(num);
    for(int i=5;i<=temp;i+=6)
        if(num%i==0||num%(i+2)==0)
            return NO;
    return OK;

判断素数就是取模，直到有一个数取模为0说明这个数不是素数，停止循环。或者直到循环结束，跳出循环也没有取模为0的，说明这个数为素数。由此可得，如果被判断数对某些数取模一定不为0那么就可以直接跳过这些数，去判断可能取模为0的数。

     我们知道只需要判断2到根号下被判断数即可得出结论，但是目前的被判断数都在6n-1和6n+1这两个数列当中，此时判断素数可以以6为单位快进，加快判断素数。原因如下：被判断数在6n-1和6n+1数列中，6n是偶数，那么6n-1和6n+1一定是奇数。对于循环中6i,6i-1,6i+1,6i+2,6i+3,6i+4,其中被判断数%(6i)，%(6i+2)，%(6i+4)一定都不为0，因为6i,6i+2,6i+4一定为偶数，而被判断数一定为奇数，所以不为0，不用判断。如果%(6i+3)!=0,原因如下：假设被判断数能被%(6i+3）==0，说明被判断数能被6i+3整数，那么也就是说被判断数能被3整出，但是被判断数属于6n-1与6n+1这两个数列之一，不可能被3整除，所以假设不成立，被判断数%(6i+3)!=0,无判断。   

举个例子：2为(6i+2)数列中的，为偶数，不可能被6n-1和6n+1数列中的数取模为0，直接跳过。3为(6i+3)数列中的数，因为6n-1和6n+1数列不可能整除3，所以直接跳过。4属于6i+4数列，为偶数，不可能被6n-1和6n+1数列中的数取模为0，直接跳过。5为6i-1数列，需要判断。6属于6i数列，为偶数，不可能被6n-1和6n+1数列中的数取模为0，直接跳过。7属于6i+1数列，需要判断。

所以循环可以直接从5开始，以6为单位加速，只需要判断(num%i)是否为0和(num%(i+2))是否为0即可，这两个又一个为0就可以判断为不是素数。如果直到循环结束也没有为0的没，说明这个数是素数。（这里的i一直属于6i-1数列，i+2一直属于6i+1数列）。

bool judge(int num){
    if(num<2)       //被判读数为1
        return NO;
    if(num==2||num==3)  //被判断数为2或3
        return OK;
    if(num%6!=1&&num%6!=5)   //被判读数不再6n-1和6n+1数列
        return NO;
    int temp=sqrt(num);
    for(int i=5;i<=temp;i+=6)     //被判断数在6n-1和6n+1数列
        if(num%i==0||num%(i+2)==0)
            return NO;
    return OK;
}