uva 10140 素数筛选（两次）

#include<iostream>
#include<cstring>
#include<cmath>
#include<cstdio>
using namespace std;

bool flag[50010];
int prime[8000];
bool r[1000010];
int num;
void Init()
{
    int i,j;
    num=0;
    memset(flag,true,sizeof(flag));
    flag[1]=flag[0]=0;
    for(i=2;i<50000;i++)
    {
        if(flag[i]) prime[num++]=i;
        for(j=0;j<num && i*prime[j]<50000;j++)
        {
            flag[i*prime[j]]=false;
            if(i%prime[j]==0) break;
        }
    }
}

int main()
{
    Init();
    long long a,b,up,mina,minb,maxa,maxb,min,max,i,j,k;
    while(cin>>a>>b)
    {
        if(a==1) a++;//a=1的情况
        min=20000000;
        max=0;
        memset(r,true,sizeof(r));
        for(i=0;i<num && prime[i]*prime[i]<=b;i++)
        {
            k=a/prime[i];
            if(k*prime[i]<a) k++;
            if(k<=1) k++;
            while(k*prime[i]<=b)
           {
                r[k*prime[i]-a]=false;
                k++;
           }
        }
        j=0;up=-1;
        for(i=0;i<=b-a;i++)
        {
            if(r[i])
            {
                if(up==-1) up=i;
                else
                {
                    if(min>i-up)
                    {
                        min=i-up;
                        mina=up;
                        minb=i;
                    }
                    if(max<i-up)
                    {
                        max=i-up;
                        maxa=up;
                        maxb=i;
                    }
                    up=i;
                }
            }
        }
        if(max==0) printf("There are no adjacent primes.\n");
        else printf("%lld,%lld are closest, %lld,%lld are most distant.\n",mina+a,minb+a,maxa+a,maxb+a);
    }
    return 0;
}