判断质数

1.暴力扫描

    直接从2扫到该数-1，思路好理解，但效率极低，时间复杂度为O(n)。
    代码如下：
     

bool isPrime_1(int num)
    {
        int tmp=num-1;
        for(int i=2;i<=tmp;i++)
          if(num%i== 0)
             return 0 ;
        return 1 ;
    }

2.优化暴力

    可以先去判断该数是否为2,3,5,7的倍数，以减少循环。
    时间复杂度为O(sqrt(n))，效率提高不少。
    代码如下：
   

bool isPrime_2(int num)
    {
         if(num%2==0||num%3==0||num%5==0||num%7==0)return 0;
         int tmp=sqrt(num);
         for(int i=7;i<=tmp;i+=2)    //从这数开始扫描
            if(num%i==0)
              return 0;
         return 1;
    }

 3.打表判断

    可用上述方法打出所需数据范围内质数表，再用bool数组判断。
    可使时间复杂度变为O(1)，但需考虑其内存大小。

 

4.质数规律判断法

 

    证明：令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 prime(int n)
{
    if(n==1)                return false;
    else if(n==2||n==3)        return true;
    else if(n%6!=1&&n%6!=5)    return false;
    int tmp=(int)sqrt(n);
    for(int i=5;i<=tmp+1;i+=6)
        if(n%i==0||n%(i+2)==0)
            return false;
    return true;
}

    该方法将大大提升判断速度，是仅次于打表速度的第二算法，而且无需有内存的顾虑。

5.Eratosthenes算法

    用一个长为N+1的数组保存信息(0为质数，1为非质数)，初始化为0，每局每一个质数并把x的倍数标为非质数。
    时间复杂度为O(nloglogn),效率接近线性。
    代码如下：
 

void primes(int n)
{
    memset(v,0,sizeof(v));
    for(int i=2;i<=n;i++)
       if(v[i])    cotinue;
       cout<<i<<endl；
       for(int j=i;j<=n/i;j++)
           v[i*j]=1;
       }
}