模板：筛素数法

参考: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:

 

#include <cmath>

#define MAX_NUM 10000


bool * arr_prime = new bool[MAX_NUM + 1];
    
for(i = 3; i <= MAX_NUM; i += 2) arr_prime[i] = true;
for(i = 4; i <= MAX_NUM; i += 2) arr_prime[i] = false;
arr_prime[2] = true;

int sqrt_mn = sqrt(MAX_NUM);

for(i = 3; i < sqrt_mn; i += 2)
{
    if(arr_prime[i])
    {
        for(j = i + i; j <= MAX_NUM; j += i) arr_prime[j] = false;
    }
}