筛素数

#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
const int maxn = 10000000;
const int maxp = 700000;
const int MOD = 1e9;
typedef long long LL;
int vis[maxn]; //vis[i] = 1是合数、vis[i] = 0是素数或者1
int prime[maxp];
// 筛素数
void sieve(int n){
    int m = (int)sqrt(n + 0.5); // 避免浮点误差
    memset(vis,0,sizeof vis);
    for(int i = 2 ; i <= m ; i++) if(!vis[i])
        for(int j = i * i ; j <= n ; j += i) vis[j] = 1;
}
// 生成素数表，放在prime数组中，返回素数的个数
int gen_primes(int n){
    sieve(n);
    int c = 0;
    for(int i = 2 ; i <= n ; i++)if(!vis[i])
        prime[c++] = i;
    return c;
}
int main(){
    int c = gen_primes(100);
    for(int i = 0 ; i < c ; i++)
        printf("%d
",prime[i]);
    return 0;
}