pku2689 素数查找

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<limits.h>
using namespace std;
__int64 prime[50000],next=0;
bool v[50000],res[10000100];
void getprime()//先找50000以内的素数
{
__int64 i,j;
memset(v,true,sizeof(v));
v[0]=v[1]=false;
for(i=2;i<50000;i++)
{
  if(v[i])
  {
   prime[next++]=i;
   if(i*i<50000)
   {
    for(j=i*i;j<50000;j+=i)
    v[j]=false;
   }
  }
}
}
int main()
{
__int64 A,B,i,j;
__int64 begin,num;
__int64 m1,m2,m3,m4,pre=-11,min,max;
getprime();//printf("%I64d %I64d\n",prime[0],prime[1]);
while(scanf("%I64d %I64d",&A,&B)!=EOF)
{
  memset(res,true,sizeof(res));
  if(A==1) A=2;//这里要特别注意，1不是素数，跳过1
  for(i=0;i<next&&prime[i]*prime[i]<=B;i++)
  {
   begin=A/prime[i]+(A%prime[i]>0);
   if(begin==1) begin++;
   for(j=begin*prime[i];j<=B;j+=prime[i])
    res[j-A]=false;
  }
  min=INT_MAX;max=0;
  for(i=A,num=0;i<=B;i++)
  {
   if(res[i-A])
   {
    num++;
    if(num!=1)
    {
     if(i-pre>max) {m1=i;m2=pre;max=i-pre;}
     if(i-pre<min) {m3=i;m4=pre;min=i-pre;}
    }
    pre=i;
   }
  }
  if(num<=1)
   printf("There are no adjacent primes.\n");
  else
  printf("%I64d,%I64d are closest, %I64d,%I64d are most distant.\n",m4,m3,m2,m1);
}
return 0;
}