[学习小记] [C] [算法] 最快的欧拉筛法与素数(复数)打印

3 篇文章 0 订阅
2 篇文章 0 订阅

[学习小记] [C] [算法] 欧拉筛与素数筛选

偶然看到了老师的期末考试题有要求打印指定范围的素数
并且按时间复杂度给分,所以在网上搜罗一天并学习后给出以下代码
此代码因该是线性时间复杂度,可能是最快的方式
在Windows下调试,时间计数器可能有问题,暂不解决
如有错误请指出

#include <stdio.h>
#include <stdbool.h>
#include <stdlib.h>
#include <time.h>

#define MaxArray 100000001                  //首先我们定义数组的最大范围,也就是可以保存多少个素数        
int Prime[MaxArray], PrimeCounter = 0;      //定义一个有MaxArray大小的数组来记录素数,再定义一个计数器来记录已经筛选出了多少个素数
bool Composite[MaxArray] = {false};         //定义一个布尔类型的集合来记录已经筛选出的合数,默认所有的数为素数

void eulerSieve (int N)                     //定义一个子程序,N是筛选的最大范围
{

    for (int NumCounter = 2; NumCounter <= N; NumCounter++)     定义一个自然数计数器,大循环,每次循环令自然数计数器加一,必须从二开始
    {
        if (Composite[NumCounter] == false)                    //判断一个数是否已经被标记为合数
            Prime[PrimeCounter++] = NumCounter;                //如果没有,则把这个数写入素数数组,并且素数计数器加一
        
        for (int Poller = 0; Poller < PrimeCounter; Poller++)     //再定义一个轮询器,第二循环,每次都从0开始轮询判断当前大循环的自然数,轮询参数传给素数组
        {
            if (NumCounter * Prime[Poller] > N)               //如果当前自然数乘以当前轮询素数大于定义的范围,则结束循环,避免无效操作
                break;

            Composite[NumCounter * Prime[Poller]] = true;    //标记一次当前素数的倍数为合数,若无上一步,此处被标记的数字将可能会大于N

            if (NumCounter % Prime[Poller] == 0)             //欧拉筛选核心步骤,判断当前数是否能被已经记录的素数整除,如果可以整除则是合数
                break;
        }
    
    }

}

int main ()
{
    int N; double Times;
    printf ("请输入素数判断的最大值,只能为正整数:");
    scanf ("%d", &N);
    clock_t start, end;
    start = clock ();
    eulerSieve (N);
    end = clock ();
    Times = (double)(end-start)/CLOCKS_PER_SEC;

    for (int Num = 0, Row = 1; Num < PrimeCounter; Num++, Row++)
    {
        printf ("%d\t", Prime[Num]);

        if (Row %10 == 0)
            {
                printf ("\n");
            }
    }

    printf ("\n计算%d位花费时间:%.50f秒\nWritten by Aloento, Bye.\n", N, Times);
    system("pause");
    return 0;
}
//Written by Aloento.

Written by Aloento - SoarTeam

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值