POJ 2689 Prime Distance

做素数方面的题目，最重要的是筛选，切记不要试除，很容易超时。

这里用到二次筛选，并且区间长度为1e6这个信息点非常重要，可以在判断l～u之间的数的时候，符合的数减去l放进数组里，那样的话数组的大小也可以开，需要使用的时候拿出来加上l即可。

//
// main.cpp
// Richard
//
// Created by 邵金杰 on 16/8/2.
// Copyright © 2016年 邵金杰. All rights reserved.
//

#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=50000;
const int maxl=1000000;
int primes[maxn],isprime[maxn];
int notprime[maxl],prime[maxl];
int getprime()
{
int pi=0;
memset(isprime,1,sizeof(isprime));
for(int i=2;i<=maxn;i++)
{
if(isprime[i]) primes[pi++]=i;
for(int j=0;j<pi&&i*primes[j]<=maxn;j++)
{
isprime[i*primes[j]]=0;
if(i%primes[j]==0) break;
}
}
return pi;
}
void work(int pi,int l,int r)
{
memset(notprime,0,sizeof(notprime));
for(int i=0;i<pi&&(long long)primes[i]*primes[i]<=r;i++)
{
int s=l/primes[i]+(l%primes[i]>0);
if(s<2) s=2;
for(int j=s;(long long)j*primes[i]<=r;j++)
{
if((long long)j*primes[i]>=l) notprime[j*primes[i]-l]=1;
}
}
}
int main()
{
int pi=getprime();
int l,r;
while(scanf("%d%d",&l,&r)!=EOF)
{
int p=0,x1=0,x2=0,x3=0,x4=0,min,max;
work(pi,l,r);
for(int i=0;i<=r-l;i++) if(notprime[i]==0) prime[p++]=i+l;
max=-maxn,min=maxn;
for(int i=0;i<p-1;i++)
{
if(prime[i]==1) continue;
if(prime[i+1]-prime[i]<min) {min=prime[i+1]-prime[i];x1=prime[i];x2=prime[i+1];}
if(prime[i+1]-prime[i]>max) {max=prime[i+1]-prime[i];x3=prime[i];x4=prime[i+1];}
}
if(max==-maxn&&min==maxn) printf("There are no adjacent primes.\n");
else printf("%d,%d are closest, %d,%d are most distant.\n",x1,x2,x3,x4);
}
return 0;
}