POJ 2689/ ZOJ 1842: Prime Distance[两次筛法/筛选法平移]

这道题是学习素数筛法的经典，应用到了区间筛素数。具体思路是先筛出1到sqrt(2147483647)之间的所有素数，然后再通过已经晒好素数筛出给定区间的素数。

题目中的U,L最大值可为整型上限，用纯粹的暴力筛法肯定要超时？怎么办，用二次筛法。U和L之间的合数，质因子不超过O（L^0.5），于是用筛法选出50000内的素数即可，因为50000的平方大于整形上线了。再用这些素数去筛出U-L之间的合数，剩下的就是U-L之间的素数了，边筛边计算两个素数之间的差，就OK了。另外，尤其小心的是U为1的情况，做素数题往往1是个坑爹的东西。。。

唔，在网上对比了很多人的解法找到一个最简洁的，学习学习！真的代码与代码之间差别很大，有的让人看不下去，有的一目了然；有的繁杂冗余，有的简洁明了~比如这题有的人写了差不多100行，这里这种才40来行，真酷！

要提高自己啊。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
using namespace std;
const int N=50000;
int r[1000000],a[N+100],b[N+100],z;

int main()
{
int a0,b0,i,j;
for(i=2;i<=N;i++)//1st prime sieve
if(!a[i]){
b[++z]=i;
for(j=i*2;j<=N;j+=i) a[j]=1;
}
while(~scanf("%d%d",&a0,&b0)){
memset(r,0,sizeof(r));
int t=0,dis,max=-1,min=9999999,m1,m2;
for(i=1;i<=z;i++){//prime sieve again
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;
}

int k=-1;
for(i=0;i<=b0-a0;i++)
if(!r[i]){
if(k!=-1){
dis=i-k;
if(dis>max){
max=dis;
m1=i+a0;
}
if(dis<min){
min=dis;
m2=i+a0;
}
}
if(i+a0!=1) k=i;
}

if(max<0) printf("There are no adjacent primes.\n");
else printf("%d,%d are closest, %d,%d are most distant.\n",m2-min,m2,m1-max,m1);
}
return 0;
}