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
    评论
根据引用\[1\]和引用\[2\]的内容,题目P1835要求统计一个范围内的素数个数。这个范围可以在int的任意范围内,但区间的大小小于1000000。为了解决这个问题,可以使用线性法来预处理素数列表,并用这些质因数标记\[L,R\]范围内的合数。 首先,使用线性法来预处理素数列表。具体步骤如下: 1. 初始化isprime数组,将所有元素标记为1,表示都是素数。 2. 从2开始遍历到50000,如果isprime\[i\]为1,则将i加入prime数组,并将isprime\[i\]标记为0。 3. 对于每个素数prime\[i\],从1开始遍历到now,并且i*prime\[j\]小于50000,将isprime\[i*prime\[j\]\]标记为0。如果i能整除prime\[j\],则跳出内层循环。 接下来,使用res数组来标记\[L,R\]范围内的合数。具体步骤如下: 1. 初始化res数组,将所有元素标记为1,表示都是素数。 2. 对于每个素数prime\[i\],从max(2,(l-1)/prime\[i\]+1)开始遍历到r,将res\[j*prime\[i\]-l\]标记为0。这里的j至少为2才有意义,否则j=0或者1的时候会出错。同时,要注意(long)j*prime\[i\]<=r,因为如果都是int会溢出。 3. 如果l<=1,将res\[1-l\]标记为0。 4. 遍历res数组,统计标记为1的个数,即为素数的个数。 根据以上步骤,可以得到解决题目P1835的代码。 #### 引用[.reference_title] - *1* *2* [【洛谷】P1835 素数密度](https://blog.csdn.net/u014156276/article/details/108294974)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值