模板：筛素数法

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