1e11以内小于等于这个数的素数的个数

最新推荐文章于 2022-07-09 19:20:44 发布
最新推荐文章于 2022-07-09 19:20:44 发布
阅读量789 收藏
点赞数
分类专栏： 素数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xiaolonggezte/article/details/72810086
版权
hdu 5901

求1e11以内的小于等于这个数的个数：


#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
#include<map>
#include<set>
#include<queue>
using namespace std;
const int maxn = 200 + 10;
#define INF 0x3f3f3f3f
#define clr(x,y) memset(x,y,sizeof x )
typedef long long ll;
#define eps 10e-8
const int Mod = 500 + 7;
typedef pair<int, int> P;

#define MAXN 100
#define MAXM 10001
#define MAXP 40000
#define MAX 400000
//#define clr(ar) memset(ar, 0, sizeof(ar))
#define read() freopen("lol.txt", "r", stdin)
#define dbg(x) cout << #x << " = " << x << endl
#define chkbit(ar, i) (((ar[(i) >> 6]) & (1 << (((i) >> 1) & 31))))
#define setbit(ar, i) (((ar[(i) >> 6]) |= (1 << (((i) >> 1) & 31))))
#define isprime(x) (( (x) && ((x)&1) && (!chkbit(ar, (x)))) || ((x) == 2))

using namespace std;

namespace pcf{
    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;
    }
}

long long solve(long long n){
    int i, j, k, l;
    long long x, y, res = 0;

    for (i = 0; i < pcf::len; i++){
        x = pcf::primes[i], y = n / x;
        if ((x * x) > n) break;
        res += (pcf::Lehmer(y) - pcf::Lehmer(x));
    }

    for (i = 0; i < pcf::len; i++){
        x = pcf::primes[i];
        if ((x * x * x) > n) break;
        res++;
    }

    return res;
}

int main(){
    pcf::init();
    long long n, res;
    int ii = 1;
    while (scanf("%lld", &n) != EOF){
        //res = solve(n);
//        ll l = 2,r = 5e7;
//        while(l < r)
//        {
            cout << l << " " << r << endl;
//            ll mid = (l + r) >> 1;
//            if(n > pcf::Lehmer(mid))
//            {
//                l = mid + 1;
//            }
//            else r = mid;
//        }
        printf("%lld\n",pcf::Lehmer(n));
        //printf("%lld\n", res);
    }
    return 0;
}

1489: L先生与质数V4

Time Limit: 1 Sec   Memory Limit: 16 MB
Submit: 344   Solved: 82
[ Submit][ Status][ Web Board]

Description

在解决了上一个质数问题之后，L先生依然不甘心，他还想计算下更多范围内的质数，你能帮助他吗？（没错这题题面和V3一毛一样）

Input

有多组测试例。（测试例数量<70)
每个测试例一行，输入一个数字n（0<n<=3000000），输入0表示结束。

Output

输出测试例编号和第N个质数。
Case X: Y

Sample Input

1
2
3
4
10
100
0

Sample Output

Case 1: 2
Case 2: 3
Case 3: 5
Case 4: 7
Case 5: 29
Case 6: 541
题意：求出第几个素数，；n <=3e6

思路： 利用上面结论进行二分，找到刚好使得小于等于他的素数个数为n的时候，这个数就是第n个数；

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
#include<map>
#include<set>
#include<queue>
using namespace std;
const int maxn = 200 + 10;
#define INF 0x3f3f3f3f
#define clr(x,y) memset(x,y,sizeof x )
typedef long long ll;
#define eps 10e-8
const int Mod = 500 + 7;
typedef pair<int, int> P;
 
#define MAXN 100
#define MAXM 10001
#define MAXP 40000
#define MAX 400000
//#define clr(ar) memset(ar, 0, sizeof(ar))
#define read() freopen("lol.txt", "r", stdin)
#define dbg(x) cout << #x << " = " << x << endl
#define chkbit(ar, i) (((ar[(i) >> 6]) & (1 << (((i) >> 1) & 31))))
#define setbit(ar, i) (((ar[(i) >> 6]) |= (1 << (((i) >> 1) & 31))))
#define isprime(x) (( (x) && ((x)&1) && (!chkbit(ar, (x)))) || ((x) == 2))
 
using namespace std;
 
namespace pcf{
    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;
    }
}
 
long long solve(long long n){
    int i, j, k, l;
    long long x, y, res = 0;
 
    for (i = 0; i < pcf::len; i++){
        x = pcf::primes[i], y = n / x;
        if ((x * x) > n) break;
        res += (pcf::Lehmer(y) - pcf::Lehmer(x));
    }
 
    for (i = 0; i < pcf::len; i++){
        x = pcf::primes[i];
        if ((x * x * x) > n) break;
        res++;
    }
 
    return res;
}
 
int main(){
    pcf::init();
    long long n, res;
    int ii = 1;
    while (scanf("%lld", &n) != EOF && n){
        //res = solve(n);
        ll l = 2,r = 5e7;
        while(l < r)
        {
//            cout << l << " " << r << endl;
            ll mid = (l + r) >> 1;
            if(n > pcf::Lehmer(mid))
            {
                l = mid + 1;
            }
            else r = mid;
        }
        printf("Case %d: %lld\n",ii ++,l);
        //printf("%lld\n", res);
    }
    return 0;
}
 
/**************************************************************
    Problem: 1489
    User: Pass_by
    Language: C++
    Result: Accepted
    Time:164 ms
    Memory:11296 kb
****************************************************************/




Lazines_by
关注
专栏目录
1e11内质数的个数
百般回首曾经的相遇，留下的只是一缕残忆
09-27 880
#include<cstdio> #include<cmath> using namespace std; #define LL long long const int N = 5e6 + 2; bool np[N]; int prime[N], pi[N]; int getprime() { int cnt = 0; np[0] = np[1] = true; pi[0] = pi[1] = 0; for (int i = 2; i < N
HDU --- 5901 Count primes 【求1e11内的素数模板题】
Anxdada -- 我等风来也等你
09-06 422
传送门 //题意: 求最大1e11内的素数个数. //刚开始做这道题时, 想得可以解决区间1e6以内的素数个数, 然后特判一下输入的数, 分段打表输出答案. 但是死活过不了. 于是搜题解说这个是1e11内的素数个数模板题??? 可以用什么Meisell-Lehmer方法做. 具体证明情况都是给的论文, 我能怎么说, 只能说还是把这个模板留下来, 免得在遇见(手动滑稽). 有两个版本, 一个复杂
1e11内素数总数
xuwnehao的博客
07-09 215
大范围内素数个数
1到n里面素数的个数n<=1e11
I good vegetable a!
09-19 1523
2016沈阳网络赛的最后一题，记录一下这个模板，来自某某大神，并不懂这个的原理，先放着当模板。 // //Created by just_sort 2016/9/12 17:01 //Copyright (c) 2016 just_sort.All Rights Reserved // #include #include #include #include #include
素数筛（1e7内包含4的素数有多少个）
qq_43279324的博客
05-02 838
题目描述 对一个素数，若其含有4，则称其为四素数，如41，149就是四素数，问1e7以内四素数有多少个。 输入描述: 无 输出描述: 1e7以内四素数的个数 用到知识点：素数筛 const int maxnum = 1e7; //prime[]数组用来存素数（下标从1开始） //prime[0]用来存素数的个数 int prime[maxnum]; void getPrime(){ for(int i=2; i<maxnum; i++){ if(prime[i]==0){ prime[
用c语言求1 n的素数个数 给出两种解法,【题目】求n以内的素数个数
weixin_42500739的博客
05-23 7429
最近在leetCode上刷提，还是满锻炼人的，为以后面试打基础吧。不多说下面开始。问题：求[2,n]之间的素数的个数。来源：leetCode OJ提示：Let's start with aisPrimefunction. To determine if a number is prime, we need to check if it is not divisible by any numbe...
求1到n的素数个数C语言,求 1~n 之间素数的个数
weixin_39936558的博客
05-23 6600
1. 筛选法筛选掉偶数，然后比如对于 3，而言，筛选掉其整数倍数；(也即合数一定是某数的整数倍，比如 27 = 3*9)int n = 100000000;bool flag[100000000+1]; // flag[0] 无用的空间；int count() {int cnt = 0;flag[2] = 1;for (int i = 3; i < n; ++...
204  统计小于n的质数的个数（埃氏筛法、欧拉筛法）
liupc的学习笔记
08-05 1239
点击此处返回总目录 一、题目 二、数学基础 1. 质数...
算法竞赛——进阶指南——NOIP2009 acwing200. Hankson的趣味题
bjfu170203101的博客
02-19 352
x一定是b的倍数，d的约数。 由于1e9内的所有数中，约数最多的数 的约数为1536个。 所以我们可以枚举d的约数，判是否成立。 但2000组数据 会T 我们可以从组成约数的本质来枚举，预处理出sqrt（2e9）的质数。 找出所有d的质因子（有可能会有大于sqrt(d)的质因子，但只可能有一个）判断下 。 然后dfs枚举就行 最坏情况2000*1536 #i...
大于或小于100万，1000万，1亿，10亿，1000亿，万亿，亿亿，10亿亿，100亿亿上下的10个质数（素数)...
weixin_33851429的博客
05-12 4381
2019独角兽企业重金招聘Python工程师标准>>> ...
求1e11以内的素数个数
__water
09-18 2158
#include<cstdio> #include<cmath> using namespace std; #define LL long long const int N = 5e6 + 2; bool np[N]; int prime[N], pi[N]; int getprime() { int cnt = 0; np[0] = np[1
求1-1e11内的素数个数(HDU 5901 Count primes )
taotao 的大学墓志
01-14 1703
找了份代码，看不懂，当模板存一波#include <bits/stdc++.h> #define ll long long using namespace std; ll f[340000],g[340000],n; void init(){ ll i,j,m; for(m=1;m*m<=n;++m)f[m]=n/m-1; for(i=1;i<=m;++i)g[i]=i-1;
HDU 5901 Count primes （求1e11内素数个数）
卡尔
09-18 5335
死亡之题…… 比赛打表打了3小时，想的是分段打表一些边界数据，类似这题：把1e9分成5000个20W大小的区间，打表出区间边界（1，200000，400000...）的答案。这样最坏查找最多也是20W-1 结局是三小时也没打出来，在第四小时的时候发现就算表打出来复杂度也在1e11左右，亲妈爆炸！ 最最后的结局是，论文题：Wiki 或者 Lehmer快速求1e11以内质数个数
超高速计算n以内素数个数（百亿内3毫秒解决）
梦想起飞的地方
03-12 7639
判断n以内素数个数有很多算法，最简单的是循环直接判断，这个效率不用说，n稍大就不行了。最流行的是筛选法，原理就是定义一个素数标志位表，初始为1，遇到一个数如果对应标志位为1判断这个数是不是素数，是将该为置1，不是放0，然后将他的倍数位置全部置0，然后继续。。这个效率还是比较快的，但是计算到10^8时候需要3s左右了，对于一般要求基本够了，但是对于ACM里面对时间要求很严还是不够。可以对帅选法进行优
(Relax 数论 1.5)POJ 1811 Prime Test(MillRabin模板题+Pollard模板题:判定大素数&&合数分解)
黄俊东的专栏
12-11 1440
#include #include #include using namespace std; typedef long long LL; const int maxn = 10000; LL factor[maxn];//用来保存质因子 int tot;//用来记录质因子的个数 const int S = 20; LL muti_mod(LL a, LL b, LL c) {
求n很大时求，<= n的素数的个数(时间复杂度o（n ^ 2 / 3）模板
AI_xue
09-19 1252
//Meisell-Lehmer #include #include using namespace std; #define LL long long const int N = 5e6 + 2; bool np[N]; int prime[N], pi[N]; int getprime() { int cnt = 0; np[0] = np[1] = true; pi[
求出小于等于n的所有素数。
肖碧文的博客
02-21 1372
求出小于等于n的所有素数。 import com.sun.org.apache.xerces.internal.impl.xs.SchemaSymbols; import java.util.ArrayList; import java.util.List; public class PrimeNumber { public static void main(String[] args){ ...
NEFU117 素数个数的位数【素数定理】
程序员充电站（itcharge）
02-02 1131
题目大意： 给你一个整数N(1 < N < 1000000000)，如果小于10^N的整数中素数的个数为π(N)， 那么问题来了：求π(N)的位数是多少。 思路： 素数的个数π(N)有素数定理可得：π(N) = N/ln(N)。本题中π(10^N) = 10^N/ln(10^N)。 问题就转换为：求N^10*ln(N^10)共有多少位。设共有x位，可得 10^x = 10^N/ln(10^N)。 对两边同时取对数log10，得： 10^x = N^10 / ln(N^10) l
c++输入一个数n，求所有小于等于这个数的素数
最新发布
07-14
这段代码定义了两个函数，`isPrime`函数用于判断一个数是否为素数，`findPrimes`函数用于找到小于等于给定数n的所有素数。在`main`函数中，先输入一个数n，然后调用`findPrimes`函数找到小于等于n的素数，并输出结果...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值