Eratosthenes筛法 快速求素数

Eratosthenes筛法，快速求素数。

时间复杂度 O(nlogn)。

---

思想

---

对于每个不超过n的非负整数p，删除2p，3p，4p，…，当处理完所有数后，还没有被删除的就是素数。

---

代码

---

#include <iostream>
#include <ctime>
using namespace std;

/**
 * Eratosthenes筛法，求1~n以内的素数
 */
int* Eratosthenes(int n) {
    int *vis = new int[n+1];
    memset(vis, false, n+1);
    int m = sqrt(n + 0.5);
    
    for(int i=2; i<=m; i++) 
        if(!vis[i])
            for(int j=i*2; j<=n; j+=i) 
                vis[j] = true;
    
    return vis;
}

int main(int argc, char *argv[]) {
    clock_t start,finish; // 定义变量
    
    int n;
    cin >> n;
    start = clock();
    int *a = Eratosthenes(n);
    for(int i=2; i<=n; i++)
        if(a[i] == false) 
            cout << i << ' ';
    cout << endl;
    
    finish = clock();
    double time=(double)(finish-start) / CLOCKS_PER_SEC;//计算运行时间
    printf("Used time %lf s\n", time);//输出运行时间
    
    return 0;
}