区间素数筛选 poj2689

题意：输出区间[a,b]中相邻素数差最大的和最小的   比如【2，17】素数有 2 3 5 7 11 13 17 输出最小2 3  最大13 17  

思路：套区间筛素数模板，遍历一次prime数组就好，注意模板里是【l，r)

模板原地址：http://www.cnblogs.com/nowandforever/p/4515612.html

模板的注释思路很清楚，慢慢学习

代码如下

#include <cstdio>#include <cstring>#include<iostream>#include <algorithm>using namespace std;typedef long long ll;const ll maxn = 1000005;bool is_prime[maxn];bool is_prime_small[maxn];ll prime[maxn];//区间内第I个素数 ll prime_num=0;//区间内素数个数 //对区间[a,b)内的整数执行筛法，is_prime[i-a]=true  ---  表示i是素数 注意这里下标偏移了a，所以从0开始。void segment_sieve(ll a,ll b) {if(a==1)a++; //1不是素数     for(ll i=0;i*i<b;++i) is_prime_small[i]=true; //对[2,sqrt(b))的初始化全为质数    for(ll i=0;i<b-a;++i) is_prime[i]=true; //对下标偏移后的[a,b)进行初始化    for(ll i=2;i*i<b;++i) {        if(is_prime_small[i]) {            for(ll j=2*i;j*j<b;j+=i) is_prime_small[j]=false;  //筛选[2,sqrt(b));            //(a+i-1)/i得到最接近a的i的倍数，最低是i的2倍，然后筛选            for(ll j=max(2LL,(a+i-1)/i)*i;j<b;j+=i) is_prime[j-a]=false;//筛【a,b】        }    }    for(ll i=0;i<b-a;++i)  //统计个数        if(is_prime[i]) prime[prime_num++]=i+a;}int main(){    ll a,b;    while(~scanf("%lld%lld",&a,&b))    {    b++;        prime_num=0;        memset(prime,0,sizeof(prime));        segment_sieve(a,b);        if(prime_num==0)        {        printf("There are no adjacent primes.\n");}else if(prime_num==1){printf("There are no adjacent primes.\n");}else{ll min=100000000;ll max=0;ll minl=0;ll minr=1;ll maxl=0;ll maxr=1;for(ll i=1;i<prime_num;i++){if(prime[i]-prime[i-1]>max){max=prime[i]-prime[i-1];maxl=i-1;maxr=i;}if(prime[i]-prime[i-1]<min){min=prime[i]-prime[i-1];minl=i-1;minr=i;}}printf("%I64d,%I64d are closest, %I64d,%I64d are most distant.\n",prime[minl],prime[minr],prime[maxl],prime[maxr]);}           }    return 0;}