码迷,mamicode.com
首页 > 编程语言 > 详细

C语言最优化求素数

时间:2015-06-16 22:57:23      阅读:205      评论:0      收藏:0      [点我收藏+]

标签:

   初学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;

}

C语言最优化求素数

标签:

原文地址:http://www.cnblogs.com/bfoyun/p/4581973.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!