埃氏筛法(快速筛选n以内素数的个数)

 1 /*
 2     |埃式筛法|
 3     |快速筛选素数|
　　　　|15-7-26|
 4 */
 5 #include <iostream>
 6 #include <cstdio>
 7 using namespace std;
 8 const int SIZE = 1e7;
 9 
10 int prime[SIZE];            // 第i个素数
11 bool is_prime[SIZE];    //true表示i是素数
12 
13 int slove(int n)
14 {
15     int p = 0;
16     for(int i = 0; i <= n; i++)
17         is_prime[i] = true;             //初始化
18     is_prime[0] = is_prime[1] = false;      //0,1不是素数
19     for(int i = 2; i <= n; i++)
20     {
21         if(is_prime[i])                //这里比较巧妙, 我只是意会
22         {
23             prime[p++] = i;             //计算素数的个数,也记录下了素数
24             for(int j = 2 * i; j <= n; j += i)      // 除掉了i的倍数的数字
25                 is_prime[j] = false;
26         }
27     }
28     return p;
29 }
30 
31 int main()
32 {
33     int n;
34     while(cin >> n)
35     {
36         int res = slove(n);
37         cout << res << endl;
38         for(int i = 0; i < res; i++)
39             cout << prime[i] << endl;
40     }
41 }