poj2689：素数筛

题目大意，给定l和u，求区间[l,u]内的素数中，相邻两个差最大和最小的素数其中 u的范围达到了2e9本质上需要找出n以内的所有素数，使用筛法。先保存50000(大于sqrt(2e9))内的所有素数，然后再去筛出区间[l,u]内的素数（题上给定l-u<=1000000）所以可以存下所有素数又由素数分布定理 素数个数s(n)~n/lnn并不是很大，所以找到所有素数保存在prime[]中扫一遍即可得到答案代码如下

#include<stdio.h>
#include<string.h>
#include <algorithm>
using namespace std;
#define N 100007
bool isprime0[50000];
int prime0[50000];
long long prime[100000];
bool isprime[1000010];
int num0;
int num;
long long x,y;
void setprime()
{
    num=0;
    for(int i=2;i<=50000;i++)
    {
        if(!isprime0[i])
            prime0[num0++]=i;
        for(int j=0;j<num0&&prime0[j]*i<=50000;j++)
        {
            isprime0[i*prime0[j]]=1;
            if(!(i%prime0[j]))
                break;
        }
    }
}
void setprime1()
{
    memset(isprime,0,sizeof(isprime));
    for(int i=0;i<num0;i++)
    {
        long long j=x/prime0[i];
        while(j*prime0[i]<x)
            j++;
        for(j=j*prime0[i];j<=y;j+=prime0[i])
            if(j/prime0[i]>1)
                isprime[j-x]=1;
    }
    if(x==1)
        isprime[0]=1;
    num=0;
    for(long long i=0;i<=y-x;i++)
    {
        if(!isprime[i])
            prime[num++]=x+i;
    }
}
int main()
{
    setprime();
    while(scanf("%I64d%I64d",&x,&y)!=EOF)
    {
        setprime1();
        long long a,b,c,d;
        long long mi=10000000,ma=0;
        if(num<2)
        {
            puts("There are no adjacent primes.");
            continue;
        }
        for(int i=0;i+1<num;i++)
        {
            long long p=prime[i+1]-prime[i];
            if(p<mi)
            {
                a=prime[i];
                b=prime[i+1];
                mi=p;
            }
            if(p>ma)
            {
                c=prime[i];
                d=prime[i+1];
                ma=p;
            }
        }
        printf("%I64d,%I64d are closest, %I64d,%I64d are most distant.\n",a,b,c,d);
    }
    return 0;
}