素数笔记

m不必被2~m-1之间的每一个整数去除，只需被2~根号m之间的每个数去除就可以了。例如判别17是否为素数，只需使2~4之间的每一个整数去除。为什么可以做如此简化呢？因为如果m能被2~m-1之间任意整数整除，如果这个数大于根号m，那这个数必定对应的还有一个比根号m小的因子（以16为例，2、8是它的因子，8大于4，2小于4）

1、时间复杂度为O（sqrt(n)）

对于每个数n，其实并不需要从2判断到n-1，我们知道，一个数若可以进行因数分解，那么分解时得到的两个数一定是一个小于等于sqrt(n)，一个大于等于sqrt(n)，据此，上述代码中并不需要遍历到n-1，遍历到sqrt(n)即可，因为若sqrt(n)左侧找不到约数，那么右侧也一定找不到约数。C++代码如下：

[cpp]  view plain  copy

1. bool isPrime_2( int num )  
2. {  
3.      int tmp =sqrt( num);  
4.      for(int i= 2;i <=tmp; i++)  
5.         if(num %i== 0)  
6.           return 0 ;  
7.      return 1 ;  
8. }  

2、更快速的方法

首先看一个关于质数分布的规律：大于等于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个为单元快进，即将方法（2）循环中i++步长加大为6，加快判断速度，代码如下：

[cpp]  view plain  copy

1. bool isPrime_3( int num )  
2. {  
3.                  //两个较小数另外处理  
4.                  if(num ==2|| num==3 )  
5.                                  return 1 ;  
6.                  //不在6的倍数两侧的一定不是质数  
7.                  if(num %6!= 1&&num %6!= 5)  
8.                                  return 0 ;  
9.                  int tmp =sqrt( num);  
10.                  //在6的倍数两侧的也可能不是质数  
11.                  for(int i= 5;i <=tmp; i+=6 )  
12.                                  if(num %i== 0||num %(i+ 2)==0 )  
13.                                                  return 0 ;  
14.                  //排除所有，剩余的是质数  
15.                  return 1 ;  
16. }  

原文：http://blog.csdn.net/huang_miao_xin/article/details/51331710