模板 - 数学 - 数论 - 简单数论 - 区间素数筛

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
 
 
const int MAXN = 1e6+1e3;   //待筛的区间[L,R]长度
const int N = 50001;//保证大于(2^31-1)的算数平方根
bool prime[MAXN];
bool seive[N];
typedef long long ll;
int L,R,len;
 
void seg_seive(ll L,ll R)   //区间筛法
{
    len=R-L+1;
    for(int i=0;i<len;i++) prime[i]=1;
    if(1-L>=0) prime[1-L]=0;   //易错因为1不是素数也不是合数,这也是区间筛的一个易错bug
    for(ll i=2; i*i<=R ;i++)
    {
        if(seive[i])
        {
            for(ll j=max((ll)2,(L-1+i)/i)*i;j<=R;j+=i)  //第二个易错点,j必须从大于1,因为L可能小于i,但是seive[i]是素数。
                prime[j-L]=false;
        }
    }
}
int main()
{
 
    for(int i=2;i<N;i++) seive[i]=1;
    for(int i=2;i*i<N;i++)  //预处理
        if(seive[i])
            for(int j=2*i;j<N;j+=i)
                seive[j]=false;
        
    while(~scanf("%d%d",&L,&R))
    {
        seg_seive(L,R);
 
        int lmax,rmax,lmin,rmin;
        int mmax=-1,mmin=(1<<30);
        int t=-1;
 
        for(int i=0;i<len;i++)
            if(prime[i])
            {
                if(t>=0)
                {
                    if(mmax<i-t) mmax=i-t,lmax=t+L,rmax=i+L;
                    if(mmin>i-t) mmin=i-t,lmin=t+L,rmin=i+L;
                    t=i;
                }
                else t=i;
            }
 
        if(mmax>0) printf("%d,%d are closest, %d,%d are most distant.\n",lmin,rmin,lmax,rmax);
        else puts("There are no adjacent primes.");
    }
    return 0;
}

 

posted @ 2019-03-09 19:38  韵意  阅读(140)  评论(0编辑  收藏  举报