质数筛

 

 

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define ll long long
using namespace std;
const int mx = 1e5;
int prim[20000];
bool vis[mx];

void da(){
  memset(vis,0,sizeof(vis));
  int tol = 0;
  for(int i = 2; i < mx ; i++){
    if(!vis[i])
      prim[tol++] = i;
    for(int k = 0; k < tol && prim[k]*i < mx; k++){
    //  if(vis[prim[k]*i] ==1)
    //    cout<<"zha"<<endl;
      vis[prim[k]*i] =1;
      if(i%prim[k] == 0)
        break;
    }
  }
}

int main(){
  da();
for(int i = 0; i < 10 ; i++)
  cout<<prim[i]<<endl;
  return 0;
}

 

 

比较省空间的方法

 

/*=================================*\ 
  素数打表 
  该函数执行后在prim[]数组中存入 
  从2开始的从小到大的numOfPrim个素数 
\*=================================*/  
const int numOfPrim = 1000;  
int prim[numOfPrim] = {2,3};  
void prime(){  
    int tally=2;  
    bool flag;  
    for(int i=5 ; tally<numOfPrim ; i+=2){  
        flag = true;  
        for(int j=0 ; prim[j]*prim[j]<=i ; j++)  
            if( i%prim[j]==0 ){ flag = false; break; }  
        if( flag ){  
            prim[tally]=i;  
            tally++;  
        }  
    }  
}