求小于N的素数个数(N>1e9)

本文介绍了一种高效的计算小于等于给定数值的素数个数的方法,使用了动态规划和递归策略,通过预处理和缓存中间结果来加速计算过程。文章详细解释了算法的设计思路和实现细节,并提供了完整的C++代码实现。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

#include<bits/stdc++.h>
#define il inline
#define pb push_back
#define ms(_data,v) memset(_data,v,sizeof(_data))
#define SZ(a) int((a).size())
using namespace std;
typedef long long ll;
const ll inf=0x3f3f3f3f;
const int N=320005;
//il int Add(int &x,ll y) {return x=x+y>=mod?x+y-mod:x+y;}
//il int Mul(int &x,ll y) {return x=x*y>=mod?x*y%mod:x*y;}
ll phi[10005][105], p2[N], ans[N];
int len, vis[N];
void init() {
	len = 0;
	for(int i=2; i<N; i++) {
		if(!vis[i]) {
			for(int j=i+i; j<N; j+=i) vis[j]=1;
			p2[len++] = i;
			ans[i] = ans[i-1]+1;
			continue;
		}
		ans[i] = ans[i-1];
	}
	for(int i=0; i<=10000; i++) {
		phi[i][0] = (ll)i;
		for(int j=1; j<=100; j++) {
			phi[i][j] = phi[i][j-1] - phi[i/p2[j-1]][j-1];
		}
	}
}
ll solve_phi(ll m, ll n) {
	if(!n) return m;
	if(p2[n - 1] >= m) return 1;
	if(m<=10000 && n<=100) return phi[m][n];
	return solve_phi(m, n-1) - solve_phi(m/p2[n-1], n-1);
}
ll solve_p2(ll m) {
	if(m < (ll)N) return ans[m];
	ll y = (int)cbrt(m*1.0);
	ll n = ans[y];
	ll sum = solve_phi(m, n) + n -1;
	for(ll i=n; p2[i]*p2[i]<=m; i++)
		sum = sum - solve_p2(m/p2[i])+solve_p2(p2[i])-1;
	return sum;
}
int main(){
    init();
    scanf("%lld", &n);
    printf("%lld\n",solve_p2(n));
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值