[模板][洛谷P1835]素数密度(区间筛)

题目←

每个数n的最小质因子一定不超过 n
预处理出 m 内的素数表,删去每个素数在n-m区间内所有倍数,剩下的一定全为素数
del数组采用下标偏移,i实际表示数n+i是否已被删去
del数组大小为0~m-n

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#define LL long long
using namespace std;
const LL MAXN = 1000000 + 50;
LL prime[MAXN];
bool sh[MAXN],del[MAXN];
int cnt;
void solve1(LL n){
    for(LL i = 2;i <= n;i ++){
        if(!sh[i])prime[++ cnt] = i;
        for(LL j = 1;j <= cnt&&i*prime[j] <= n;j ++){
            sh[i*prime[j]] = true;
            if(i%prime[j] == 0)break;
        }
    }
}
int tot;
void solve2(LL n,LL m){
    for(LL i = 1;i <= cnt;i ++){
        LL j = n/prime[i];
        if(j*prime[i] < n)j ++;
        if(j == 1)j ++;
        for(j;j*prime[i] <= m;j ++){
            if(!del[j*prime[i] - n])tot ++;
            del[j*prime[i] - n] = true;
        }
    }
}
LL n,m;
int main(){
    scanf("%lld%lld",&n,&m);
    LL t = sqrt(m) + 1;
    solve1(t);
    solve2(n,m);
    printf("%d",m - n + 1 - tot);
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值