素数/质数的三种筛法:暴力、埃氏、欧拉

#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时:
在这里插入图片描述

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值