模板：篩素數法

標籤：style   blog   http   color   for   ar   art   div   log   

參考:http://blog.csdn.net/liukehua123/article/details/5482854

 

1.開一個大的bool型數組prime[]，大小就是n+1就可以了.先把所有的下標為奇數的標為true,下標為偶數的標為false.

    2.然後：

      for( i=3; i<=sqrt(n); i+=2 )

      {   if(prime[i])

          for( j=i+i; j<=n; j+=i ) prime[j]=false;

      }

    3.最後輸出bool數組中的值為true的單元的下標，就是所求的n以內的素數了。

    原理很簡單，就是當i是質(素)數的時候，i的所有的倍數必然是合數。如果i已經被判斷不是質數了，那麼再找到i後面的質數來把這個質

數的倍數篩掉。

 

Code:

 

 1 #include <cmath> 2  3 #define MAX_NUM 10000 4  5  6 bool * arr_prime = new bool[MAX_NUM + 1]; 7      8 for(i = 3; i <= MAX_NUM; i += 2) arr_prime[i] = true; 9 for(i = 4; i <= MAX_NUM; i += 2) arr_prime[i] = false;10 arr_prime[2] = true;11 12 int sqrt_mn = sqrt(MAX_NUM);13 14 for(i = 3; i < sqrt_mn; i += 2)15 {16     if(arr_prime[i])17     {18         for(j = i + i; j <= MAX_NUM; j += i) arr_prime[j] = false;19     }20 }21     

 

模板：篩素數法