poj2689Prime Distance 素数测试

朴素素数测试是O(x^1/2)的，每一个数都测试下来就炸了

然而如果全部预处理的话才是做大死，时间空间各种炸（大约有1亿个数）

所以怎么平衡一下呢？

其实在预处理的时候可以只处理一半：把2147483647^1/2内的质数全部预处理出来（这些就是要用的全部质数），然后用这些质数线性筛一筛就能得到正解

= =

没了？

没了。

还是要吐槽一下数论题目虽然代码、题解都很好写，但我不相信我能在赛场上想到正解。。。

#include <cstdio>
#include <iostream>
#define ma 65536 
using namespace std;
bool b[ma];
int a[ma];
bool f[1000001];
int main()
{
    int l,r,j=0;
    for(int i=2;i<=ma;i++)
    {
        if(!b[i])
            a[++j]=i;
        for(int k=1;k<=j && a[k]*i<=ma;k++)
        {
            b[a[k]*i]=true;
            if(i%a[k]==0)
                break;
        }
    }
    int n=j;
    while(~scanf("%d%d",&l,&r))
    {
        for(int i=0;i<=r-l;i++)
            f[i]=true;
        for(int i=1;i<=n;i++)
            for(int j=max(2,(l-1)/a[i]+1);j<=r/a[i];j++)
                f[a[i]*j-l]=false;
        int last=-1;
        int min1=-1,min2,max1=-1,max2;
        if(l==1)
            f[0]=false;
        for(int i=0;i<=r-l;i++)
        if(f[i])
        if(last==-1)
            last=i;
        else
        {
            if(min1==-1)
            {
                min1=last;
                min2=i;
            }
            if(i-last<min2-min1)
            {
                min1=last;
                min2=i;
            }
            if(max1==-1)
            {
                max1=last;
                max2=i;
            }
            if(i-last>max2-max1)
            {
                max1=last;
                max2=i;
            }
            last=i;
        }
        if(max1==-1)
            printf("There are no adjacent primes.
");
        else
            printf("%d,%d are closest, %d,%d are most distant.
",min1+l,min2+l,max1+l,max2+l);
    }
    return 0;
}

其实压一压行可以在30行左右搞定的