博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
线性素数筛(欧拉筛)(超级好的MuBan)
阅读量:4961 次
发布时间:2019-06-12

本文共 1057 字,大约阅读时间需要 3 分钟。

Problem:找出小于等于n的所有素数的个数。

#include 
using namespace std;const int maxn = 1e6;int prime[maxn]; // 欧拉线性素数筛,O(n)bool vis[maxn]; // 标记int Prime(int n) { memset(vis,false,sizeof(vis)); int cnt = 0; vis[0] = vis[1] = true; for(int i = 2; i <= n; i ++) { if(!vis[i])prime[cnt++] = i; for(int j = 0; j < cnt && i*prime[j] <= n; j ++) { vis[i*prime[j]] = true; if(!(i%prime[j])) break; } } return cnt;}int main(){ int n; cin >> n; int ans = 0; ans = Prime(n); cout << ans << endl; return 0;}

if(i % prime[j] == 0) break;

解释:

      首先,任何合数都能表示成多个素数的积。所以,任何的合数肯定有一个最小质因子。我们通过这个最小质因子就可以判断什么时候不用继续筛下去了。

      当i是prime[j]的整数倍时(i % prime[j] == 0),i*prime[j+1]肯定被筛过,跳出循环。

      因为i可以看做prime[j]*某个数, i*prime[j+1]就可以看做 prime[j]*某个数*prime[j+1] 。而 prime[j] 必定小于 prime[j+1],

所以 i*prime[j+1] 必定已经被 prime[j]*某个数 筛掉,就不用再做了√

      同时我们可以发现在满足程序里的两个条件的时候,prime[j]必定是prime[j]*i的最小质因子。这个性质在某些题里可以用到。

      这样就可以在线性时间内找到素数啦~\(≧▽≦)/~

解释转自

转载于:https://www.cnblogs.com/lcchy/p/10139546.html

你可能感兴趣的文章
在SQLServer 2005附加SQLServer 2008数据库异常处理
查看>>
网易新闻侧滑抽屉效果(利用父子控制器实现)
查看>>
Ceph:pg peering过程分析
查看>>
4.高级特性(1)
查看>>
table添加正确的样式
查看>>
求斐波那契数列的第n项
查看>>
一些程序的整理
查看>>
Ubuntu 10.04 compiz的安装
查看>>
angular项目引用第三方公共js文件
查看>>
内置函数
查看>>
JAVA时区处理(转)
查看>>
sqlserver 字符串相关函数
查看>>
C#_02.10_基础一_.NET框架
查看>>
apk编译/反编译工具baksmali… 分类: Android安装及...
查看>>
递归趣文【Matrix67】
查看>>
转载: android 学习架构
查看>>
POJ-3176 Cow Bowling
查看>>
JAVA中保留小数的多种方法
查看>>
TCP的拥塞控制 (三)
查看>>
Windows Phone7页面导航
查看>>