其他分享
首页 > 其他分享> > M - Help Hanzo LightOJ - 1197 (大区间素数筛法)

M - Help Hanzo LightOJ - 1197 (大区间素数筛法)

作者:互联网

题解:素数区间问题。注意到a和b的范围是1<<31,所以直接暴力打表肯定不可以。如果一个数是合数,他的两个因子要么是两个sqrt(x),要么就分布在sqrt(x)两端,所以我们可以根据sqrt(n)之前的数来把sqrt(n)之后的素数给筛出来。

首先进行1e6的素数打表。然后对每个l,r,首先找到第一个大于等于l的数,然后根据刚才的筛出来素数,一直累加,直到r为止,这样就可以把因子在1e6范围的合数给标记上了。注意我们保存的时候数组中要减去l,这样数组只要开到1e5就可以了。注意当l=1时,1不是素数,要特判一下。

code:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll N=1E6+7;
bool primes[N];
ll  pre[N];
ll pos=0;
ll mp[100000+7];
void inint(){
    ll maxxn=1e6;
    primes[1]=1;
    primes[2]=0;
    for(ll i=2;i<=maxxn;i++){
        if(!primes[i]) pre[++pos]=i;
        for(ll j=1;j<=pos&&i*pre[j]<=maxxn;j++){
            primes[i*pre[j]]=1;
            if(i%pre[j]==0) break;
        }
    }
}
int main(){
    inint();   
    ll t;
    cin>>t;
    int time=0;
    while(t--){
        memset(mp,0,sizeof mp);
        int l,r;
        cin>>l>>r;
        for(ll i=1;i<=pos;i++){
            ll c=l/pre[i];
            if(l%pre[i]!=0) c++;
            for(ll j=c*pre[i];j<=r;j+=pre[i]){
                if(j==pre[i]) continue ;
                mp[j-l]=1;
            }
        }
        ll ans=0;
        for(ll i=l;i<=r;i++){
            if(!mp[i-l]) ans++;
        }
        if(l==1) ans--;
        printf("Case %d: %lld\n",++time,ans);
    }
    return 0;
}

 

标签:long,Help,ll,LightOJ,1197,素数,mp,1e6,primes
来源: https://www.cnblogs.com/Accepting/p/12601524.html