[JZOJ] 5835. Prime

weixin_33725807

于 2018-10-19 16:39:00 发布

阅读量110

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/ghostcai/p/9817281.html

本文介绍了一种利用区间筛法解决特定数论问题的算法实现。通过预处理素数并标记合数，算法能在给定区间内快速找出满足特定条件的数。文章详细展示了算法流程，包括素数生成、区间筛选及最终结果计算。

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

显然的区间筛法

#include<algorithm> 
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;

typedef long long ll;

ll L,R,k,n;

const int MAXN = 10000005;

int prime[MAXN],tot;
bool isx[MAXN];
void mkprime(){
    for(int i=2;i<=n;i++){
        if(!isx[i]){prime[++tot]=i;}
        for(int j=1;i*prime[j]<=n&&j<=tot;j++){
            int v=prime[j];
            isx[i*v]=1;
            if(i%v==0)break;    
        }
    }
}

bool vis[MAXN];
int main(){
    freopen("prime.in","r",stdin);
    freopen("prime.out","w",stdout);
    cin>>L>>R>>k;
    n=sqrt(R);
    n=min(n,k);
    mkprime();
    for(int i=1;i<=tot;i++){
        int v=prime[i];
        ll cur=1ll*v*(L/v);
        if(cur<L)cur+=v;
        if(cur==v)cur+=v;
        for(;cur<=R;cur+=v)vis[cur-L]=1;    
    }
    ll ans=0;
    int up=R-L,i;
    for(i=0;i<=up;i++){
        if(!vis[i]) ans^=(i+L);
    }
    cout<<ans;
    return 0;
}

转载于:https://www.cnblogs.com/ghostcai/p/9817281.html