高效素数打表

大家所知的素数打表时间复杂度几乎都是n2。

就是这种：

void init_prime() { int i, j; for(i = 2;i <= sqrt(1000002.0); ++i) { if(!prime[i]) for(j = i * i; j < 1000002; j += i) prime[j] = 1; } j = 0; for(i = 2;i <= 1000002; ++i) if(!prime[i]) prime[j++] = i; }

在网上搜到了另一种方法，效率提高了不少。巧妙之处还在研究中，等理解差不多就仔细分析一下。

下面是代码，大家可以发表一下自己的看法，交流交流。

代码如下：

/* 遇到素数需要打表时，先估算素数的个数： num = n / lnx; num为大概数字，越大误差越小（只是估计，用于估算素数表数组大小） 这个打表法效率貌似很高，网上说几乎达到了线性时间（不知道是真是假=。=） */ #include<iostream> #include<cstring> #include<cstdio> #include<ctime> #include<algorithm> using namespace std; int n; bool visit[10100000]; int prime[10000000]; void init_prim() { memset(visit, true, sizeof(visit)); int num = 0; for (int i = 2; i <= n; ++i) { if (visit[i] == true) { num++; prime[num] = i; } for (int j = 1; ((j <= num) && (i * prime[j] <= n)); ++j) { visit[i * prime[j]] = false; if (i % prime[j] == 0) break; //点睛之笔 } } } int main() { memset(prime, 0, sizeof(prime)); int count = 0; cin>>n; init_prim(); for(int i = 0; i <= n; ++i) if(prime[i]) { cout<<prime[i]<<" "; count++; } cout<<endl; cout<<"素数个数为："<<count<<endl; }