素数线性筛选
2015-07-21 10:21
176 查看
因为线性筛选的原理是利用每个合数的最小质因子去筛选
这是让复杂度达到O(N)的关键一句。原理如下:
因为 i%prime[k] == 0 所以 prime[k] 是 i 的最小质因子
所以 i 可以分解成某个数乘以 prime[k] 即 i = n * prime[k]
我们下一个要筛选的素数是i * prime[k+1] ,然后我们分解 i 得到
n * prime[k] * prime[k+1] , 很明显这个素数可以由较大的数乘以prime[k]筛选掉
还有除了 2, 3 以外只有 6N+1,或者 6N+5 才会是素数。
6N,6N+2, 6N+3, 6N+4分别至少能被2, 2, 3, 2整除,所以不可能是素数。
[code]if(i % prime[k] == 0) break;
这是让复杂度达到O(N)的关键一句。原理如下:
因为 i%prime[k] == 0 所以 prime[k] 是 i 的最小质因子
所以 i 可以分解成某个数乘以 prime[k] 即 i = n * prime[k]
我们下一个要筛选的素数是i * prime[k+1] ,然后我们分解 i 得到
n * prime[k] * prime[k+1] , 很明显这个素数可以由较大的数乘以prime[k]筛选掉
还有除了 2, 3 以外只有 6N+1,或者 6N+5 才会是素数。
6N,6N+2, 6N+3, 6N+4分别至少能被2, 2, 3, 2整除,所以不可能是素数。
[code]int dis[2] = {4, 2}; bool isprime[100000005]; int prime[6000000], top; void Prime() { int i, k, pos = 0; top = 0; prime[top++] = 2; prime[top++] = 3; for(i = 5; i <= 100000000; i += dis[pos ^= 1]){ if(!isprime[i]) prime[top++] = i; for(k = 0; i*prime[k] <= 100000000 && k < top; k++){ isprime[i*prime[k]] = true; if(i % prime[k] == 0) break; } } }
相关文章推荐
- 北京数字认证股份有限公司-数字证书认证管理解决方案
- android常用对话框封装
- SQLServer 数据库镜像+复制切换方案
- 在ubuntu下安装sourceinsight
- 7月18日全球域名注册商(国际域名)保有量及市场份额
- 静态变量的初始化
- 浅析pragma once与include guard
- android从资源文件中读取文件流并显示的方法
- Java面试-Struts2
- #TeX括号#将一种格式的括号变为有左右双引号的区分(白皮书字符串)
- 封装网络请求,下载,上传功能
- 相关列的基数计算
- Java_getProperties以记之
- .netframe work的操作系统支持版本
- 转载priority_queue用法总结
- sudo kill -9 pid
- 如何系统地学习Node.js?
- Unity3d iOS 内购详细流程总汇
- iOS编程第四版第10章 UINavigationController
- C语言之字符串模糊查询方法的实现