素数（Prime）

素数的判断：

 #include<math.h>
 bool IsPrime(int n)
 {
     if(n <= )
         return false;
     int sqr = (int)sqrt(1.0*n);
     for(int i=; i<=sqr; i++)
         if(n%i==)  return false;
     return true;
 }

获取素数表：

 const int M=;
 vector<int> prime;
 bool isPrime[M]={};

 void FindPrime()
 {
     for(int i=; i<M; i++)
         if(IsPrime(i))
         {
             prime.push_back(i);
             isPrime[i] = ;
         }
 }

埃氏筛法（Eratosthenes）：

 #include<cstdio>
 #include<algorithm>
 using namespace std;
 const int M=;
 int isPrime[M];

 void Eratosthenes()
 {
     fill(isPrime, isPrime+M, );
     isPrime[]=;
     isPrime[]=;
     for(int i=; i<M; i++)
         if(isPrime[i] == )
         {
             printf("%d ", i);
             for(int j=*i; j<M; j+=i)
                 isPrime[j] = ;
         }

 }

 int main()
 {
     Eratosthenes();
     return ;
 }

• 时间复杂度：O（ Nlog(logN) ）
• 第一个测出地球周长的男人-。-

欧氏筛法（Euler）：

 #include<cstdio>
 #include<vector>
 using namespace std;
 const int M=;
 int isPrime[M];

 void Euler()
 {
     fill(isPrime, isPrime+M, );
     isPrime[]=;
     isPrime[]=;
     vector<int> prime;
     for(int i=; i<M; i++)
     {
         if(isPrime[i])
         {
             prime.push_back(i);
             printf("%d ", i);
         }
         for(int j=; j<prime.size(); j++)
         {
             if(i*prime[j] > M)
                 break;
             isPrime[i*prime[j]] = ;
             if(i%prime[j] == )
                 break;
         }
     }
 }

 int main()
 {
     Euler();
     return ;
 }

• 时间复杂度：O（N）
• 欧拉待过的俄法德，都曾是世界最强的国家，如果当初欧拉来中国，是不是咱们就年年ACM总冠军了-。-