[学习小记] [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