初学C语言,看到求素数的代码,不是太满意,自己设计了一个计算流程图.在贴吧里也问了一下,和"埃拉托色尼(Eratosthenes)筛法”有些撞车,当然在下不敢与数学巨匠们争长短,不过还是想自己试一下.
百度了一下有人已经用这类方法做了代码,不过和我的想法有出入。
思路：
1要避免冗余的求余计算的计算次数,这样在数字足够大的时候更能体现更高的计算效率,
1)分母排除偶数,可以减少一半的计算量.具体做法是改+1递增为从第二个素数3开始+2递增;
2)筛选求余计算的分子,而根据素数的定义,只要不被小于分母的素数整除即可,?设置最大求余分子计算至分母平方根+1;如此可以避免对每个数进行求余,可以减少大量的冗余计算量.
3)快速淘汰合数,从最小素数开始计算,能更快更有效率
2 找出一个素数,然后赋值,可以避免存储无用的数值,利用找出的素数作为分子进行计算,可以进一步提高效率.
经过2个晚上的摸索,终于成功了,(初学者就这样吧,还是不要要求太多了......)
下面是代码：
——————————————————————————
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <conio.h>
int main()
{
int a[1001];//数组
int N;//需要输出素数的个数
int M;//计算中间数
printf("亲,要多少个素数呢?不要大于1000哦.\n");
scanf("%d",&N);
a[0]=2;//赋值给两个特殊的质数
a[1]=3;
int m1=1;//循环2控制变量
int M1;//M1:求余分子最大数
int j;
/*-------------------循环体---------------------*/
for(int i=1;i<N;i++)//for1:从最小以此找出素数
{
m1=1;
for(M=a[i]+2;m1;M=M+2)//for2:分母+2直至符合素数条件,分子
{
M1=(int)sqrt(M)+1;//计算求余运算的最大分子
for(j=1;j<=i;j++)//for3:从最小的素数开始求余
{
if(a[j]<M1)//判断求余分子是否足够大
{//不够大时,计算是否能被整除
if(M%a[j]==0)//如果被整除
{
j=i+1;//促使for3 循环结束
}
else//不能被整除
{
}
}
else
{
a[i+1]=M;//足够大赋值给数列
j=i+1;//for3跳出
m1=0;//for2跳出
}
}
}
}
for(int iii=0;iii<N;iii++){printf("%7d,",a[iii]);}
printf("\n");
getch();
return 0;
}