首页 > 代码库 > poj 2689 Prime Distance 【数论】【筛法求素数】
poj 2689 Prime Distance 【数论】【筛法求素数】
题目链接: 传送门
题目大意: 给你L和R两组数,L和R的范围是2^32,其间隔(即R-L最大为1,000,000.) 。让你求出L和R之间素数的最大间隔和最小的间隔。
比如 2 17。之间的最小素数间隔是2 3,最大的素数间隔是11 17。
要是直接进行一个2^32次方筛法然后在判断是会T的。
我们这样来想,筛法求素数的原理是什么:
/**vis数组标记为0则说明是素数*/ int vis[10005]; void getPrimevis(int n) { int m=sqrt(n+0.5); memset(vis,0,sizeof(vis)); for(int i=2; i<=m; i++) for(int j=i*i; j<=n; j+=i) vis[j]=1; }我们只用找到m=sqrt(n)内的素数和合数,然后对于1~n中的数如果输1~m中的素数的倍数,则筛掉。
则判断L到R中素数的个数就也不是很难了。
对于这个题目来说
先把m=sqrt(n)中的素数和合数都求出来,对于L~R中的数进行一个判断。比如让你求1000到2000的结果。
先把1000/c[i]的值记为s,在吧2000/c[i]的值记为t,然后对于s到t之间我们用相应的素数c[i]乘s 到 相应的素数c[i]乘t 得到的值减去1000(这样有助于减少数组的大小)。
这样就把L~R的数组转换成了r数组,下面就是一个暴力的搜的过程了。
主要是一些细节的处理,比如1的处理,边界情况的处理,自己向来这种题目写的蛋疼......
#include<queue> #include<map> #include<set> #include<vector> #include<stack> #include<ctime> #include<cstring> #include<iostream> #include<stdio.h> #include<cstdlib> #include<functional> #include<cmath> #define LL long long using namespace std; const int N=50000; int Max=0xfffffff; int a[50100]; int r[1000000],b[N+100],z=0; int main() { int a0,b0,i,j; for (i=2; i<=N; i++) if (!a[i]) { b[++z]=i; for (j=i*2; j<=N; j+=i) a[j]=1; } while (cin>>a0>>b0) { memset(r,0,sizeof(r)); int t=0,dis,mmax=-1,mmin=Max,m1,m2; for (i=1; i<=z; i++) { int s,t; s=(a0-1)/b[i]+1; t=b0/b[i]; for (j=s; j<=t; j++) if (j>1) r[j*b[i]-a0]=1; //利用1~50000中的素数,把合数筛掉,转换成r数组 } int k=-1; for (i=0; i<=b0-a0; i++) if (!r[i]) //若是质数 { if (k!=-1) { dis=i-k; if (dis>mmax) mmax=dis,m1=i+a0; if (dis<mmin) mmin=dis,m2=i+a0; } if (i+a0!=1) k=i; //跳过1 } if (mmax<0) cout<<"There are no adjacent primes."<<endl; else { cout<<m2-mmin<<','<<m2<<" are closest, "; cout<<m1-mmax<<','<<m1<<" are most distant."<<endl; } } return 0; }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。