Prime Distance POJ - 2689 线性筛
Posted guangheli
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了Prime Distance POJ - 2689 线性筛相关的知识,希望对你有一定的参考价值。
一个数 $n$ 必有一个不超过 $sqrt n$ 的质因子。
打表处理出 $1$ 到 $sqrt n$ 的质因子后去筛掉属于 $L$ 到 $R$ 区间的素数即可。
Code:
#include<cstdio> #include<cstring> #include<algorithm> #include<iostream> using namespace std; const int Range=50000; const int N=1000000+233; int f[N],vis[Range+233],prime[Range]; int cnt; void get_prime(){ for(int i=2;i<=Range;++i){ if(!vis[i])prime[++cnt]=i; for(int j=1;j<=cnt&&prime[j]*i<=Range;++j){ vis[prime[j]*i]=1; if(i%prime[j]==0)break; } } } int main(){ get_prime(); int L,U; while(scanf("%d%d",&L,&U)!=EOF){ memset(f,0,sizeof(f)); if(L==1)L=2; for(int i=1;i<=cnt;++i){ int a=L%prime[i]==0?L/prime[i]:L/prime[i]+1; int b=U/prime[i]; for(int j=a;j<=b;++j)if(j>1)f[j*prime[i]-L]=1; } int p=-1,x1,x2,maxans=-1,minans=N,y1,y2; for(int i=0;i<=U-L;++i) if(f[i]==0){ if(p==-1){p=i;continue;}; if(maxans<i-p){maxans=i-p,x1=p+L,x2=i+L;} if(minans>i-p){minans=i-p,y1=p+L,y2=i+L;} p=i; } if(maxans==-1) cout<<"There are no adjacent primes."<<endl; else cout<<y1<<","<<y2<<" are closest, "<<x1<<","<<x2<<" are most distant."<<endl; } return 0; }
以上是关于Prime Distance POJ - 2689 线性筛的主要内容,如果未能解决你的问题,请参考以下文章
POJ-2689-Prime Distance(素数区间筛法)