大家所知的素数打表时间复杂度几乎都是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;
}