17.leetcode题目204： Count Primes(不是最快的！)

前言：这道题目一开始觉得挺简单的，编写一个判断每一个数是否是质数的函数isPrime：

</pre><pre name="code" class="cpp"> bool isPrime(int m){
    if(m<=1){ 
    return false; 
  } 
//for(int i=2;i<m/2;i++){ //根据 hint 2修改 
for(int i=2;i<=sqrt(m);i++){ //根据 hint 3修改 
      if(m%i==0){ 
       return false; 
    } 
  } 
   return true; 
}

</pre><p>然后对于小于n的每一个数调用isPrime函数，算法时间复杂度O(n^1.5)。</p><p>这肯定是能算出来的，但是运行超时。</p><p>所以，一步步看hint。</p><p></p><pre name="code" class="cpp">class Solution {
   public: 
     int countPrimes(int n) { 
          int count=0; 
          bool *isPrime=new bool[n]; 
          memset(isPrime, true, sizeof(bool) * n);//初始化语句，memset()
       for(int i=2;i*i<n;i++){ //i<=sqrt(n);只用判断前sqrt(n)个数，因为大于sqrt(n)的数如果是某个数的因数，则该数一定已经被标记 
         if(isPrime[i]==true){ 
      for(int j=i*i;j<n;j+=i){ //标记i的倍数，从第i*i开始标记，因为比如i=3，i的倍数为6、9、12 ...6为2的倍数，已经被2的倍数标记为false 
       isPrime[j]=false; 
    } 
   } 
 } 
     for(int i=2;i<n;i++){ 
       if(isPrime[i]==true){ 
        count++; 
      } 
} 
    delete[] isPrime; 
   return count; 
    }
 };

进一步改进：除了2以外，只有奇数才有可能是质数，偶数不可能为质数，所以计数时只考虑奇数。

class Solution {
    public: 
    int countPrimes(int n) {
    if(n<=2){ 
       return 0; 
    } 
    int count=1;
    bool *isPrime=new bool[n]; 
    memset(isPrime, true, sizeof(bool) * n);//初始化语句，memset() 
   for(int i=3;i*i<n;i+=2){ //i<=sqrt(n);只用判断前sqrt(n)个数，因为大于sqrt(n)的数如果是某个数的因数，则该数一定已经被标记 
  if(isPrime[i]==true){  
    for(int j=i*i;j<n;j+=2*i){ //标记i的倍数，从第i*i开始标记，因为比如i=5，i的倍数为10、15、20、25、30、35 ...15为3的倍数，已经被3的倍数标记为false,20为偶数,不计数.30为偶数,我们从35开始表示,且只标记奇数,即j+2*i标记为false 
    isPrime[j]=false; 
   } 
 } 
}
    for(int i=3;i<n;i+=2){ //只对奇数进行判断
       if(isPrime[i]==true){ 
         count++; 
   } 
 } 
    delete[] isPrime; 
   return count; 
    }
 }; //运行时间28ms

再改进：去掉上一程序中最后一个循环即统计true的个数的循环

class Solution {
    public:
    int countPrimes(int n) {
    if(n<=2){
     return 0; 
    } 
    int fal=0; 
    bool *isPrime=new bool[n];
    memset(isPrime, true, sizeof(bool) * n);//初始化语句，memset() 
   for(int i=3;i*i<n;i+=2){ //i<=sqrt(n);只用判断前sqrt(n)个数，因为大于sqrt(n)的数如果是某个数的因数，则该数一定已经被标记 
      if(isPrime[i]==true){
       for(int j=i*i;j<n;j+=2*i){ //标记i的倍数，从第i*i开始标记，因为比如i=3，i的倍数为6、9、12 ...6为2的倍数，已经被2的倍数标记为false 
        if(isPrime[j]==true){ //为避免重复计算fal,这里要判断if(isPrime[j]==true),因为比如63,被i=3时标记,但i=7时,7*7+2*7=63,如果没有这条判断语句的话,会重复计数 
         fal++; 
         isPrime[j]=false; 
      } 
     } 
    } 
   } 
    int count=n/2-fal; 
    delete[] isPrime; 
    return count; 
    }
   }; //运行时间12ms

看到有0ms的，不知道是真假。

总结：数学分析很重要，算法设计很重要！！！