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;
}