POJ 2689 Prime Distance

题目：

两个数L,U，(1<=L< U<=2,147,483,647)

对于区间[L,U]内的所有相邻的素数，找到距离最小的那一对和距离最大的那一对

INPUT:

一行两个数L,U，相差不超过1,000,000.

OUTPUT:

按照示例输出最小相邻距离和最大相邻距离对应的一对数

2,3 are closest, 7,11 are most distant.

如果区间内的素数凑不成对，则输出

There are no adjacent primes.

分析：

由于[L,U]长度有限，类似hash的思路只用一个大小为1,000,000的数组存储，用筛法先找出50000以内的素数，再把[L,U]内的素数筛出来，然后遍历一遍，找出最小相邻距离和最大相邻距离对应的一对数

代码：

#include<cstring>
#include<iostream>
#include<cmath>
using namespace std;
typedef long long ll;
int prime[500010],a[10000010];
int main(){
	memset(prime,0,sizeof(prime));
	prime[0]=prime[1]=1; 
	for(int i=2;i<=50000;i++){
		if(prime[i]==0){
			for(int j=i*2;j<=50000;j+=i){
				prime[j]=1;
			}
		}
	}
	ll m,n;
	while(cin>>m>>n){
		memset(a,0,sizeof(a));
		for(int i=2;i<=50000;i++){
			if(prime[i]==0){
				ll k=(m/i)*i;
				if((m/i)<=1)k=2*i;
				for(;k<=n;k+=i){
					if(k>=m&&k<=n){
						a[k-m]=1;
					//	cout<<k<<' ';
					}
				}
			}
		}
	//	cout<<endl;
		if(m==1)a[0]=1;
		int flag=1;
		ll last=-1,cmin=n-m+1,cmax=-1,minlast=-1,minnext=-1,maxlaxt=-1,maxnext=-1;
		for(int i=0;i<=n-m;i++){
			if(flag==1&&a[i]==0){
				last=i;flag=0;
			}
			else if(a[i]==0){
				if(i-last<cmin){
					cmin=i-last;minlast=last;minnext=i;
				}
				if(i-last>cmax){
					cmax=i-last;maxlaxt=last;maxnext=i;
				}
				last=i;
			}
		}
		/*for(int i=0;i<=n-m;i++){
			cout<<a[i]<<' ';
		}
		cout<<endl;*/
		if(last==-1||cmin==n-m+1||cmax==-1){
			cout<<"There are no adjacent primes."<<endl;
		}
		else{
			cout<<minlast+m<<","<<minnext+m<<" are closest, "<<maxlaxt+m<<","<<maxnext+m<<" are most distant."<<endl;
		}
	}
	return 0;
}