Count primes (Meissel-Lehmer算法)
来源:互联网 发布:单端口多用户什么意思 编辑:程序博客网 时间:2024/06/05 23:48
Easy question! Calculate how many primes between [1…n]!
Input
Each line contain one integer n(1 <= n <= 1e11).Process to end of file.
Output
For each case, output the number of primes in interval [1…n]
Sample Input
2
3
10
Sample Output
1
2
4
[分析]
模板,有必要记录一下。
[代码]
#include <iostream>#include <cstdio>#include <algorithm>#include <cstring>#include <vector>#include <cmath>#define MAXN 100 // pre-calc max n for phi(m, n)#define MAXM 10010 // pre-calc max m for phi(m, n)#define MAXP 40000 // max primes counter#define MAX 400010 // max prime#define setbit(ar, i) (((ar[(i) >> 6]) |= (1 << (((i) >> 1) & 31))))#define chkbit(ar, i) (((ar[(i) >> 6]) & (1 << (((i) >> 1) & 31))))#define isprime(x) (( (x) && ((x)&1) && (!chkbit(ar, (x)))) || ((x) == 2))using namespace std; long long dp[MAXN][MAXM]; unsigned int ar[(MAX >> 6) + 5] = { 0 }; int len = 0, primes[MAXP], counter[MAX]; void Sieve() { setbit(ar, 0), setbit(ar, 1); for (int i = 3; (i * i) < MAX; i++, i++) { if (!chkbit(ar, i)) { int k = i << 1; for (int j = (i * i); j < MAX; j += k) setbit(ar, j); } } for (int i = 1; i < MAX; i++) { counter[i] = counter[i - 1]; if (isprime(i)) primes[len++] = i, counter[i]++; } } void init() { Sieve(); for (int n = 0; n < MAXN; n++) { for (int m = 0; m < MAXM; m++) { if (!n) dp[n][m] = m; else dp[n][m] = dp[n - 1][m] - dp[n - 1][m / primes[n - 1]]; } } } long long phi(long long m, int n) { if (n == 0) return m; if (primes[n - 1] >= m) return 1; if (m < MAXM && n < MAXN) return dp[n][m]; return phi(m, n - 1) - phi(m / primes[n - 1], n - 1); } long long Lehmer(long long m) { if (m < MAX) return counter[m]; long long w, res = 0; int i, a, s, c, x, y; s = sqrt(0.9 + m), y = c = cbrt(0.9 + m); a = counter[y], res = phi(m, a) + a - 1; for (i = a; primes[i] <= s; i++) res = res - Lehmer(m / primes[i]) + Lehmer(primes[i]) - 1; return res; }int main(){ init(); long long int n; while(scanf("%lld",&n)!=EOF) { printf("%lld\n",Lehmer(n)); } return 0;}
阅读全文
0 0
- Count primes (Meissel-Lehmer算法)
- HDU 5901 Count primes (Meisell-Lehmer算法 -- 模板)
- HDU5901 Count primes(Meisell-Lehmer算法+模板)
- HDU 5901 Count primes 2016年沈阳网络赛 (Lehmer素数计数)
- hdu 5901 Count primes(Meisell-Lehmer 统计n(很大)以内的素数个数)
- [hdu5901 Count primes]Meisell-Lehmer求质数个数PI(X)
- 算法练习(29):Count Primes
- Meissel Lehmer Algorithm快速统计前n个元素中所有的质数以及改编(论文题)
- 【HDU】5901】【模板题】Count primes 【Meisell-Lehmer求质数个数】
- Leetcode(204) Count Primes
- Leetcode(59)-Count Primes
- 204. Count Primes(unsolved)
- 204. Count Primes(埃拉托色尼)
- Count primes (模板题)
- 【算法作业7】LeetCode 204. Count Primes
- 15算法课程 204. Count Primes
- Count Primes
- Count Primes
- chrome developer tool 调试技巧
- timezone
- DSP基础知识
- Tomcat-----Tomcat设置启动参数
- 编译caffe时遇到的undefined reference to的问题
- Count primes (Meissel-Lehmer算法)
- PHP 知识点
- 最大子数组问题--分治法的思想
- CSS3自定义滚动条样式 -webkit-scrollbar
- ORACLE中用rownum分页并排序的SQL语句
- 关键词多空格处理
- JAVA通过JDBC链接数据库获取数据(三)PreparedStatement预编译SQL
- 【007】Java语言(环境变量配置)
- 我是如何对网站CSS进行架构的