#include <iostream>
#include <ctime>
using namespace std;
const int N=1e7+5;//定义常量方法一
#define M 1e7+5//定义常量方法二
int a[N];//c++初始化便是0,不需要再初始化。
int b[N];
int c[N];
int prime[N];
int cnt=0;
void eum_sieve(int n)//暴力筛素数
{
int i,j;
for(i=2;i<n;i++)
{
for(j=2;j<=i/j;j++)//试除法
{
if(i%j==0)
break;
}
if(j>i/j)
a[i]=1;//为素数
}
}
void er_sieve(int n)//埃氏筛
{
int i,j;
for(i=2;i<=n/i;i++)
{
if(!b[i])//此处为b[i]==0为素数
{
for(j=i*i;j<n;j+=i)
b[j]=1;//b[j]=1为合数
}
}
}
void el_sieve(int n)//欧拉筛
{
int i,j;
for(i=2;i<n;i++)
{
if(!c[i])//c[i]=0为素数
prime[cnt++]=i;//存储素数i。
for(j=0;j<cnt&&i*prime[j]<n;j++)
{
c[i*prime[j]]=1;//为合数:prime[j]的i倍为合数
if(i%prime[j]==0)//保证每个合数只筛一次,由最小质因数筛出。
break;
}
}
}
int main()
{
int i,cnt1,cnt2,cnt3;
clock_t s1,s2,s3,f1,f2,f3;
cnt1=cnt2=cnt3=0;
s1=clock();//返回cpu当前时钟单元数
eum_sieve(N);
f1=clock();
s2=clock();
er_sieve(N);
f2=clock();
s3=clock();
el_sieve(N);
f3=clock();
for(i=2;i<10000;i++)
{
if(a[i]==1)
{
cnt1++;
cout<<i<<ends;
if((cnt1+1)%10==0)
cout<<endl;
}
}
for(i=2;i<10000;i++)
{
if(b[i]==0)//b[i]=0为素数
{
cnt2++;
cout<<i<<ends;
if((cnt2+1)%10==0)
cout<<endl;
}
}
while(prime[cnt3]<10000)
{
cout<<prime[cnt3++]<<ends;
if((cnt3+1)%10==0)
cout<<endl;
}
cout<<endl<<"cnt1="<<cnt1<<endl;
cout<<"eum_sieve time cost:"<<double(f1-s1)/CLOCKS_PER_SEC<<endl;
cout<<endl<<"cnt2="<<cnt2<<endl;
cout<<"er_sieve time cost:"<<double(f2-s2)/CLOCKS_PER_SEC<<endl;
cout<<endl<<"cnt3="<<cnt3<<endl;
cout<<"el_sieve time cost:"<<double(f3-s3)/CLOCKS_PER_SEC<<endl;
return 0;
}
当N=1e7+5时:
当N=1e8+5时: