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

[算法]素数筛法

时间：2014-12-10 16:26:08      阅读：221      评论：0      收藏：0      [点我收藏+]

标签：算法   素数   

【方法一】

bubuko.com,布布扣

bubuko.com,布布扣

【代码一】

[cpp] view plaincopy
  1. //判断是否是一个素数  
  2. int IsPrime(int a){  
  3.     //0,1，负数都是非素数  
  4.     if(a <= 1){  
  5.         return 0;  
  6.     }  
  7.     //计算枚举上界，为防止double值带来的精度损失，所以采用根号值取整后再加1，即宁愿多枚举一个，也不愿少枚举一个数  
  8.     int bound = (int)sqrt(a) + 1;  
  9.     for(int i = 2;i < bound;i++){  
  10.         //依次枚举这些数能否整除x，若能则必不是素数  
  11.         if(a % i == 0){  
  12.             return 0;  
  13.         }  
  14.     }  
  15.     return 1;  
  16. }  

bubuko.com,布布扣

【方法二】

bubuko.com,布布扣

【代码二】

[cpp] view plaincopy
  1. #define MAXSIZE 10001  
  2.   
  3. int Mark[MAXSIZE];  
  4. int prime[MAXSIZE];  
  5.   
  6. //判断是否是一个素数  Mark 标记数组 index 素数个数  
  7. int Prime(){  
  8.     int index = 0;  
  9.     memset(Mark,0,sizeof(Mark));  
  10.     for(int i = 0;i < MAXSIZE;i++){  
  11.         //已被标记  
  12.         if(Mark[i] == 1){  
  13.             continue;  
  14.         }  
  15.         else{  
  16.             //否则得到一个素数  
  17.             prime[index++] = i;  
  18.             //标记该素数的倍数为非素数  
  19.             for(int j = i*i;j < MAXSIZE;j += i){  
  20.                 Mark[j] = 1;  
  21.             }  
  22.         }  
  23.     }  
  24.     return index;  
  25. }  

【方法三】

这种方法比较好理解，初始时，假设全部都是素数，当找到一个素数时，显然这个素数乘上另外一个数之后都是合数

把这些合数都筛掉，即算法名字的由来。但仔细分析能发现，这种方法会造成重复筛除合数，影响效率。

比如10，在i=2的时候，k=2*15筛了一次；在i=5，k=5*6 的时候又筛了一次。所以，也就有了快速线性筛法。

【代码三】

[cpp] view plaincopy
  1. int Mark[MAXSIZE];  
  2. int prime[MAXSIZE];  
  3.   
  4. //判断是否是一个素数  Mark 标记数组 index 素数个数  
  5. int Prime(){  
  6.     int index = 0;  
  7.     memset(Mark,0,sizeof(Mark));  
  8.     for(int i = 2; i < MAXSIZE; i++)  
  9.     {  
  10.         //如果未标记则得到一个素数  
  11.         if(Mark[i] == 0){  
  12.             prime[index++] = i;  
  13.         }  
  14.         //标记目前得到的素数的i倍为非素数  
  15.         for(int j = 0; j < index && prime[j] * i < MAXSIZE; j++)  
  16.         {  
  17.             Mark[i * prime[j]] = 1;  
  18.             if(i % prime[j] == 0){  
  19.                 break;  
  20.             }  
  21.         }  
  22.     }  
  23.     return index;  
  24. }  
利用了每个合数必有一个最小素因子。每个合数仅被它的最小素因子筛去正好一次。所以为线性时间。
代码中体现在：
if(i%prime[j]==0)break;
prime数组 中的素数是递增的,当 i 能整除 prime[j]，那么 i*prime[j+1] 这个合数肯定被 prime[j] 乘以某个数筛掉。
因为i中含有prime[j], prime[j] 比 prime[j+1] 小。接下去的素数同理。所以不用筛下去了。

在满足i%prme[j]==0这个条件之前以及第一次满足改条件时,pr[j]必定是pr[j]*i的最小因子。





参考博文：点击打开链接


习题练习：点击打开链接







[算法]素数筛法

标签：算法   素数   

原文地址：http://blog.csdn.net/sunnyyoona/article/details/41846513
踩
(0)
赞
(0)
   
举报
评论 一句话评论（0）
登录后才能评论！
分享档案
更多>
2021年07月29日 (22)
2021年07月28日 (40)
2021年07月27日 (32)
2021年07月26日 (79)
2021年07月23日 (29)
2021年07月22日 (30)
2021年07月21日 (42)
2021年07月20日 (16)
2021年07月19日 (90)
2021年07月16日 (35)
周排行
mamicode.com排行更多图片
更多
友情链接
兰亭集智   国之画   百度统计   站长统计   阿里云   chrome插件   新版天听网
关于我们 - 联系我们 - 留言反馈
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码！