poj2689两次筛法

我们都知道一次筛法求1-n之间的素数，这个筛法的算法复杂度为O(N)，但题目中的U,L最大值可为整型上限，用纯粹的暴力筛法肯定要超时？怎么办，用二次筛法。U和L之间的合数，质因子不超过O（L^0.5），于是用筛法选出50000内的素数即可，因为50000的平方大于整形上线了。再用这些素数去筛出U-L之间的合数，剩下的就是U-L之间的素数了，边筛边计算两个素数之间的差，就OK了。另外，尤其小心的是U为1的情况，做素数题往往1是个坑爹的东西。。。附代码：

#include <iostream>
#include <cmath>

using namespace std;

const int N=50000;
const int Max=0xfffffff;

int r[1000000],a[N+100],b[N+100],z;

int main()
{
    int a0,b0,i,j;
    for (i=2;i<=N;i++)
        if (!a[i])
        {
                  b[++z]=i;
                  for (j=i*2;j<=N;j+=i) a[j]=1;         
        }
    while (cin>>a0>>b0)
    {
          memset(r,0,sizeof(r));
          int t=0,dis,mmax=-1,mmin=Max,m1,m2;
          for (i=1;i<=z;i++)
          {