素数打表比较快的方法是利用素数的倍数不是素数的性质来找的,算法将近为线性时间O(n)
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
第一遍先删去二的倍数
3 5 7 9 11 13 15 17 19 21 23 25 27 29
再删去3的倍数
2 3 5 7 11 13 17 19 23 29
再删5的倍数
2 3 5 7 11 13 17 19 23 29
因为7*7=49>30了所以我们到7既可以了
那么剩下的都是30以内的素数了
const int maxn=1e6+5;
bool vis[101000000];
int prime[maxn];
void init()
{
memset(vis,true,sizeof(vis));
int num=0;
for(int i=2;i<maxn;i++)
{
if(vis[i]==true)
{
prime[num++]=i;
}
for(int j=0;(j<=num)&&(i*prime[j]<maxn);j++)
{
vis[i*prime[j]]=false;
if(i%prime[j]==0) break;
}
}
}