P1835 素数密度(大数区间筛)

素数密度 - 洛谷https://www.luogu.com.cn/problem/P1835

注意一点求合数的时候只能从大于等于l的开始

#include <iostream>
#include <cstdio>
#include <string>
#include <algorithm>
#include <vector>
#include <queue>
#include <stack>
#include <cstring>
#include <set>
#include <unordered_map>
#include <cmath>
#include <map>
#include <cctype>
#include <cstdlib>
#include <deque>
typedef long long ll;
typedef unsigned long long ull;
using namespace std;
const int MN = 65005;
const int MAXN = 1e6 + 5;
const int INF = 0x3f3f3f3f;
#define IOS ios::sync_with_stdio(false)
#define lowbit(x) ((x)&(-x))

int prime[MAXN];
bool isp[MAXN];
int cnt;

inline void init() {
	for (int i = 2; i <= 50000; i++) {
		if (!isp[i]) {
			prime[++cnt] = i;
			for (int j = 1; j <= cnt && i * prime[j] <= 50000; j++) {
				isp[prime[j]*i] = true;
				if (i % prime[j] == 0)
					break;
			}
		}
	}
}


int main() {
	ll l, r;
	init();
	scanf("%lld %lld", &l, &r);
	l=l==1?2:l;
	memset(isp, 0, sizeof(isp));
	for (int i = 1; i <= cnt; i++) {
		ll p = prime[i];
		ll st=(l+p-1)/p*p>2*p?(l+p-1)/p*p:2*p;
		for (ll j = st; j <= r; j += p) {
			isp[j - l+1] = true;
		}
	}
	int ans = 0;
	for (int i = 1; i <= r - l+1; i++) {
		if (!isp[i])
			ans++;
	}
	printf("%d", ans);
	return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值